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Introduction 



1.1 Particle physics, it is amazing! 

Progress in the understanding of elementary particles is amazing. For more than a cen- 
tury the smallest building blocks of nature have been studied, and discoveries are still 
being made today. While studying the ingredients of nature, fundamental and inspiring 
theories have been developed making this field of increasing interest. 

At the beginning of the last century, Einstein studied the concept of time in order to 
explain the discrepancy between the theories of Newton and Maxwell. This led to the 
publication of his theory of special relativity in 1905 |1]. Combining this theory and 
the quantum theory, Dirac predicted the existence of the antiparticle of both the electron 
and proton |2). The antiparticle of the electron, the positron, was discovered in 1933 by 
Anderson yl] and the antiproton was found in 1955 by Chamberlain et al.|4| 

Around the 1930's, explanations were sought for j8-decay, which is one particular 
form of nucleus-disintegration. Experimental studies of this phenomenon seemed to 
show that the energy before and after the decay were not the same: some energy was 
missing. To circumvent the potential violation of energy conservation (Newton's law), 
Pauli suggested between 1930 and 1933 at several conferences a new kind of particle', 
which would be produced during radioactive decay without notice. This new particle. 



'Pauli publicized this new particle at several conferences among which the Solvay Congresses in 1930 
and 1933. 
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called neutrino by Fermi, was observed in 1956 by Reines and Cowan |5j. 

Around 1960 particle accelerators discovered new kinds of hadronic matter. In 
order to classify the observed hadrons Gell-Mann in Ref. |6] and Zweig in Ref. J7| 
introduced, independently, a substructure with three types of quarks. Since then several 
other quark-types have been discovered and just a decade ago the last quark with a mass 
of almost 200 times the proton mass, the top quark, was discovered at Fermilab |8, 9]. 
Since this quark and its mass were already predicted on the basis of data taken by the 
large electron-positron collider at CERN, this was once again a stunning success for 
particle physics. 

In the last century particles have been found which were predicted by theory and 
theories have been developed on the basis of experimental observations. It is expected 
that during this century some of the predicted particles, such as the ones responsible 
for spontaneous symmetry breaking (the Higgs-sector), will be observed. The interplay 
between experiment and theory in this field is a guaranteed success to explore what 
nature will offer us next. 



1.2 QCD and single spin asymmetries 

Quantum chromodynamics (QCD) describes the interactions between quarks and glu- 
ons and is constructed from powerful theories and concepts. The main ingredients are: 
the theory of relativity, the quantum theory, and the concept of gauge invariance. The 
first gauge theory was developed more than a century ago. Around 1865 Maxwell 
wrote down his equations describing the interactions between electromagnetic fields 
and matter The equations are a set of differential equations which also raised some 
questions, one of which was that the potentials obeying the equations are not unique. 
This point became clarified towards the end of the nineteenth century; it was considered 
as a mathematical symmetry which was apparently left in the equations. This mathe- 
matical symmetry allowed for a set of transformations of the potentials which would 
not aff'ect physical observables. Nowadays this symmetry is named gauge invariance 
and the potentials are often called the gauge fields. 

In the beginning of the twentieth century gauge invariance was considered more 
seriously. While incorporating the symmetry in quantum mechanics, Fock discovered 
that, besides the gauge fields, the wave function of the electron should transform as well 
to maintain consistency with the theory of relativity. In order to preserve invariance 
of observables under gauge transformations, the wave function of the electron must 
obtain a space-time-dependent phase. However, the question remained whether the 
gauge fields were to be considered as fundamental fields or whether they just alleviated 
complex calculations. For a review on the historical roots of gauge invariance the reader 
is referred to Jackson, Okun lioll . 

In the second half of the twentieth century the question on whether potentials are 
more fundamental than electric and magnetic fields was finally addressed. Aharonov 
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Figure 1.1: The schematic setup of an Aharonov-Bohm like experiment. S represents 
a source of electrons producing an interference pattern on the screen owing to the two 
slits. The interference pattern shifts in a particular direction if the solenoid in B, point- 
ing out of the plane, is given a current. If the screen is far away from the two slits 
then the shift is proportional to ^dx ■ A(x). The path of the integral is the closed path 
formed by the two dashed lines and A(x) is the potential field. One can call this shift 
of interference pattern an asymmetry because the direction of the shift depends on the 
direction of the magnetic field. 



and Bohm apparently^ rediscovered that an electron can obtain a phase shift from its 
interaction with the potential even if it only travels in regions in which there is no 
electric or magnetic field flSJ. The experiment carried out by Chambers showed that 
instead of the electric and magnetic fields, the non-uniquely defined potentials should 
be considered as the fundamental fields in quantum electrodynamics 1 14]. A schematic 
setup of the experiment is given in Fig. II. II 

As compared to electrons and photons, the situation of quarks and gluons is much 
more involved. In contrast to electrons and photons, free quarks and gluons have for 
instance never been observed. They only seem to exist in a hadron (confinement) which 
indicates a strong interaction. On the other hand, perturbation theory turns out to pro- 
vide a satisfactory description for collisions involving hadrons at high energies, showing 
that the interaction at high energies must be weak. This particular scale dependence of 
the interaction strength confronted physicists with a big challenge. 

The solution came from a quantum field theory. In quantum field theories infinities 
often appear. In the 1940's Dyson, Feynman, Schwinger, and Tomonaga showed that in 
quantum electrodynamics such infinities can be handled by renormalizing the observ- 
ables. In contrast to the general expectation, 't Hooft and Veltman showed in 1972 that 
this was also possible for non-Abelian gauge theories llallTll . Using the machinery of 



^It seems that Ehrenberg and Siday already pointed out that enclosed magnetic fluxes could cause phase 
shifts. Their work 1,11.1 has been cited in the subsequent paper of Aharonov and Bohm [IZl . 
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Figure 1.2: An illustration of the interactions between the quark-jet and the target- 
remnant (Sivers effect). These interactions, which are on the amplitude level and lead 
to phases as we will see in chapter|3l give rise to interference contributions in the cross 
section and could produce single spin asymmetries in the idealized jet-production in 
semi-inclusive deep-inelastic scattering. 
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Figure 1.3: A nonzero transverse single target-spin asymmetry for n^ and n^ in 
electron-proton scattering measured by the HERMES collaboration. As will be ex- 
plained in chapter 13 the asymmetry is in the scaling limit proportional to the phase 
picked up by the outgoing quark and is also called the Sivers effect. The plot was taken 
fromRef.lllSll. 



1.2 QCD and single spin asymmetries 1 1 

't Hooft and Veltman, Gross, Politzer, and Wilczek derived in Ref. I18lll9ll that a scale- 
dependent interaction strength appears when quarks and gluons are characterized with 
a color The constructed non-Abelian gauge theory, called quantum chromodynamics, 
has a vanishing interaction strength at large momentum transfers - as those occurring in 
high energy collision experiments - which is called asymptotic freedom. 

Not being able to apply perturbation theory, most of the low energy regime of quan- 
tum chromodynamics is at present not calculable from first principles. Bjorken and 
Feynman introduced the idea of absorbing the nonperturbative part in parton probabil- 
ity functions which could be measured in several experiments. These functions describe 
how quarks are distributed in hadrons (distribution functions) or how they "decay" into 
a hadron and accompanying jet (fragmentation functions). The functions introduced by 
Feynman depend only on the longitudinal momentum fraction because at high energies 
the transverse momenta of quarks in hadrons can often be neglected. This somewhat 
ad hoc procedure of absorbing the nonperturbative part in functions, called the parton 
model, can be translated into more rigorous QCD and is very successful in describing 
various kinds of data. 

One observation, studied in this thesis, cannot be explained by the parton model, 
namely the observation of single spin asymmetries. In single spin asymmetries, one of 
the participating particles in a scattering process carries or acquires a certain polariza- 
tion. If the scattering cross section depends on the direction of this polarization, one 
has a single spin asymmetry. Large single spin asymmetries in inelastic collisions were 
discovered in hyperon-production in hadron-hadron scattering at Fermilab |201. Since 
then, single spin asymmetries have been observed in various processes. 

Several explanations for single spin asymmetries at large scales were developed 
over the last two decades. One of the most important ideas, proposed by Sivers in 
Ref. 1211 1221, was to allow for a nontrivial correlation between the transverse momen- 
tum of the quark and its polarization. After incorporating the transverse momenta of 
quarks in an extended version of the parton model, it is at present understood that there 
are two sources for single spin asymmetries. The first is the presence of interactions 
within a fragmentation process (see Collins |23]). The other source, appearing in the 
idealized single jet-production in semi-inclusive deep-inelastic scattering, see Fig. II. 21 
turns out to be the phase which the struck quark picks up due to its interaction with 
the target-remnant (also exists for fragmentation). This particular phase shows up as a 
gauge link, which is a mathematical operator, in the definition of transverse momentum 
dependent distribution functions. Since these functions are defined in terms of non- 
local operators inside matrix elements, the presence of this gauge link is also needed 
for invariance under local gauge transformations. Note that the obtained phase of the 
outgoing quark has similarities with the phases of the electrons in the Aharonov-Bohm 
experiment. Having the same origin of the effect for both asymmetries is surprising. 
The first nonzero experimental data which directly measures this phase was obtained 
by the HERMES collaboration in 2004 and is given in Fig. II. 31 
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1.3 Outline of the thesis 

The appearance and treatment of phases in several hard scattering processes will be 
studied in this thesis. In 2002 this particular topic became popular after Brodsky, 
Hwang, and Schmidt showed that such phases, leading to single spin asymmetries, 
could be generated within a model calculation |24). The obtained phases were at- 
tributed by Collins |25], Belitsky, Ji, and Yuan |26] to the presence of a fully closed 
gauge link in the definition of parton distribution functions. In this thesis these ideas 
are implemented in a diagrammatic expansion which is a field theoretical description 
of hard scattering processes. The effect of the gauge link is studied in several hard pro- 
cesses like hadron-hadron and lepton-hadron scattering. Although only QCD is studied 
in this thesis, gauge links also appear in other gauge theories like quantum electrody- 
namics. It is therefore to be expected that gauge links could provide a description of 
single spin asymmetries in pure electromagnetic scattering as well. 

For a full appreciation of the contents of this thesis familiarity with particle physics 
and quantum field theory is needed. Some excellent textbooks or reviews have been 
written by Anselmino, Efremov, Leader [27,]. Barone, Ratcliffe I28i1. Halzen, Mar- 
tin m. Leader ^, Peskin, Schroeder JmRvder IJl . WeinberglssllJl. 

Chapter 13 will begin with an introduction of the kinematics of some processes in 
which the hard scale is set by an electromagnetic interaction. A discussion of the di- 
agrammatic approach will be given together with the definitions of distribution and 
fragmentation functions. This chapter contains some new results from Ref. 1 35, 36]. 

In chapter 13 the diagrammatic approach will be applied to obtain cross sections 
of some electromagnetic processes assuming factorization. The gauge Unk inside the 
definition of parton distribution and fragmentation functions will be derived, showing 
the consistency of the applied approach at leading order in as (tree-level). The pres- 
ence of the gauge link will lead to the interesting prediction of Collins 112 511 that T-odd 
distribution functions in the Drell- Yan process appear with a different sign compared to 
semi-inclusive deep-inelastic scattering. This chapter is based on Ref. 135. 361 . 

Chapterl^will begin by considering gauge links in more complicated processes (and 
beyond tree-level). Besides the gauge links which are found in the electromagnetic pro- 
cesses new gauge links will be encountered which is quite a surprise. We will also see 
that the appearance of these new structures is an essential ingredient in the discussion of 
factorization. A set of tools will be developed which allows for a quick determination 
of the gauge link for arbitrary scattering processes. This chapter contains unpublished 
material of which some results were given in Ref. I37ll38ll39ll . 

In chapter|5lthe physical effect of the new gauge links will be illustrated in almost 
back-to-back hadron-production in hadron-hadron scattering. A new observable is con- 
structed which is directly sensitive to the intrinsic transverse momenta of partons. In the 
same way as T-odd distribution functions change sign in the electromagnetic processes, 
the T-odd distribution functions receive a gauge link dependent factor in the studied 
asymmetries of hadron-hadron scattering. This chapter is based on Ref. 14011 . 
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The formalism will be introduced which was initiated by Collins, Ralston, and Soper in 
Ref. |41, 42, 43, 44|. The formalism carefully considers the role of intrinsic transverse 
momenta of partons in hard scattering processes. Although the first part of the chapter 
is already present in the literature (and partly based on Ref. r45T), it remains worthwhile 
to look at some parts in more detail to elaborate upon the approximations and the phi- 
losophy behind certain approaches. Since this part is meant as an introduction for the 
non-experts, the reader can skip those sections which are familiar to him or her. 

Starting from section |231 the second part contains new elements which have been 
developed over the last few years. These new elements, which are one of the highlights 
in QCD-phenomenology, originate from the presence of Wilson lines or gauge links 
in parton distribution and fragmentation functions. These functions will be defined and 
they turn out to provide valuable information of partons inside hadrons and parton decay 
into hadrons. As we will see in the following chapters, the presence of these Wilson 
lines in the definitions of these functions lead to interesting predictions. A summary 
will be given at the end of this chapter. 
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2.1 Kinematics of electromagnetic scattering processes 

Physicists use often the Lorentz-invariance of the theory to choose the most conve- 
nient frame for their purposes. This has produced several frame definitions and frame- 
dependent interpretations. In principle the definitions and results can be compared by 
making the appropriate coordinate transformations but in practice this has often led to 
confusion. In this section several scattering processes will be introduced and their kine- 
matics will be set up such that theoretical predictions and experimental results can be 
compared frame-independently. 

As was advocated at the Transversity workshop in Trento 2004 (see also Ref. 14611 '). 
one can, in order to clarify this situation, express all results in easy-to-compare frame- 
independent observables which is possible owing to the Lorentz-invariance of the the- 
ory. For example, the variables in the invariant cross section are usually the momentum 
and spin vectors which have specific transformation properties. A much better choice 
would be to express the cross section in terms of the possible invariants. The invariants 
are frame-independent and can therefore be directly calculated in any frame! 

A drawback of this approach is that equations become rather lengthy and that we are 
not used to think in a frame-independent manner. To aid our intuition and to support 
the frames which are already in use, a Cartesian basis will be employed which will 
serve as an interface. Such bases were already introduced before, see for instance Lam, 
Tung 1 47] and Meng, Olness, Soper |48]. This will result in short expressions while 
maintaining manifest frame-independence. 

In this thesis the metric tensor of Bjorken and Drell 14911 will be employed, reading 

/° - -^" = -g"^ = -g'' = 1 , g'' = for / ^ ;, (2.1) 

and the antisymmetric tensor is normalized such that e"'^^ = 1. The Einstein summation 
convention will also be used, meaning that if a certain index appears twice in a product 
it is automatically summed over all its values unless stated otherwise. Furthermore, 
natural units with Ti - c - I will be used. 



2.1.1 Semi-inclusive deep-inelastic lepton-hadron scattering 

In the deep-inelastic scattering (DIS) process, a lepton with momentum / and mass m^, 
strikes with a large momentum difference (Z ■ P » M^) a hadron (sometimes called tar- 
get), with momentum P and mass M. The interaction, mediated through the exchange 
of a highly virtual photon with momentum q (and -M^ » q^), causes the hadron to 
break up into all kinds of particles being most often other hadrons. The measurement 
is called: inclusive if only the scattered electron is measured, semi-inclusive if an addi- 
tional particle (or more) with momentum Ph and mass Mh is measured, and exclusive if 
all (but one) final-state particles are detected. The semi-inclusive process is illustrated 
in Fig. l2.1h . 
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(a) 



(b) 



Figure 2.1: Semi-inclusive DIS. Figure (a) illustrates the process in leading order in 
(He.m.- Figure (b) presents the definition of the Cartesian basis, except for the vectors gy 
and Pt, all vectors lie in the plane shown. 



Traditionally the following Lorentz-invariants have been introduced to characterize 
experimental events 



Q' - -q\ 


W^ = (P + qf. 


s = {P + if. 


""'^IP-q 


^ PI' 


-2Ph ■ q 


PPh 
' Pq' 


(2.2) 




several Lorentz-invariants in semi-inclusive DIS 



where z = z/,(l + 0{M^IQf)) and y ^ Q^{1 + 0(M^/Q^))/(xbs). 

The Cartesian basis is defined through an orthogonal set of basis vectors e, (see 
Fig. 12.1b ). The space-like vector e, is chosen such that it is pointing in the opposite 
direction of q. The time-like vector e, is constructed from P subtracting its projection 
along q, and the transverse directions are fixed by choosing e ^ along the components of 
the sum of the lepton momenta which are perpendicular to e, and e,. The definition of 
e,, follows from demanding a right-handed coordinate system. This Cartesian basis is 
by construction frame independent and is mathematically defined as 



^^ 



-q" 



qi" + IxbP^" 



2' '"eV^T^' 

^l^f^-^,e] + ^,e\, A^H/7A,(foranyA), ^l ^ e'^^^^e.^e,^, 

Cartesian basis for semi-inclusive DIS 



II + 17 

|L+ii| 
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V / 






P,/ / 


>^^..^^Aaaaai 


'^ py\x 



(b) 



Figure 2.2: Two examples of fixed frames. Only the momentum P/, has a nonzero 
y-component and Ues in the dashed planes. 



Note that the antisymmetric tensor, ^, has been defined. 

In a Cartesian basis a general vector can be easily decomposed into a linear combi- 
nation of the basis vectors having frame-independent coeflicients, for example 



Ph - (Ph ■ e,) e, - (Ph ■ Cx) ex - (Pi, ■ e,.) e,. - (Ph ■ e.) e. 



(2.4) 



By using such decompositions head-on cross sections' can be written in terms of 
invariants only. Those cross sections can be indicated with the following notation: 
cr(Inv: p\, p2, • ■ •) = cr( all possible invariants of p\,p2, ■ ■ .). 

The advantage of this approach becomes clear by the following example. Suppose 
we are interested in azimuthal asymmetries in semi-inclusive DIS and we would like to 
predict or measure the quantity 



A= fd^P;,^ Pft^-e, 



dV(Inv:/,Z',P,5,n) 



(2.5) 



If our frame would be defined such that our proton is moving along the H-z-axis and our 
-Hjc-axis is proportional to /■'" + /'""" (see Fig. 12.2b ). then A would read in that frame 



frame a: A = I d^P/,^ |P/,^| 



sin (pv„. 



d-o-(lm:l,l',P,S,Ph) 



d^Pftx 



(2.6) 



while if our z-axis would lie in the opposite direction and keeping the same x-axis (see 
Fig. 12.2b ). A would read 



frame b: A = - fd^P,,^ |P,,^| 



sin^P,,^ 



d^o-(lm:l,l',P,S,Ph) 



d^Phi. 



(2.7) 



Although the appearance of the expressions diff'ers by a sign (due to a different def- 
inition of the azimuthal angle (frp,,^), the quantity A is the same in the two different 
coordinate systems. 

'Cross sections are generally not invariant under Lorentz transformations in contrast to head-on cross 
sections (see for instance Ref. |,31J). In the latter the initial particles are aligned. 
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(a) 



Figure 2.3: The Drell-Yan process. Figure (a) illustrates the process in leading order in 
«e.m.- Figure (b) presents the Cartesian basis in which only /, /', and Cy are out of the 
drawn plane. 



2.1.2 The Drell-Yan process 

In the Drell-Yan process two hadrons with momentum Pi and P2 collide and produce a 
virtual photon with a large invariant momentum squared q^ » M^. This virtual photon 
decays into an antilepton and lepton which are measured in the final state. The process 
has been illustrated in Fig. l2.3h . 

Characterizing invariants for this process are 



Q^ = q\ s = (Pi + P2f, XI = 



Q' 



X2 = 



Q' 



l-p, 



y = 



(2.8) 



2Pi-q' 2P2-q " qPi 

several Lorentz-invariants in Drell-Yan 



and the Cartesian basis is chosen to be (see also Fig. 12.3b ') 



^^ 



Q' 



^ = 



^P, 



(2.9) 



flV 



U V , U V 

e^e, + e^e^, 






^-iPl+P2)l 



^s-M\[l--^)-Ml{l-f) 
A^=g'[X(foranyA), 

Cartesian basis for Drell-Yan 



Since q is now time-like, e, is chosen along q. The vector e, is chosen perpendicular 
to q and such that Pi has a large positive component. The sum of the incoming hadron 
momenta perpendicular to e, and e, sets the direction of ex, and e,, follows from e j. 
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(a) 



(b) 



Figure 2.4: The electron-positron annihilation process. Figure (a) illustrates the process 
in leading order in ae.m.- Figure (b) presents the Cartesian basis in which only f ;,, and 
e,. are out of the plane. 



2.1.3 Semi-inclusive electron-positron annihilation 



In electron-positron annihilation an electron and positron collide with a large momen- 
tum difference, producing at leading order in as two jets (see Fig. l2.4h ). We will assume 
here that in both jets one hadron is detected and that Q^ > M^. 

The variables which characterize this process are (see Fig. 12.4b ') 



e^ = q\ s = (h + hf 



Zh\ = 



2Ph\ ■ q 



Zhl = 



2Ph2 ■ q 



(2.10) 



Q2 ' -- Ql ■ 

several Lorentz-invariants in electron-positron annihilation 



The Cartesian basis is defined through 






^ = 



Zl,2 hi 1 



e^/l-^ 



// = ^^ - ^,e] + 4e:, A'^ ^ //A, (for any A), e^" ^ e^'^^^'e.^e,,, 

'' " "" -l^e,,p^e^ey,e'^. (2.11) 



^n? 



^y — ^j_ ^xv> 



Cartesian basis for electron-positron annihilation 
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2.2 Cross sections 

In this section the cross section formula for semi-inclusive DIS will be derived and 
results for Drell-Yan and electron-positron annihilation will be stated. First, some con- 
ventions will be given. 

The helicity of a parton with momentum p and spin s is defined here to be 

A=^^. (2.12) 

Is -pi 

Dirac spinors and particle states are normalized such that 

u{k,A)u{k,A') — 2m6AA', (2.13) 

<P,^ I P',A') = 2Ej. (27t)^ 6\P' -P) 6m'- (2.14) 

The standard cross section for semi-inclusive DIS is (see for example Ref. 13111 ') 
1 d^Pi, dH' 



do- - — 



F(2nf2Ef,_ {2jif2Ex' 

^X { 7^^^^\M?i2nf6\l + P-Px-Ph-r)- (2.15) 

As we can see from this equation, the cross section is built up out of: several phase- 
space factors, a sum over all possible final states, an invariant amplitude, a delta- 
function which expresses momentum conservation, and a flux factor F which is given 
by 

F = AE\Ev\v,-vp\, (2.16) 

where v/ and vp are the velocities. The phase-space factors together with the delta- 
function are Lorentz invariant. Since the invariant amplitude is also Lorentz invariant 
the transformation properties of the cross section are set by the flux. In this thesis only 
head-on collisions will be considered, meaning that the motion of the initial particles is 
aligned. In that case the flux takes the form 

F ^2s{\+0[m-IQ^)). (2.17) 

For nucleons, consisting of strongly interacting quarks and gluons, the interaction 
between the electrons and the hadrons is at lowest order in a = e^ jiAn) mediated 
through the exchange of a virtual photon. Owing to the presence of charged quarks, 
the virtual photon feels the electromagnetic current between the incoming and outgo- 
ing hadrons. The process is illustrated in Fig. l2.5l and leads for the invariant amplitude 
to 

iM = {.-ie)u{l',A.')ypu{l,A) ^ ,APh,Px\(.-ie)FmP,S)u,,, + 0{e\ (2.18) 
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where the subscript c indicates that only connected matrix elements should be consid- 
ered. The blob in Fig. l2.5l expresses that all kinds of interactions are present. Math- 
ematically this means that the currents in Eg. 12. 181 are in the Heisenberg picture. A 
derivation of this equation is often omitted in textbooks but gives considerable insight 
into the approximations made. Therefore, a schematic derivation for the interested 
reader is provided in appendix l2.AI 



The square of the amplitude, needed for the cross section, 
can be written as a contraction between the leptonic tensor 
and a hadronic part in leading order of a, giving (neglecting 
lepton masses) 

|M|2 = ^L(;f)//f,;,(l+0(.2)), (2.19) 

4'f' = 6u' {2l/,+2lyi;-Q^gf,y+2iA e^yp^cff) , (2.20) 

-liv 
ilH) 




H'^L = ,n{P,s\no)\Px;Ph,s,,}out,c 



X ou.<Px; Ph,Sh\rmP,Sh,,„ (2.21) Figure 2.5: The virtual 

photon coupling to the 
where the A is the heUcity of the incoming electron and A' quark current in semi- 
is the helicity of the outgoing electron. Defining now the inclusive DIS. 
hadronic tensor W to be 



2^ ^(;h) - (2^ Z / iS^J'^''''^' -P.-P.-P-D K.^, (2.22) 
enables us to write the cross section as 



EnEv ^ttttt:^ = - — 4f Km (1 + Oia)) . 



d^l'd^P,, S 24 .V "(/„) 



(2.23) 
cross section for semi-inclusive DIS 



The interesting information on the distribution and fragmentation of quarks is captured 
in the hadronic tensor Wf^L. 

The inclusive cross section can be obtained by summing over all observed final- 
state hadrons and integrating over their phase space. This leads for the hadronic tensor 
to 



2M W'^ = — > I — 

^^^DK) 2;r4^j (2;r)32£p, 



(27r)'^d\l' +Px-P-l) 



X in{p, s\j''mPx)out,c out{Px\rmp, s )i„,c, (2.24) 



2 . 3 Operator product expansion 2 1 



and gives for the cross section 



Er^^-^-^fA"'^^o,s>(^^0(a)). (2.25) 

cross section for inclusive DIS 



The cross sections for Drell-Yan and electron-positron annihilation can be derived 
similarly. One obtains for Drell-Yan in terms of 

4^*^ = Saa' {2l/y + 21J'^ - Q%v + 2iAe^yp^qPr) , (2.26) 

H'^DY) = in{PA,SA;PB,SB\J''mPx}oul,c out{Px\rmPA, S a', Pb, S B}in,c, (2.27) 

j3 



<^^ = (5^ Z J ji^J^^^'^'^'^ ^P.-l-r- P.WZr. (2.28) 



the following cross section (a factor 2 was included for summing over lepton polariza- 
tions) 



^.^>' ^ = ^.^^''^^'m (1 + Oia)) . (2.29) 

cross section for Drell-Yan 



For electron-positron annihilation the cross section for producing two almost back-to- 
back hadrons reads (a factor i was included for averaging over lepton polarizations) 



dV 



a' 



,2 



£p,„ Ev,, ..^ . = -v^A ^-'Wr;,_, (1 + 0{a)) , (2.30) 



""- d^PM d^P/,2 42" "^ <«^-) 



cross section for electron-positron annihilation 



where (with |Q) representing the physical vacuum) 

4v '"' = 5.U' {2l^K + 2lyl'^ - Q^g^y + 2iA e^ypJPr) , (2.31) 

jflV 



HZ,-, = mnO)\Px;PhuSl,uPh2,S,,2)o 



X outiPx; Phu SkuPh2,Sh2\rm^\, (2.32) 
Wfl _, = ^ r "^f^ (27r)V(^ - Px - Pm - Ph2) Hf^ ,. (2.33) 

2.3 Operator product expansion 

There are two methods to gain more information from the hadronic tensor. In 1968 the 
first method was proposed by Wilson in Ref. L50.1 and is called the operator product 
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expansion. The second method, the diagrammatic expansion, was proposed by Politzer 
in Ref. pi'] in 1980 and will be introduced in the next section. 

The operator product expansion is useful for inclusive measurements. As an illus- 
tration let us consider the inclusive DIS process. Having no hadrons observed in the 
final state, the sum over all final QCD-states is complete. Together with the fact that 
the proton is a stable particle one can rewrite the hadronic tensor in Eq. l2.24l into^ 

2^<fl/s) = ^ / 'l'-^ ^"■'<^' ^ I ^■^'^^^' ■^''(O)J 1^' ^ >- ^2.34) 

According to the Einstein causality principle the commutator of two physical operators 
should vanish for space-like separations. In our case this means that only the area x^ > 
gives a contribution. Under the assumption that the hadronic tensor is well behaving 
for x^ > one can show that the main contribution comes from x^ » in the Bjorken 
limit (fixed xb and Q — > oo), implying light-cone dominance. 

The idea of Wilson, which was later proven in perturbation theory in 1970 by Zim- 
merman"^, is that for small separations one can make a Taylor expansion for operators. 
This expansion is called the operator product expansion and reads 

Oa(x)Ob(0) ^ Vq^aCO). (2.35) 

x^O ^ — ' 
n 

This above relation also holds for commutators and by using dispersion relations the 
short distance expansion can be applied for inclusive DIS. 

In Drell-Yan the sum over final QCD-states is complete which enables one to obtain 
a product of current operators. One can gain insight in the various structure functions 
in which the cross section can be decomposed but since these structure functions will 
depend on the two hadrons they are inconvenient for the study of the structure of a 
single nucleon. In addition, the process is not light-cone dominated which complicates 
the application of the operator product expansion. 

Another situation is encountered in electron-positron annihilation. When summing 
over all final states the commutator can be obtained, but if we are interested in how 
quarks decay into hadrons we would not be able to sum over such a complete set. The 
use of the operator product expansion is therefore limited here as well. 

Summarizing, the operator product expansion is a useful approach that finds appli- 
cations in inclusive DIS and electron-positron annihilation. At the same time, it is also 
limited to those processes. Parton distribution functions which are measured in DIS 
cannot be compared to more complex or less inclusive processes and quark decay can- 
not be studied within this approach. We will proceed by applying an extended form of 

^A necessary condition to have a stable particle is that the sum over energies of its possible decay prod- 
ucts is larger than the energy of the considered particle. Together with the fact that the zeroth momentum 
component of the virtual photon in DIS is positive, it can be shown that the hadronic tensor vanishes if q is 
replaced by -q. This enables one to obtain the commutator. 

''For reference, please consider chapter 20 of Ref. 1341 . 
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Feynman's parton model. In the original parton model it is assumed that the underlying 
process is a partonic scattering process multiplied by distribution and fragmentation 
functions. These functions describe the probability of finding on shell constituents in a 
hadron or how a quark decays into a particular hadron. In the next section the diagram- 
matic approach as an extension of this model will be discussed. 



2.4 The diagrammatic expansion and the parton model 

Background of the diagrammatic expansion 

The operator product expansion is of limited use for the study of the nucleon's structure. 
In the case of Drell-Yan we were not able to study the structure of a single nucleon 
although one could imagine that the chance of producing a virtual photon should just 
be proportional to the chances of finding a quark in the nucleon and an antiquark in the 
other nucleon. This idea of expressing the cross section in terms of probability functions 
which are then convoluted with some parton scattering cross section was suggested by 
Feynman and is nowadays called the parton model. 

The parton model had already lots of successes. It gave for instance an intuitive 
explanation for the approximate Bjorken scaling which was observed at SLAC. In an 
QCD-improved version of the parton model one could even predict the scaling vio- 
lation with a set of equations called evolution equations (for example see Altarelli, 
Parisi |52|). Another success of the parton model is the observation of jets. A jet is 
a set of particles of which their momentum differences can be characterized with a 
hadronic size. By assuming that these jets are produced by partons which "decay" into 
these jets one is able to predict the number of jets appearing in scattering processes. 
However, the appearance of jets also creates a problem with color Since partons carry 
color charges, they should somehow loose this color when decaying into a set of color- 
less hadrons. It appears that this issue does not influence the scattering cross sections 
at large momentum transfers. 

The success of the parton model relies on the asymptotic freedom property of 
QCD II8II19I1 . This property allows one to apply perturbation theory for elementary 
particle scattering in strong interaction physics in the presence of large scales. Strictly 
speaking, it remains, however, to be proven whether one can apply perturbation theory 
in hadronic scattering processes as well. The present idea is that suitable hadronic scat- 
tering processes can be described in terms of short-distance physics, the hard scattering 
part, and the long-distance nonperturbative physics which is captured in probability and 
decay functions. Since the latter are nonperturbative in nature, the approach should at 
least be self-consistent to all orders in perturbation theory. This description in separated 
terms is csWtA factorization. 

The diagrammatic approach is an extended form of Feynman's parton model. Orig- 
inating from field theory, the approach includes the possibility that several parton-fields 
from a hadron can participate in a scattering process (involving multi-parton correla- 
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tors), whereas the parton model only considers the possibility of hitting a single parton. 
The approach agrees with the operator product expansion when applicable. In 1980 
it was suggested by Politzer in Ref. 1 5 1 ] in order to describe the subleading orders in 
MjQ, involving "higher twist" operators in matrix elements (to be discussed in chap- 
ter|5}. Subsequently, it was applied by Ellis, Furmanski, and Petronzio in Ref. 153115411 . 
Although similar assumptions as in the parton model are made, the starting point is 
more general because it allows for more possible interactions. As we will see later, 
some of these interactions will provide an explanation for single spin asymmetries (see 
the work of Qiu and Sterman 15 5l l56l 15711 '). The diagrammatic approach was further 
developed and used in several applications, some of which to be discussed later in this 
thesis. 

The assumption in the diagrammatic approach is that in- 
teractions between the incoming hadrons and outgoing jets 
can be described in perturbation theory and hence can be di- 
agrammatically expanded with in the hard part a sufficiently 
small coupling constant. With respect to asymptotic free- 
dom this requires the incoming hadrons and outgoing jets 
to be well separated in momentum space. We will there- 
fore impose that the products of external momenta are large 
{Pi ■ Pj » M' for / ^ i) and assume that interactions between 
outgoing jets can be neglected. Non-perturbative physics in- 
side the jets and hadrons is maintained. Together with the 
assumption of adiabatically switching on and off the interac- 
tions, the applied assumptions are sufficient to describe gen- 
eral QCD-scattering processes. 

As an example, in the case of semi-inclusive DIS a hadron 
is detected which is well separated from the incoming nu- 
cleon in momentum space. Hence, there must have been a 
partonic scattering. The virtual photon has struck a quark 

which after interacting with the photon decays in a separate jet including the observed 
hadron. This process is illustrated in Fig. 12.61 There is also the chance that more jets 
are being produced which may not be observed. However, that possibility is expected 
to be subleading in a,. In section IZ21 an expression for the invariant amplitude was 
obtained in Eg. 12. 181 Making the assumption of adiabatically switching on and off the 
interactions more explicit (see also Fig. l2.6t . this result can be rewritten into (see also 
appendix l2.At 




Figure 2.6: The dia- 
grammatic approach il- 
lustrated for two jet- 
production and an ob- 
served hadron in semi- 
inclusive DIS. 



{e',Px,Ph\iT\P,e) 



u(k'^, A'){-ie)yuu(ke, A) — 

q 



: lim Id 



X lim ^A''ze-'''\,,{Ph,Px\U{U,,tf) 

tf—)O0 



U(tf,ti)(-ie)J';iz) 



U(ti,to)\PS}i„^„ (2.36) 
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where the subscript / denotes the interaction picture and to defines the quantization 
plane. 

The assumption of the diagrammatic expansion can now 
be used to expand the bracketed term in the above equa- 
tion. In general a complete expansion will connect an ar- 
bitrary number of lines to the several jets. Such a complete 
expansion, however, is not necessary. The interactions be- 
tween lines which are connected to one jet as illustrated in 
Fig. l2.7l can be absorbed in the matrix elements. Therefore, 
only those parts will be expanded which cannot be absorbed 
in one of the participating jets. In general this leads to ma- 
trix elements in which the interaction picture fields become 
Heisenberg fields. 




Figure 2.7: An interac- 
tions which can be ab- 
sorbed in the jet defini- 
tion 



Applying the diagrammatic expansion 

Using the diagrammatic approach the cross section for a general scattering process can 
be calculated in a number of steps. An outline of its derivation for two examples is 
given in aDDendix l2.BI 

1 . Write down all squared Feynman diagrams with an arbitrary amount of external 
parton-Hnes and connect them in all possible ways to the external jets and parti- 
cles. Each external parton-line carries an independent momentum variable (for 
example pi or kj). Any interaction which can be absorbed in one of the partici- 
pating jets should not be included. 



2. Replace the external spinors or polarization vectors of step 1 by an appropriate 
correlator as will be defined below. For instance, uu — > <£), Me'^S — » <E)^, etc. 

3. Integrate over all parton momenta and impose total momentum conservation by 
adding (2;7r)'*5'*(incoming - outgoing parton momenta). 



4. If there is any QED part in the diagram, calculate that part with ordinary Feynman 
rules. 



Divide by the flux factor and multiply by the phase-space factors of the produced 
particles, d^k; /((27r)^2£'^). 
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Some of the correlators which appear in cross sections are (see also Fig. l2.8l and Fig. l2.9> 



'^aabip)^ [^ e'^HP, S |^,-,,(0) if^aOlP, S >e, (2.37) 

_L_^ e'Pie'P'^''-^\P, S |^,;,(0) Af (/,) ^,;,(^)|P, S )„ (2.38) 



•^:^;;::X.>''— '^") 



/ 



d g d 7/1 . . . d 7;„ ,x , •;„(;„ -f) ,>„(;;„ -f) 
(2;r)4(«+i) 

X <p,5i^,;,(0) a;;'(77i) . . .a;'"(77„) ^,,„(^)|p,5>e, (2.39) 



sot7ie parton distribution correlators containing two quark-fields 



tin n j4> 



X out{Ph,Si,;Px\4'j,bm^)c, (2.40) 

X out<f /„ 5/,; Pxl^,>(0)Af (77)|Q)c, (2.41) 



^vX„./^'^"-'^«> 



y f dPx rd^d_77i_^^^^,,^^_,,,,^^^^^^_,,_,,__ 



■^J (27r)32£p,J (2;r)4(«+i) 

x<Q|i^;,„(^)m,5^;Px>out,c 

X ,u,{Ph,SH;Px\ikj,b(o)A';;im) ■ ■ .A;'"(/7„)iQ>e. (2.42) 



sot7ie parton fragmentation correlators containing two quark-fields 



All fields, so iff and A, carry a color index (a, fe e 1, 2, 3, Z, e 1, . . . , 8) over which it is 
summed in the cross section (no averaging for initial states), and the indices {/, j} denote 
Dirac indices. Time-ordering is not present because the matrix elements are connected 
and the correlators will be integrated over the small parton's momentum component, 
putting the fields on the light-front (see Jaff'e 1.5 8.. .5 9.1 . Diehl, Gousset L60.I '). 
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Figure 2.8: Some of the distribution correlators containing two quark-fields 




Ph^.. -.^Ph 




Figure 2.9: Some of the fragmentation correlators containing two quark-fields 

To illustrate these rules we return to our example in which there are only two jets 
produced (see the amplitude diagram in Fig. l2.6> . In step 1 we have to write down for 
the cross section all squared Feynman diagrams which contribute to this process. Let us 
consider one of the contributing diagrams as given in Fig. 12.1 Oh . Expressing the colors 
very explicitly by giving the spinors and polarization vector a color charge (superscript) 
and a color index (subscript), one finds for the diagram (note also that m^' ~ dq^a) 



Yu'f{k)iyjj 



?1,?2, 



(k-p\Y-m?-+ie 



-yuf\p - p,)gef'^yp,) X -d'?'\k)y6,,uf\p) 



= Tr' 



'^•'^\\uc{k)ua{k)\iya 



Li 



{k—pip--nP-+ie 



J [Mfc' 



(p-Pi)ge-f{pi)uc{p)]r\ (2.43) 



where Tr^'' stands for a trace in color and Dirac space, and Ucik) = Y^ u^ and similarly 
for e;. According to step 2 we replace the bracketed terms by A and O^ and in step 3 we 
integrate over p, k, and pi, and multiply this by {27T)*S{p + q-k). In step 4 we multiply 
the result of step 3 with the leptonic tensor and the photon propagator Applying step 5 
this diagram contributes to the cross section as 



da- = 



d^Pft d^l' "4 



jjlH) 

2s (2nf2Ev, (2ny2Ex ''" 



X Tr^'"^ [A(fc)] iyjii 



4 n 

^ J dV d'*^ dVi (27r)V(p + q-k) 



{k-p\y-—m?-+ie 



/[i'X(/^,Pi)]/ + 



(2.44) 



where the dots denote contributions from other diagrams. 

The previous example illustrates that the diagrammatic approach enables one to 
express any cross section in a set of correlators. The result is in general an infinite 
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outgoing hadron 
and its parent jet 

k 



P-Pi 




incoming proton 
and its outgoing remnant 

(a) 



Ph// V'^ 




(b) 



Figure 2.10: Application of the diagrammatic approach in semi-inclusive DIS. Fig- 
ure (a) represents the two jets and the interactions between them, and in figure (b) the 
correlators describing the jets are included. 



sum of partonic scattering diagrams connected to all kinds of correlators. That result 
is fairly exact; it relies on the possibility of applying perturbation theory between the 
external jets and the possibility of defining the jets. However, the cross section is not 
yet a product of probability functions as it is in the parton model. In order to obtain this 
product similar assumptions as in the parton model will be made here. 

These assumptions are that the parton-lines connecting the soft blobs are approx- 
imately on their mass-shell and collinear with their parent hadron, the hadron which 
is connected to the blob of the parton. This assumption is less strict than the assump- 
tions of the successful parton model in which partons were treated as in essence free, 
collinear, and on the mass-shell. Since the assumptions made here are similar to those 
made in the parton model, we will refer to them as the parton model assumptions. 

To exploit these assumptions a set of light-like vectors will be chosen such that 
P ~ Qn++ 0(M^/Q) «_ and Ph~ Qn_+ 0(M^/Q) n+ and 



«_■«+ = 1, 



and 



(2.45) 



where the bar on a vector denotes reversal of the spatial vector components. This intro- 
duction of vectors, called a Sudakov-decomposition, has the advantage that it seems that 
our target-hadron is now moving very fast in a particular direction without actually hav- 
ing boosted the target. The frame in which the proton is moving very fast is also called 
the infinite momentum frame and if its momentum is opposite to the momentum of the 
photon one also refers to it as the Breit frame. Using the parton model assumptions for 
the integration over p ■ «_ and k ■ n+ and applying the relation Aah{k) = 6ah Tr*' A{k)/3, 
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one obtains for Eq. l2.44l fp^ = p ■ n^, etc.) 



xTr" 



i Tr- ( r d^- A(^)) Tr- (iyjA ^ " f ' "^ " . |,,.or 



X Jdp-dprcDX(A/^i))/j ^^__^^{\+0{m^IQ^)) + ..., (2.46) 

k-=q- 

where the subscript T denotes transverse components with respect to n_ and n+. 

The applied Sudakov-decomposition can be used for general scattering processes 
as long as the scalar product of the observed momenta is large. For each observed 
hadron one can introduce a light-like vector along which the hadron is moving. Since 
the parton-lines are approximately collinear and on shell, one of the components of the 
parton momenta, p ■ n, must appear to be very small. In general one should be able 
to neglect these components in the hard scattering part such that one can integrate the 
considered correlator over this variable. In the next section we will see that correlators 
which are integrated over the small momentum components are probabilities in leading 
order in M/g- 



2.5 Quark distribution functions for spin-i hadrons 

The various functions for spin-| hadrons will be introduced and their relevance will be 
pointed out. For spin-1 targets the reader is referred to Bacchetta, Mulders |61]. To 
define the parton distributions a set of light-like vectors is constructed such that 



1 = n_ • n^. 




ri- ~ n+, 


P-2P-"-^^'"'' 




4' ^ ef"^^^n,,n^^. 


//^g^^-z/X-nX, 




A% = g^^Av, for any A. (2.47) 


the basis in 


which 


parton distribution functions are defined 



For any vector A we also define A^ = A ■ n^, which means that P^ is defined to be a 
Lorentz-invariant. To describe the spin of the hadron one usually introduces 

5 =-5z.^«-+5z.^n++5r, with 5^ + S^j. = !> (2-48) 

which satisfies the necessary constraints: PS = 0, 5" = if P = 0, and 5^ = -1. 
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In the next subsections we will parametrize an expansion in M/P^ of quark-quark 
correlators, but M/P^ does not have to be small. However, in order to make use of 
the truncated expansion, calculations in the next chapters will be performed such that 
P^ » M. Since the light-like vectors «_ and n+ are defined up to a rescaling {n+^>an+, 
ri-^a^hi-), this does not put any constraint on P or the frame. In fact, for a target at 
rest one has for instance n+ = (M/2P+) (1,0, 0, 1) and n_ = (P+/M) (1,0, 0, -1) where 
P^ :s> M can still be chosen. 

As discussed in the previous section, one encounters in the diagrammatic expan- 
sion an infinite set of correlators which can all be integrated over the small momentum 
components. As we will see in the next chapter, this infinite set can be rewritten into a 
single new correlator containing the gauge link (to be defined below). In the discussed 
electromagnetic processes basically two kinds of correlators appear in the final result. 

The first kind appears in cross sections which are not sensitive to the transverse 
momenta of the constituents and is the so-called integrated correlator. Including a 
Wilson line operator £,, it reads 

<l),v(x,P,5) = f^-e'''''^'{PS\ifrj(0)£''-^\0-,n^i(^)\P,S}c\e=o^ (2.49) 

where over the color indices was summed and where x is the longitudinal momentum 
fraction of the quark with respect to its parent hadron, x = p^ jP^ . The Wilson line 
operator, or also called gauge link, is a 3 x 3 color-matrix-operator and makes the 
bilocal operator ij/ibiO) ^I'UaiO invariant under color gauge transformations. A Wilson 
line along a path "BFiA) with H'^(O) - a^ and H''(l) - // is defined as 

1 
£(a,b) =l-igJdA ^A,(E(A)) 



+ i-igfJdAi ^A,(E(Ai)) JdA2 ^A,(E(A2)) + ..., (2.50) 

^1 

where A^ = A^f'. In Eg. 12. 491 an abbreviation was introduced for links along straight 
paths. In the abbreviation it is indicated which variables are constant along the path (O7- 
and ^^) and which coordinates are running (the minus components). The path for this 
case is illustrated in Fig. 12.11b . The integrated correlator will be parametrized in the 
next subsection. 

The second kind of correlator is encountered in cross sections which are sensitive 
to the transverse momenta of the quarks. Calling it the unintegrated correlator, it also 
contains a Wilson line operator and reads 



(I)M(^, pT,P,S)^ J "^^^^^^f e'P^P, 5 I0-/O) £[^1(0, n H^m S )c 



,..0 ' (2.51) 
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Figure 2.1 1; The paths of the various gauge Unks which connect the two quark-fields 
in the correlator. 



where 



£^(0,^) = i:""' ^*(0", +co-)i:=^°°" ^\Qt,^t)-C^^' ^'(+co-,d. 



(2.52) 



and similarly for {^^,^^, oo") «-> {^-,^+, oo+). In these correlators the link runs via a 
nontrivial path like the ones drawn in Fig. 12. lib and Fig. 12. lib . After an integration 
over pt the unintegrated correlator reduces to the integrated correlator. The nontrivial 
paths in the unintegrated correlator are a source for interesting phenomena like single 
spin asymmetries. These correlators will be parametrized in the second subsection. 

As a last remark, in chapter|2and|5]we will also encounter a gauge link of which 
its path is closed, see Fig. 12. lib . This "closed" gauge link or Wilson loop is defined as 



2.5.1 Integrated distribution functions 



(2.53) 



The specific form of the integrated distribution correlator (Eg . 12.491 has the following 
analytical properties 

cD'(x,P,5) = /<D(ji:,P,5)/ 

<^*(x,P,S) = iysC <D(x,/',5) iy^C 

where C denotes charge conjugation and 75 = iy^y^y^y^. Note that P ^ P implies 
that «_ «-> «+. Using these constraints the correlator has been decomposed on a basis 



(hermiticity). 


(2.54) 


(parity). 


(2.55) 


(time-reversal). 


(2.56) 
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of Dirac structures i44il62L 16311 . In the notation of Ref. 16311 this gives'* (for conventions 
on names see for instance Ref. I64ll ') 



<D(x, P, 5) = -(/i ^+ + Sl 81 75?^+ + hi 75^ TVi+) (twist 2) 

+ j^\e + gT75^T+SLhiy5 — - — I (twist 3) 

+ higher twist. (2.57) 

parametrization of the integrated correlator 



where all functions depend on x and on a renormalization scale. The scale dependence 
can be calculated by applying evolution equations (for example see Ref. I65l 16611 \ 
The above expansion in MjP^ is not to all orders. Calculations of cross sections 
should therefore be constructed such that MjP^ is small in order to employ the above 
parametrization. 

In the parametrization we have indicated the twist of the distribution functions. The 
twist of a distribution function defines at which order in MjP^ the function appears in 
the parametrization. This definition of "operational twist" for nonlocal operators was 
introduced by Jaffe in Ref. |67]. The standard definition of twist, which counts the 
dimensions of local operators, agrees with this definition for local operators appearing 
in the Taylor expansion of the nonlocal matrix element. 

As will be shown in the next subsection, the leading twist functions allow for a 
probability interpretation. In order to obtain a parton interpretation for the functions, 
the Heisenberg field operators, which are present in the correlators, can be expanded 
in creation and annihilation operators at particular point in time. Therefore, the func- 
tions which are integrated over x and are local in time allow for a parton interpretation. 
Moreover, in those matrix elements the gauge link has vanished, simplifying the inter- 
pretation. The function f\ describes the chance of hitting an unpolarized quark (red, 
green, or blue) in an unpolarized nucleon, the function g\ describes longitudinally po- 
larized quarks - or actually chirally left or right handed quarks - in a longitudinally 
polarized nucleon, and hi or transversity describes it for transverse polarizations. This 
last function is at present the only leading twist integrated function which has not yet 
been measured. Recently, it has been suggested in Ref. 1681 by Bakker, Leader, and 
Trueman that this function appears together with orbital angular momentum of the 
partons in a simple sum-rule. This would provide for the first time in semi-inclusive 
DIS access to the orbital angular momentum of partons inside a nucleon, making the 
transversity function of increasing interest. For more information on the interpretation 
of these functions the reader is referred to Ref. I28ll69ll. 



^T-odd integrated distribution functions are in this approach zero and have therefore been discarded. 
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2.5.2 Transverse momentum dependent distribution functions 

Gauge Invariant correlators and T-odd behavior 

In electromagnetic processes one encounters at first sight two kinds of unintegrated 
correlators in the final result of the diagrammatic expansion. They are defined as 



0M(x, pT,P,S)^ J ^^2^''"'^^^' ^ '"^^'^^^ ^'^'^°' ^^ ^'^^^'^' ^ ^' 



,..0 ' (2.58) 



where the two different paths of the gauge links X'"^^ are indicated in Fig. l2.1lb .c. In 
general also other paths will appear in link operators, but those will be discussed in 
chapter|3]and|5] The analytical structure of the correlator has the following properties 

<t>^^^\x,pT,P,S) = y° <:>^^\x,pT,P,S)y° (hermiticity), (2.59) 

^^^\x,pT,P,S) = y° <:>^^\x,-pT,P,-S) / (parity), (2.60) 

<I)M'(x,/57,f,5) = /ysC ^^''\x,-pT,P,S) ijsC (time-reversal). (2.61) 

The last equation shows that the time-reversal operation relates the two different paths 
of the gauge links. This enables us to decompose O in two classes, one which is the 
average and is called T-even, and one which is the difference and is called T-odd (T-odd 
does not mean breaking of time-reversal in QCD, the name is similar to P-odd) 

<l)t'^-'=^™l(x,/?7-,P,5)= i(oW(x,7?7-,i',5)-H<l)H(x,/?7-,i',5)), (2.62) 

^^'^'°'^'^\x,pT,P,S)^ ]-{^'''^^{x,pT,P,S)-(b''-\x,pT,P,S)). (2.63) 

Although the phenomenology of T-odd distribution functions was studied, see for 
instance Anselmino, D'Alesio, Boglione, Murgia 17(1 1? iL 17211 and Boer, Mulders 17311 . 
T-odd distribution functions were not really believed to exist as separate distributions 
(for example see Collins in Ref. L23.] '). After Brodsky, Hwang, and Schmidt showed 
in Ref. I24ll that unsuppressed T-odd effects could be generated in a particular model, 
the existence of T-odd distribution functions was taken as a serious possibility (see also 
Colhns 1I25I] and Behtsky, Ji, Yuan |26|). At present the data of HERMES as shown in 
chapter[n Fig. II. 31 indicate that T-odd functions might really exist. 

Another mechanism for T-odd functions has also been suggested by Anselmino, 
Barone, Drago, and Murgia in Ref. ^4J based on nonstandard time-reversal. This 
mechanism is not taken into account here but could be included at a later stage. If 
this mechanism is realized in nature it would lead to universality problems for the dis- 
tribution functions similar to those akeady appearing for fragmentation functions. The 
latter problem will be discussed in the next section. 

The T-even and T-odd parts of O'^^ have identical parity and hermiticity proper- 
ties and they can be parametrized in a set of functions. Before going to the explicit 
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parametrizations it is interesting to note that a distribution correlator can now be given 
by a T-even and a sign dependent T-odd part 



This means that T-odd distribution functions enter with a sign 125L I75r depending on 
the path of the gauge link. In the next chapters we will see that this path is set by the 
process or subprocess. 

The introduction of the distribution functions 

The correlator can be parametrized as follows I3a.l63ll76ll77ll78ll79ll 
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(twist 2) 

(twist 2) 

(twist 3, unpolarized) 

(twist 3, polarized) 

(twist 3, polarized) 

(twist 3, polarized) 

(twist 2) 

(twist 3, unpolarized) 
(twist 3, polarized) 
(twist 3, polarized) 
(twist 3, polarized) 



(2.65) 
parametrization of the distribution quark-quark correlator 



All functions have the arguments x and p^ and also depend on a renormalization scale. 
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In contrast to the integrated distribution functions the scale dependence is not known 
for transverse momentum dependent functions (see for instance Henneman f6CT]). 

In the T-odd correlator the new functions g-'-, Cj, and f^ are included. The function 
g-^ (as defined in Ref. 13611 ') and the existence of the others were discovered in Ref. 13611 . 
Subsequently, a complete parametrization was given by Goeke, Metz, and Schlegel in 
Ref. |79]. The fact that ^Spr <I>'^"'**(x,/?7-) = leads to constraints for the T-odd 
functions h, fj, e^, and fj (see for example Ref. l35ll79ll V 

The first transverse moment of the correlator Q> and some function fi is defined as 



%{x, P,S)= f d^PT Pt 'I>(^, Pt,P,S). 



(2.66) 
(2.67) 



The introduced functions describe how the quarks are distributed in the nucleon. 
The leading twist functions (twist 2) are again probability functions and therefore con- 
tain valuable information. For instance, the functions fi(x,pj) and giL(x,Pj) are gen- 
eralizations of /i(x) and gi{x). For more information on the interpretation of T-even 
functions the reader is referred to Ref. |64, 80|. 

For the interested reader the proof for the probability interpretation of /i is given 
here (see for instance Ref. 1159118 ill , and Ref. ll82ll for related work) 



/i = iTr^_(D 



[T-even] 



V2 '? =0 



lf+=0 



>0. 



(2.68) 



For the other leading twist functions (T-even and T-odd) the proof is analogous. The 
fact that the leading twist functions are probabilities leads to several positivity bounds, 
see for instance Soffer [831, and Bacchetta et al. f81'|. 

The leading twist T-odd functions (such as the Sivers function f^j) do not have 
a parton interpretation in terms of quarks although they are probabilities as well. To 
obtain an interpretation we consider the transverse moment of a T-odd correlator (see 
also Ea. l2.66> and rewrite the quark transverse momentum as a derivative acting on the 
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gauge links alone 

X [idi [£}^\o, n--c^'^(o, n]) u^m s }c\e=o ■ (2.69) 

p*=.xP* 

Using identities of Ref. |3S] (the second identity can be proven by using that G^" ~ 
[/D+, iD"l id^£^*-^'-(T]-,^-) = X.^"-^Hi]-,^-)iD+(^), and shifting /D+ to the side) 



id-^j:^---"\0T,^T) 




= X^"""' ''\QT.^T)iD"j{±^,a\^T), 




^■z)^(^,fl^^r)£^^•"'(^,^) 








a\^T)Z^^'''*{T]-,n^ (2.70) 


r 


identities concerning gauge links 



where a^ is some constant and iDj{^ ,a^,^T) = idZ + gAji^ ,a^ ,^t), the derivative 
on the difference of the links can be written as 

<* 2J 27r 

xX°^'f'(0-,77-)Gr('7)X°^'^'(?7^r)'AK^)l^,5)c|„-r=o- (2.71) 

This matrix element is also called a gluonic pole matrix element. In 1991 such a matrix 
element was suggested by Qiu and Sterman in Ref. 1.84. .85.1 as an explanation for single 
spin asymmetries in photon-production in hadron-hadron scattering and it has been 
subsequently studied in several other articles |57, 86, 87, 88, 89, 90|. We see here that 
the same matrix element appears in the two separately suggested mechanisms for T-odd 
effects 13511 (the soft gluon effects of Qiu and Sterman and the gauge link). 

At first sight it seems strange that the nontrivial path is the origin of T-odd functions. 
One might think that such functions should vanish, since in the light-cone gauge, A^ = 
0, the T-odd distribution functions become proportional to transverse gauge links at 
infinity (see Eg. 12. 691 and Fig. 12. lib , d) and shouldn't those matrix elements vanish? 
Since we are at present not able to calculate those matrix elements in QCD we do not 
really know the answer, but we do know that nonzero A -fields can give effects in areas 
where the physical fields, like G"^, are required to vanish. An example of such an effect 
is the Aharanov-Bohm experiment which was discussed in chapter^ 
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The Lorentz invariance relations and g^ 

Before the paper of Brodsky, Hwang, and Schmidt 12411 appeared, physical effects from 
the gauge Hnk were assumed to be absent which led to several interesting observations. 
Not only did T-odd effects disappear, but relations between the various functions in 
the correlator were also obtained. By arguing that the correlator <^(x,pT,P,S) could 
be written in terms of only fermion fields, the starting point then was another object 
<l>(p, P,S) which is defined by 

Oijip, P,S)^ f ^ e'P^ {P, S I iA/0) «A,(f ) 1^, S )c. (2.72) 

J (27r)^ 

This quantity <l>(p, P,S) is the fully unintegrated correlator, from which (^(x, pt,P,S) 
(without gauge link) is obtained via 



/ 



<D(x,/?7-,P,5)= dp- ^(p,P,S). (2.73) 

Without the gauge link the parametrization of the object ^{p, P,S) contains less func- 
tions than 0{x,pt,P,S). This led to the so-called Lorentz-invariance relations^ (see 
Boer, Jakob, Henneman, Mulders, Tangerman 163. 78. 94. 951 ) 

grix) = gx(x) + —g\j\x), g^{x) = -—gf^\x), hjix) = -—h^^^\x), 

hdx) - /zi(x) - ^hfl\x), h\jp\ p2 ) = hT(p\ Vh-) - hjiP^ Pr). (2.74) 

Similar relations for T-odd distribution functions were obtained as well. 

Taking the gauge link into account affects these relations. Since the gauge link 
runs in the n_-direction via infinity (see Fig. l2.1lL it is not clear how to construct an 
n_-independent '^(p,P,S) which after an integration over p^ leads to 0(x, pT,P,S) 
containing the gauge link. In 2003, this was made explicit by Goeke, Metz, Pobylitsa, 
and Polyakov in Ref. ll96ll . Taking the n_-dependence into account, they showed that 
the former proof of the Lorentz-invariance relations failed. 

It seems to be impossible to maintain the Lorentz-invariance relations for T-odd 
functions since the involved matrix elements are intrinsically nonlocal (see for instance 
Ea. l2.71> . This is in contrast to the first transverse moment of a T-even function which 
is local in the light-cone gauge. This could imply that the Lorentz-invariance relations 
might still hold for the T-even functions. 

The n_-dependence of the gauge link implies not only the need to revisit the Lorentz- 
invariance relations [96], but also leads to new functions in the parametrizations as dis- 
covered in Ref. |36] and confirmed in Ref. 17911 . Since the origin of g-'- is connected to 



'Based on Lorentz invariance, relations were also derived by Bukhvostov, Kuraev, and Lipatov in Ref. l9lL 
I93.I93II . It is at present unclear why their relations are different from the relations given here. 
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the non-vahdity of the Lorentz-invariance relations, a measurement of g-'- or checking 
the Lorentz-invariance relations (Eg. 12.741 would be interesting. In the next chapter it 
will be pointed out how g-'- can be accessed. 

2.6 Quark fragmentation functions into spin-| hadrons 

The introduction of the fragmentation functions proceeds analogously to the introduc- 
tion of the distribution functions. A set of light-like vectors is introduced such that 
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— n_ ■ n+. 


n_ ~ n+. 


Ph 


Ml 


4' ^ e^-^^n+pn-,,, 


/; 


^f^-n>'y_-n\n^_. 


A^ = g^^Ay, for any A. (2.75) 




the basis in which 


parton fragmentation functions are defined 



Also here we define for every vector A the Lorentz-invariants A* = A- n^:. The spin of 
the observed hadron is decomposed as 

PJ. Ml, 2 2 

Sh^ShL—n--ShL7r^n++ShT with 5 ,,^ -H S,,^. = 1 . (2.76) 

Similarly as for distribution functions, Mi,/Pj^ does not have to be small. However, 
calculations of cross sections will be constructed in such a way that Pj^ » M/,. 

In general one encounters integrated and unintegrated fragmentation functions. The 
integrated fragmentation functions will be discussed in the next subsection. For the 
unintegrated functions we will see that the two different link structures (as for the dis- 
tribution functions) produce problems with universaUty 13511 . In the second subsection, 
these functions will be parametrized. 



2.6.1 Integrated fragmentation functions 

The integrated correlator which appears in cross sections is expressed as^ 



A,/.-.,P.,S„,4i;/^/^.'.".'-„<.,„5„:.* 



(0)|Q)e 



X (fi|£"^' ""(0+, 77+)^,(/,)|P/„ 5ft; Px>out,c|,-=o, (2.77) 

17=0 



*In order to write the gauge link in one of the matrix elements the light-cone gauge was chosen in an 
intermediate step. 
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where a factor 1/3 is introduced in the definition to average over the initial quark's 
color, and where the color indices are contracted. Its analytical structure satisfies the 
following constraints 

A'(z-\ Ph,Sh) = / A(z-' , Ph,Sh) / (hermiticity), (2.78) 

A(z-\ Ph,Sh) = / A(z-' , n, -Sh) / (parity). (2.79) 

For fragmentation functions there is no constraint from time-reversal. In the case of 
distribution functions we have that T\P)in - \P)o\ii - \P)ia, but for fragmentation func- 
tions there are several particles in the out-state which might interact with each other. 
Therefore, such an identity does not hold and integrated T-odd fragmentation functions 
appear (in contrast to the integrated distribution functions). This mechanism to gener- 
ate T-odd functions was introduced by Collins in Ref. 12311 and was shown to exist in 
model calculations by Bacchetta, Kundu, Metz, and Mulders in Ref. 11971 19811 . 
We continue by defining T-odd and T-even for fragmentation and write 

Mz-\Ph,Sh) = AlT-™l(z-',P,„5,) + A^''-°^''\z-\Ph,Sh), (2.80) 

where AtT-^^veni^^-i ^ P/„ 5/,) and AP^'°'^^\z-\ Ph,Sh) obey 

A'^''"^^^'{z-\Ph,Sh) = {-){iysC) A'^''"^^\z-\Ph,Sh) {iy.C). (2.82) 

This gives the following parametrization of fragmentation into a spin-i hadron 



A(z \Ph,Sh) = -\Di i- - S],L Gi 1/-J5 + H, ^/iT-'Z-rs) (twist 2, T-even) 



(twists, T-even) 



zP^\ 2 ) 

Ml, I no- [*^->'/+]\ 

+ —^\Dt ej jpShT^ - ShL El lys + i H — I (twistS, T-odd) 

-I- higher twist. (2.83) 

parametrization of the integrated fragmentation correlator 



where all functions depend on z and a renormalization scale. 

The functions Di, Gi, and Hi have similar interpretations as the analogous distri- 
bution functions. The function Di describes for instance how an unpolarized quark 
(being either red, green, or blue) decays into a hadron plus jet (over final state col- 
ors is summed). T-odd effects appear at subleading twist for fragmentation into spin- j 
hadrons, in contrast to fragmentation into spin-1 hadrons where T-odd effects already 
appear at leading twist (see for example Bacchetta, Mulders 19911 '). 
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2.6.2 Transverse momentum dependent fragmentation functions 

Gauge invariant correlators and T-odd behavior 

Similarly as for the unintegrated distribution correlators one encounters fragmentation 
correlators with two different links, defined through 



'^ ,1^1, i lu^ hi 



hYj\z-\kT,Ph,Sh) 



" 3 4^ J {InflEv, j (271) 



f 

X out<PA, 5,; Px\^j{Q)I?'- ""(0, ±cx,+)i:±-". ""(Or, cx,r)|Q>c 
X {Q.\J:^^^' ""(oor, 77r)i:°^' ""(+«,+, 77+)(A,(77)m, 5,,; Pz)out.c| ,r=o ■ (2-84) 

When studying the analytical structure of this correlator one finds 

AW"'"(z-',^r,^A,5,,) = / At^l(z"',fcr,n,5,,)/ (hermiticity), (2.85) 

A^^\z-\kT,Ph,SH) = / A^^\z-\-kT,Ph, -Sh) r° (parity). (2.86) 

As for the integrated correlators, the time-reversal operation does not lead to additional 
constraints. This means that the functions appearing in the two different fragmentation 
correlators, A'^^', cannot be related. That holds for the T-odd functions as well as for 
the T-even functions. Since the functions become universal after an integration over kj, 
there could be a universality relation for the unintegrated functions but at present such 
a QCD-relation is unknown. This forms a problem for the by Collins |23] suggested 
method of accessing transversity via the Collins function Hf (Collins effect). 

The problem with universality comes from the interplay of two effects, the final- 
state interactions (in and out-states) and the gauge link. If one of the two mechanisms 
would be suppressed the situation is simplified. For instance, if final-state interactions 
in the out-states are suppressed, then T-odd functions will enter with a sign depending 
on the gauge link or process similar to the situation for the distribution functions. On 
the other hand, if gauge links do not influence the expectation value of matrix elements, 
then fragmentation functions will be the same in all processes (no gauge link means no 
process-dependence fr om th at source). The latter scenario has been observed by Metz 
in a model calculation llOOll and has subsequently been advocated by Metz and Collins 
in Ref. IIOIH . In the next chapter the discussion on universality will be continued. 

The parametrizations 

Transverse moments of fragmentation functions are often encountered. For a function 

Di{z,z^kj), it is defined as 

Df\z) = z^ fdhr ^ Di(z,zhl). (2.87) 

J 2Ml 
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The correlator is decomposed as (see also Ref. I3a.l79l[l()2lll03ll ') 
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(twist 3, polarized) 
(twist 3, polarized) 

(twist 2) 

(twist 3, unpolarized) 

(twist 3, polarized) 

(twist 3, polarized) 

(twist3, polarized) 

(2.88) 



parametrization of the fragmentation quark-quark correlator 



where all functions depend on z and z^k^. Also here a renormalization scale is involved. 
In order to address the universality issue, it is important to measure T-even and T- 
odd fragmentation functions in different processes. For fragmentation functions a set 
of Lorentz invariance relations has also been put forward. For the validity of these 
relations, the same issues as discussed in the previous section play a role. 
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2.7 Summary and conclusions 

We introduced several Cartesian bases and expressed head-on cross sections in terms of 
Lorentz-invariants. It was shown how frame-independent observables can be defined, 
simplifying comparisons between theoretical predictions, experimental observations, 
and fixed frame definitions which are already in use in the literature. 

Two approaches were discussed to access the parton distribution functions. The first 
approach is the operator product expansion. Although this method has a firm theoreti- 
cal basis its applicability turns out to be limited. The other approach, the diagrammatic 
expansion, is a field theoretical extension of the parton model. Although slightly less 
rigorous than the operator product expansion, it can be applied to most hadronic scat- 
tering processes and agrees with the operator product expansion when applicable. The 
diagrammatic approach will be used in the rest of the thesis. 

It was indicated that an infinite set of diagrams, appearing in the diagrammatic 
expansion, can be rewritten such that transverse momentum dependent correlators in- 
cluding a gauge link appear. These correlators underlie the definition of parton distri- 
butions, and provide important information on the partonic structure of hadrons. The 
definition of these correlators contains a bilocal operator and a gauge link which ensures 
invariance under local color gauge transformations. Since the path of the gauge link in- 
troduces a directional dependence, new parton distribution and fragmentation functions 
were discovered at twist three (for related work, see Goeke, Metz, Schlegel 17911 '). A 
measurement of these new functions would contribute to the understanding of the the- 
oretical description. 

Using time-reversal, the functions can be divided in two classes, called T-even and 
T-odd. It was shown that the first transverse moment of T-odd distribution functions 
corresponds to the gluonic pole matrix element, which Qiu and Sterman suggested 
to explain T-odd effects I84l 18511 . As discovered by Collins |25], it was also found 
that T-odd distribution functions appear with opposite signs in the parametrization of 
correlators which have a gauge link via plus or minus infinity. In chapter0]we will see 
that other link structures in correlators can appear, giving more complex factors. 

For fragmentation functions a universality problem was encountered due to the two 
possible mechanisms to produce T-odd effects, the final-state interactions and the pres- 
ence of a gauge link. This means that the value of unintegrated fragmentation functions 
and their transverse moments can be different for different experiments, forming a po- 
tential problem for extracting transversity via the Collins effect. It is therefore important 
to compare fragmentation functions which are measured in different processes (for ex- 
ample their z-dependences). If for instance one of the two effects is suppressed, a simple 
sign relation (plus or minus) should appear between functions which have gauge links 
via plus or minus infinity. The discussion on the universaUty of fragmentation functions 
will be continued in the next chapter 
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2. A Outline of proof of Eq. |238 



For the interested reader a derivation of Ea. l2T8l will be pres ented here. The derivation 
is based on section 4.2 and 7.2 of Peskin and Schroeder 113 ill . 

The invariant amplitude can be obtained via the S-matrix 

5=1 + iT, (2.89) 

free<^X, Ph, l'\iT\P, Ofree = (2;r)V {I + P - I' - P, - P^) iM, (2.90) 

and the LSZ reduction formula 

free(^X. Ph7^S\P, /)free 

iZg iZp l/jQ IjLfj iZx 
•^ 

P - ml + ie P^- M^ + ie l'^ - m] + ie P\ - M\ + ie P\-M\ + ie 
~ fd'^xi d^xz d^i dV2 d'y, e'(''>'.+n>.+Pxy3-P.v,-/.2) 

where the Z,'s represent the field-strength renormalizations, |0) represents the physical 
vacuum, T is the time-ordering operator, and ~ means that the two expressions agree 
in the vicinity of the poles. These poles are generated by the unbounded interval of the 
integrals and correspond to the asymptotic incoming and outgoing states. The task is 
now to work out the right-hand-side and to identify these poles. 

The time-ordered product can be calculated in the interaction picture. In the inter- 
action picture the Heisenberg fields are decomposed in creation and annihilation opera- 
tors at some point in time, fo, which also defines the quantization plane. Evolving these 
operators in time by the Hamiltonian defines the interaction picture fields. Since the 
creation and annihilation operators can be interpreted as asymptotic states it is a priory 
not clear whether this expansion is valid for quarks and gluons in QCD. Since our final 
result in this section does not depend on the explicit expansion of these creation and 
annihilation operators, we will skip here this technical point. Assuming the vacuum 
structure of QCD to be simple and indicating the fields in the interaction picture with a 
superscript /, we obtain 



Eq. (IIU) ~ lim r d^xi d^xz d^yi d^z d^3 e'(''"> +^*.v2+p«'3-p., -to) 



<oi^[ ""1i!''"' ] K(y2Wx(y^)u(T, -t) 



,l,'p{xOU(P,S) 
IM 



10) 



x[<o|c/(r,-r)|0)r', (2.92) 
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where 



<D(x) = f/(x",fo)1>(x)f/^(x",fo) 



t^^o 



Uitb,Q = Texp 



n, 
-i fdt'H,{t') 



(2.93) 
(2.94) 



and where Hj is the interacting part of the Hamiltonian in the interaction picture, and |0) 
is the bare vacuum. The interaction Hamihonian contains the electromagnetic and QCD 
interactions. The electromagnetic interactions can be treated in perturbation theory 
while the QCD interactions in general cannot. 

To produce nonforward matrix elements there should be at least two electromag- 
netic interactions, once between the incoming and outgoing electron and once between 
the incoming and outgoing hadrons, giving 



tr,,{Px,PhJ'\iT\P,l)tre, 

iZg iZp 



iZg 



iZh 



iZy 



P - ml + ie P2-M^ + ie l'^ - m^ + ie P^ - Ml + ie P\-Ml + ie 



fdVdSdVdVdVae' 



x(0|r 



»(/',^')iA'Cvi) 



7-^C»(l-,£) J 

'^iiy2Wxiy^)U{T,zi)A'p{z,)efjzi)U{zuZi)Al{z2) 



X erj(z2)U(z2, -T) 



2M 






|0)[<0|f/(r,-r)|0>r 



(2.95) 



where 7^, ^(z) - Qejq il/^iq{z)y''^[i iz). Using Wick's theorem the lowest nontiivial 
order in e can be worked out, giving explicitly 

Eg. (11951 

= rd4:cid4;c2dSidVd4y3e'('''''+^''>'2+^«''-^-^>-'-'2)(-02 lim fd^zi d^zj 



x{0\T 



u(l'.A')^'jyO 



Im^ 



efjzi) 



¥Mim,X) 



Inig 



|0)<0|rA^(zi)A^(z2)|0) 



X {0\T<^',{y2Wx(y^)U{T,Zi)U{z,,Z2)eJl,{z2)Uiz2, -T) 

xmu{T,-Tmr'+o{e') 



i!/'jxi)U(P,s) 



2U 



10) 



f- 



4 ^i{Pkyi+Pxyi-Pxi) 



5 -, (-ie)u(l\ ^')ypu(l. A) fd\i dV d' 

+ le l'^ - m^ + ie J 

-i) lim fd'z2e-"i'^-^{0\T(i>i(y2W^(yi)U(T,z,)U(zuZ2)eJj(z: 



x( 
xUiz2,-T) 



4''p(Xi)U(P,S) 
2M 



\0)mU(T,-Tm]-'+O(e% 



(2.96) 
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In order to obtain the singularities in the hadron energies it is assumed that the incoming 
hadron only starts interacting after some point in time, f,-, and that the outgoing hadrons 
are well-separated wave-packets after some point in time called tf. These assumptions, 
known as adiabatically switching on and off the interactions, restricts Zj to be f,- < 
Zj < tf- By doing so, the poles of the incoming nucleon and outgoing hadrons can be 
identified after which the limits f, — > -oo and f / — > co can be taken. One obtains 



Eq. (|221 = 



P- -m\ + ie l'^ -m\ + ie 



(-ie)u(l' , A')ypu(l, A)- 



X rd'*xid%2d%3e'(^"«+^«'^-'^-">(-0 lim lim f 

J t,^-"^ 7-^(oo-,e) J 

Ij — >oo ^ 

X mT<^i(y2W^(y3)U(T, f/)[f/(f/, tdeJ-[(z2)]u(ti, -T) 

xmu{T,-Tmr'+o{e') 

i i iZp iZv 



d Z2 e 



iA^(xi )t/(P,S) 

2M 



10) 



iZi, 



P- - m\ + ie l'^ ■ 



■ ml + ie P^-M^ + ie p^ - M^ + ie PI - MI + ie 



x(-ie)u(l',A')ypU(l,A)(-i) lim 



fd^zie 



-iqzi 



—I 



x,,,{Pi„Px\TU(to,tf) 



U(tf,ti)eJ^(z2) 



U(ti,to)\P,S},,,,+0{e^). (2.97) 



Note that all the necessary poles including the field-strength renormalizations at order 
e^ (for nonforward matrix elements Z^ = 1 at this order) have been produced. After 
shifting the current, J'^izi), yielding the delta-function which expresses momentum 
conservation, we can read off the invariant amplitude M from Eg. 12.901 giving 



iM - (-ie)u(l', A')ypu(l, A) 



X Um -jout{Ph,Px\U(to,tf) 



f/(f/,f,)(-/e)7;(0) 



= {-ie)u{r,A')ypu{l,A) — ,,i{P,„Px\{-ie)nO)\P,S)n.,c +0(e') 
Note that the current, JP, is here in the Heisenberg picture. 



U{ti,to)\PS},n,c+0(e') 

(2.98) 



2.B The diagrammatic expansion 



For the interested reader a short outline of the derivation of two diagrams in the dia- 
grammatic approach will be given here. It will be pointed out how the results can be 
generalized. 
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Following the same procedure as applied for semi-inclusive DIS in the previous 
appendix, one can derive that the invariant amplitude for general processes reads 

(free states P|iT|free states Q)= lim out<f |t^(fo, tA U(tf, f,) t/(f,-, fo)l2)mc, (2.99) 

tj^f — CO -^ ^ ' 

where the f/'s are defined in Eg. 12. 941 When applying the diagrammatic expansion one 
expands the bracketed term. 

At tree-level the result is expressed in terms of fields inside matrix elements. These 
fields can always be transported to the origin, yielding momentum conservation for 
the incoming and outgoing particles. The lowest order result is the tree-level squared 
Feynman diagram with spinors or polarization vectors replaced by Heisenberg fields in 
matrix elements, multiplied by a delta-function expressing momentum conservation. 

For example, consider hadron-hadron scattering with momenta P\ and Pi produc- 
ing two hadrons with momenta K\ and K^ approximately back-to-back in the azimuthal 
plane. It is assumed that the outgoing hadrons are well separated from the incoming 
hadrons. If this is not the case then fracture functions are needed to describe the com- 
bined process of distribution and fragmentation. Although the incoming and outgoing 
partons can be quarks or gluons, we will restrict ourselves to (anti)quarks as external 
partons, the gluons can be included later on. In that case the incoming quarks or anti- 
quarks interact via a gluon giving quarks and antiquarks in the final state. Expanding 
the bracketed term of Eq. l2.99l to lowest order this contribution is 

= (/g)2rd4xd43;out<^i,^2,^|r^(x¥(x)«A(x)^(3'¥(y)(A(y)l^i,P2>i„,c+<3(5'). (2.100) 

Overall momentum conservation can be simply obtained by making an overall shift of 
X and subsequently redefining y. One finds 

{KuK2,x\iT\PuP2) = \aS ou,{KuK2,x\Tip(Qmo)Hmiymy)Hy)\PuP2)in,c 

X (2n)''6\P,+P2-Ki-K2-Y^Px:)+0{g^). (2.101) 

Since no gluons appear in the initial or final state, those fields need to be contracted 
with each other yielding the gluon propagator (Wick's theorem). Each of the quark- 
fields can be connected to one of the incoming or outgoing jets. All possibilities should 
be included. Picking one of the contributions as displayed in Fig. I2.12h and working 
out the contractions gives (k\ = Px,+Ki, ki = Px^+Ki, pi = Pj-P^.) 



Eq. (EHD ={2n)^S'(Pi+P2-Ki-K2-J]Px.)- 



' (P2 - fe)^ + ie 
X outiKuPxMmn), (ithy") out<^xil<A(0)|Pi>i„,c 

X out{K2,Pxm0m)c {ith7M) out<i^X2|tA(0)|P2)m,c. (2. 102) 
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P2 -k 




P2 - k2\ O 




Figure 2.12: Two contributions which appear in the diagrammatic approach. 



If more interactions are included these steps can be performed in the same manner. 
Each interaction in U consists of fields which are fully integrated over their independent 
coordinate variables. As an example we will consider one higher order contribution 
which will yield Fig. 12. 12b . The interaction can be expressed as 



{KuK2,X\iT\Pi,P2} 

= (ig)' fd\ d^ out{KuK2,X\T;f,(zmzmz)mm)4'm(y)^(y)i/^(y)\PuP2}u 



{2nf6\P,+P2-K,-K2- Y^ Px,) + . . . , (2.103) 



where the dots denote contributions from other interactions. We need to include all 
possible contractions but we will consider now one of the possibilities. Contracting 
ifi(z) with (/'(O) and A(0) with A{y) and considering the remaining fields to be connected 
to the jets in a certain way, we find 



One contribution of Eq. ( I2.103> 
= (2n)'6\P,+P2-K,-K2- yPx^)8' [d'zd'y ^ 



^' d^ ,,,, ,,--! 



X out<^l , Px, \^{z)\0)ciyata 



iJk + m) 
I2 - m^ + ie 



irtboul{PxMmPl)in 



X out<^2, PxMy)K(z)\^')ci7,tbout{PxM0)\P2}in,c (2.104) 

We would like the fields which already appeared at "tree-level" to be at the origin. 
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Shifting those fields and redefining the other integration variables, one finds 






(l2+kiy-m^ + le 

(2.105) 

The result is the original expression plus an additional interaction which contains in- 
tegrals over z and l2. We are actually deriving the Feyman rules. If we would have 
considered a virtual correction there would have been an additional contraction and one 
would remain with one integral over l2. 

In general one should have the following situation: the tree-level result consists of 
momentum conservation and matrix elements containing fields at the origin; if more 
external partons^ are present then the corresponding fields are integrated over their co- 
ordinates and momenta, and their interactions are described by the Feynman rules. Vir- 
tual corrections are similar except that the integral over the coordinates can be carried 
out. 

The cross section is obtained by removing the delta-function related to momentum 
conservation, taking the square, and multiplying with the removed delta-function (see 
Eq. 12.151 12.891 I2.90t . Except for the phase-space and flux factors this is the cross 
section. In general the total momentum conservation can be rewritten by giving each 
matrix element a number (/) and naming the sum of the outgoing momenta (also Pxi 
here) Ki^^i and incoming momentum /",;„. The total momentum conservation can then 
be written as 



(27r)454 



;- 1 ;_ 1 ^ .*J 



V'=i j=l 



X {2n) 6 (external parton momenta), (2. 106) 



where f, are the momenta of the external partons and where the last delta-function ex- 
presses the momentum conservation of the external partons. The first n delta-functions 
are used to shift the fields of the matrix elements by introducing another integral over 
the coordinates 

6\t„ ± P„;„ + K„,J = I iL%,'('»-^,..n-^,.o.)f„ . (2. 107) 



J (2nr 



After an integration over the unobserved momenta one obtains the correlators as intro- 
duced in sectionlZ4l 



'External partons are incoming or outgoing partons connected to correlators or decaying as a jet. 



3 

Electromagnetic 

scattering processes at 

leading order in as 




The diagrammatic approach, introduced in the previous chapter, will be employed to 
describe at leading order in as cross sections in which the hard scale is set by an elec- 
tromagnetic interaction. In the discussion on the meaning of ag it will be argued that 
even at leading order there is an infinite amount of interactions that should be consid- 
ered. Including next-to-leading order corrections in inverse powers of the hard scale, 
these leading-order-ffs -interactions will be evaluated in order to obtain the cross sec- 
tions for semi-inclusive DIS, electron-positron annihilation, and Drell-Yan in terms of 
distribution and fragmentation correlators. These correlators contain gauge links of 
which their paths are different for different processes. Using the parametrizations of 
these correlators some explicit asymmetries for semi-inclusive DIS will be given. 

In the second part of this chapter, the universality of fragmentation functions will be 
discussed and the link structures appearing in deeply-virtual Compton scattering will 
be studied. The calculation of the latter illustrates the wide-ranging applicability of the 
diagrammatic approach. A summary is provided at the end of this chapter. 



50 



3 Electromagnetic scattering processes at leading order in a^ 



Pht4 V;. 



Phf4 V/« 



Pj Vh 







(d) (e) 

Phti ^'^Ph PhU' V/^ 





(b) (c) 

P,J/ Vh Ph// Vh 

q & ^ q q / \ q 




Figure 3.1: Various interactions between the two jets and the elementary scattering part 
in semi-inclusive DIS. 



3.1 Leading order in as 

In the previous chapter the diagrammatic expansion was developed to handle scattering 
processes. In order to apply this expansion, the various jets need to be well separated. 
This can be further translated into the demand to have large momentum transfers, im- 
plying that the inner products of momenta connected to different jets are large. Besides 
an expansion in a^, this presence of large energy-scales allows for a possible expan- 
sion of physical observables in the inverse powers of that energy-scale (M/Q). Both 
expansions will be made throughout the chapter. 

Applying the diagrammatic expansion, all leading as diagrams which can con- 
tribute to the processes need to be included. For semi-inclusive DIS, the lowest order 
(in g) diagram is given in Fig. l3.lh . Interactions which can be absorbed, like the one 
in Fig. 13. lb . are already present in the correlator definition and should not be included 
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separately. Diagrams which have several jets are a, suppressed compared to diagrams 
which contain less jets. The restriction to leading order in a, leads to only two outgo- 
ing jets for semi-inclusive DIS, electron-positron annihilation, and Drell-Yan. Further, 
the interactions between the jets as illustrated in Fig. 13.1b are considered to be weak 
and will be neglected. This assumption is directly related to the question on whether 
factorization holds for the correlators. This will be discussed in the next chapter. 

Another possible contribution is that a gluon-line will connect a jet with the elemen- 
tary scattering diagram. Since there can be interferences, such as depicted in Fig. 13. IB . 
such an interaction is proportional to the coupling constant, g, and is not Os suppressed. 
We shall see in the next section that one can resum all collinear gluon-lines connecting 
the correlator and the scattering diagram to all orders in g. The other possibility of 
additional collinear quarks, such as in Fig. 13. lb . is not taken into account here. Ellis, 
Furmanski, and Petronzio |53, 54], found that such contributions are suppressed with 
M^/Q^ and are therefore beyond the aimed accuracy of this chapter. 

Besides additional parton-lines from the correlator one might wonder whether there 
are other corrections in the scattering diagram itself which need to be included at lead- 
ing order a.,. Such corrections consist of loops and gluon radiation and can partly be 
absorbed in the correlator definitions. For instance, if the loop diagram in Fig. 13. IF 
contains a far off shell gluon, then it cannot be absorbed. However, if the gluon is ap- 
proximately on shell and collinear to one of the jets, as illustrated in Fig. 13. lb . then it 
is absorbed in the correlator definition up to some scale. For gluon radiation a similar 
situation applies. If the gluon in Fig. l3.lh has a small transverse momentum with re- 
spect to the parent hadron (P), then it is already present in the correlator definition. The 
gluon in Fig. l3.1L cannot be absorbed and represents an as correction. 

Absorbing interactions in the correlator introduces a particular scale. There were 
a few attempts to calculate this scale dependence (evolution equations), but including 
effects from intrinsic transverse momentum this has not been achieved so far For a 
discussion on the scale dependence and various other difficultie s the reader is referred 
to Ellis, Furmanski, Petronzio f54T, Collins, Soper, Sterman I104i llOST . Hagiwara. 
Hikasa, Kai 1 106], Ahmed, Gehrmann 1 107], Boer ] 108], and Henneman I66il. 

The conclusion of this section is that, restricted to leading order in as, one must 
resum over all possible gluons which connect the correlators and the scattering diagram 
to obtain the cross section (generalizations of Fig. 13. li d). Such resummations will be 
performed for several processes in this chapter 
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In the first two subsections the expressions for the hadronic tensor including M/Q cor- 
rections will be derived. In the third subsection the results for the hadronic tensor will 
be used to express some asymmetries in terms of distribution and fragmentation func- 
tions. We will begin here by introducing the calculation of the hadronic tensor. 
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Applying the general rules of the diagrammatic expansion, the hadronic tensor is 
found to be a sum over various diagrams. Considering only leading order in as, all 
those diagrams (to all orders in g) will be considered that are not already absorbed in 
the correlators. In order to structure the calculation we shall treat here the lowest and 
one higher order contribution in g. Those contributions in combination with a number 
of other diagrams will be studied in detail in the next subsections. The final result, 
Eg. 13.461 will consist of correlators containing gauge links for which parametrizations 
have been given in the previous chapter (see Eg. 12. 57112. 65112.831 I2.88> . 

At lowest order in g (see also Fig. 13. Ih ). the contribution to the hadronic tensor 
is a convolution between the scattering diagram, the fragmentation correlator, and the 
distribution correlator In this section arguments connected to the parent hadron will 
be often omitted for notational convenience, so <!>(/?, f, 5) — > <l>(p). The result for the 
hadronic tensor reads 

IMW" = r dV d'^k 6\p + q-k) T^'^ [^{p)^ A{k)y] + 0{g), (3.1) 



where O and A were defined in the previous chapter (Eq. 12.371 I2.40t . Note that the 
Heisenberg fields in the correlator contain interactions to all orders in g, see for example 
the interaction in Fig. 13. lb . 

In order to evaluate the result, a Sudakov-decomposition is made. A set of light-like 
vectors ({n_, n+] with «_■«+ = ! and «_ ~ «+, where the bar denotes reversal of spatial 
components) is introduced such that 



„ xbM^ ^ Q 
eV2 XB^/2 

p ZhQ ^ Ml 
Ph = —p n- + — :r-p n+, 
V2 ZhQ^/2 

Q + 0(M^/Q) Q+0(M'-/Q) 

q^ F n- F n^+qj, (3.2) 

V2 V2 

Sudakov-decomposition for semi-inclusive DIS 



where Q^ = Q^-^q\ - Q^-\-0{M^). Using the parton model assumptions the correlators 
in Eg. 13.1 l ean be integrated over the small momentum components 



IMW" = I d^pT d^kr 6^(pt + qr - kj) 
1 



Tfifdp-^ip)] yji^i fdk^ A{k)\ : 

\J 'p+=-q* \J lk-=q- 

{\+0{M^IQ-))+0{g)+0{as), (3.3) 



x-Tr° 

3 
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where the relation ^at - Sat Tr'' (5/3 was applied. 

At higher orders in g one needs to include gluonic diagrams in which gluon-lines 
are connected to the correlators. If there is more than one soft correlator present, as is 
studied here, the inclusion of such interactions becomes more complex. A simplifica- 
tion can be obtained by assuming the expressions for the cross section to be color gauge 
invariant, allowing for a suitable gauge choice. A convenient gauge turns out to be the 
light-cone gauge with retarded boundary conditions. In such gauges, see for example 
Ref. 1 26 1, one of the light-cone components is set to zero together with the transverse 
polarizations at light-cone minus infinity, ox A^iif) - Ajij]^, -oo, tjt) - 0. 

In the calculation of the hadronic tensor, we shall employ the equations of motion. 
It has been put forward by Politzer 1 5 1 ] and subsequently by Boer 1 64] that the classical 
equations of motion, (i^-m)t//{x) - 0, hold within physical matrix elements. They lead 
to the following identity for the correlators A and A^ (see also Eg. 12.401 and Eg. 12.411 
and where A^ = A^ f/) 

A(k)()t -m)^-g fd'ki Al(k, k,)y,. (3.4) 



We will briefly come back to the validity of these equations in the next chapter (subsec- 
tion l432l Drell-Yan). 

When including gluons from the fragmentation and distribution correlator, the stan- 
dard treatment shows that gluons which are backwardly polarized (S^ ~ p^) lead to 
0(m^/Q^) suppressed matrix elements and can therefore be safely neglected. This 
means that the A ^-gluons from the fragmentation correlator and the A -gluons from 
the distribution correlator can be discarded. Together with the chosen light-cone gauge 
this means that only the transversely polarized gluons of the fragmentation correlator 
and the longitudinally and transversely polarized gluons from the distribution correlator 
need to be considered. 

We shall continue the calculation by considering a simple higher order interaction. 
Taking a gluon from the distribution correlator as displayed in Fig. 13. IB gives the fol- 
lowing contribution to the hadronic tensor 

IMW" = 2MW^;g i^uj^ + other diagrams, (3.5) 



xTr' 



D,C 



' (k - piY - m^ + le 



(3.6) 



where O^ is now the quark-gluon-quark correlator as defined in Eq. 12.381 In the next 
equation the above expression is rewritten to indicate how the various parts of the quark- 
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propagator in combination with the gluon polarization contribute to the cross section 

2^^Rg mirf = / dV d'k 6\p + q-k) 

xilMv/", +2Mm/"' . +2M\\^\ , +2Mv<'' , V (3.7) 

2A^<'/ = f— ^-^^^^^Tr«'^[<l.X(p,pi)/A(^)r-f,0^ + m)y''], (3.8) 
L,)(+m J (k-pxY-m-^+ie L / J 

2^^"/ = r r. "^'^^2 ^"?^- Ti-^'^ [<^Up,Pi)rmy-t, (-//ir)y1 , (3.9) 
^•"''iT' J (K-z^O'^-m^+ie L ' J 

2M^r = [-IZliz^Tr^.cb^ -( )yA(^)y„f, (^-/i^+m) y"! , (3.10) 

2M< _= f-iZlll^ Tr^.c [a)^,-(p, pO/A(fe)r„f/ (-/^Ir") /] . (3. ll) 

As we will see later in this section, the various terms contribute to the following or- 
ders in MjQ: IMv/^ , contributes at leading order to the longitudinal gauge link, 



IMv/'' , contributes at subleading order, IMv/'' , , contributes at leading order 
to the transverse gauge Unk, and 2My/^ _ + _ contributes at subleading order. 

In the following subsections we will evaluate these terms and include higher order 
gluon insertions from the distribution correlator We will begin by analyzing the leading 
order in detail which leads to the gauge link; the next-to-leading order in MjQ will be 
discussed in the second subsection in which transversely polarized gluons from the 
fragmentation correlator contribute as well. 

3.2.1 Leading order in MIQ 

Longitudinal gauge link 

We will study here in detail the contribution of a longitudinally polarized gluon inserted 
in the diagram. This follows the leading order in MjQ calculations of Bjorken, Kogut, 
Soper 1 109]; Efremov Radyushkin 1 1 10]; and ColHns, Soper, Sterman ] 1 IL .llZ.lll . 
In Boer, Mulders 111411 those calculations were extended by including MjQ corrections 
in the diagrammatic approach. In that paper explicit calculations were given to order ^ 
which were generalized by using arguments based on Ward identities. 

When studying this problem Ward identities should be handled with care ]39]. The 
considered gluon (with momentum /?i) is approximately longitudinally polarized and is 
inserted at various places in an amplitude, seducing one to use p'lAipipi) = 0. However, 
in the limit of pi — > the Ward identity does not contain any information (0 = 0). 
Since this issue especially becomes relevant in more complicated diagrams (chapter|3 
explicit calculations will be performed in this chapter. 
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Figure 3.2: Various contributions to the hadronic tensor The interactions in the figures 
cannot be absorbed in a single jet and therefore need to be considered. Figure (c) is not 
of the type of figure ITTI d because the interaction is connected with the hard part. In 
figure (a) only the indicated part of the numerator of the quark propagator is considered. 
The part which is not indicated will be treated in the next subsection. 



Inserting the gluon on the left-hand-side of the cut, see also Fig. l3.2h . we need to 
evaluate Eg. 13.81 and Eg. 13.91 The contribution from Eg. 13.91 will be calculated in the 
next subsection and turns out to appear at subleading order To perform the integral 
over p\, the denominator {{k-pi)^-trp-+iey^ in Eg. l3.8l is simplified into (-Ip^k'^ + 
/e) '. This simplification, called the eikonal approximation, is justified when making 
the parton model assumptions which were discussed in section 12.41 The integral over 
pt in IMv/^j (Eg. l3.8> takes now the following form 



L/ 



/ 



Ap\ 



Jp^i'i-^) 



-2p'^k + ie 



AliT]). 



(3.12) 



Assuming k to be positive, the integral can be performed by calculating the residue. 
This leads to a Heaviside function {6) 



-0{vi--f)A^{rf). 



(3.13) 



The Heaviside function expresses that there will be only contributions to the hadronic 
tensor when ri^>^^, meaning that A^(?7"— >oo) contributes, but 4^(7/"— >-oo) does not. 
Using the equations of motion one finds for 2Mvv'"' , (y^(/t+m)-2k^-(^-m)y^) 

L,jt+m 



2M^/; = rAl^e'>^<P,5|^(0)/Tr' 
x[\+0{m^IQ')) 



A(k) 



rX-ig) 



r 

)jdj]-A^{r]-,e. 



JdViTr^''' ^X{p,Pi)rJd''k, A^^{k,ky)ypt, 



—2k p| + ie ' 



^T)iif(mp,s)c 



(3.14) 
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where the relation A„b - 6„b Tr'' A(k)/3 was used. The first term above contributes 
at leading twist and is exactly the first order gauge link expansion running via infinity. 
The second term in Eg. 13. 141 coming from applying the equations of motion, is canceled 
when the other diagrams of Fig. l3T2l are included. Although this cancellation' occurs 
in the chosen light-cone gauge, where the gluon propagator has the numerator (see for 
example Ref. L26.1 ') 

d^'iD^g^^--^-^, (3.15) 

one does obtain another term proportional to A^^i. We see that by considering the other 
diagrams the A^ term got replaced by a A^^i term. This cancellation was observed in 
Boer, Mulders 1 114]. Since this issue only plays a role when considering fragmentation 
correlators (a free outgoing quark does not produce Aa terms), this kind of cancellation 
is expected to hold to all orders. One should find by including higher order diagrams 
that the Aaa term gets replaced by a Aaaa term and so on. This cancellation, which 
is expected to hold at each order in g, has not been proven to all orders and deserves 
further investigation. 

Before continuing it should be pointed out that the cancellation also occurs when 
the Feynman gluon propagator is taken in the hard part. This suggests the hard part and 
the correlators to be separately gauge invariant, which is a minimal requirement for a 
factorized description. 

In order to generalize the above result. Eg. 13. 141 we consider « longitudinally polar- 
ized gluon insertions on the left-hand-side of the cut which contribute to the hadronic 
tensor as 

IMW" = r (fp d'^k S\p + q-k) 2Mv/\, + other diagrams, (3.16) 

2A^<I^ = fdVi ■■■dV Tr°'^[(D+;;t^J/7,pi-/72,...,p„-i-p„,/7„)/A(^) 

xi-g) (ti.j )-, — Tj — r—---'^'h7)7, — T2 — rr^y ■ ^^-^^^ 

(k-pn) -m^+i£ (k—piY-m^+ie 

Performing all the ^/-integrals one straightforwardly obtains an ordered product 

2M<I^ = / ^/PHP,Sm)r^r^ ^ fi-igf Jd/^r Jd;72 . . . Jd,?,; 

xAUT]-i,e,^T).-A^(T]-,t,^Tm^)\P,S\(l +0{MIQ)), (3.18) 
which is the n*-order longitudinal gauge link expansion. 



'in order to achieve this cancellation it was assumed that the gluon connecting the fragmentation correlator 
is outgoing and approximately on its mass-shell (k^ > 0). 
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Figure 3.3: Various contributions to the hadronic tensor. In each of the diagrams it is 
indicated which part of the numerator of the quark propagator is considered. 



To summarize, by including longitudinally polarized gluons from a distribution cor- 
relator which interact with an outgoing quark, one obtains the longitudinal gauge link 
£^T' f (oo ,^ ) (for definition see section IZSl . Inserting longitudinally polarized glu- 
ons from the distribution correlator on the right-hand-side of the cut also produces a 
longitudinal gauge link. Its result can be obtained by taking the complex conjugate of 
the result we just obtained (Ea. l3.18i and interchanging /i with v. For those insertions 
one finds the gauge link X"'^' ^ (0", oo"). 



Transverse gauge link 

We continue by studying the leading contributions from transversely polarized gluons. 
For a long time transversely polarized gluons were thought to appear always suppressed 
in cross sections (see for instance Mulders, Tangerman LllSJ '). In 2002 Belitsky, Ji, and 
Yuan showed in Ref. 12611 that this idea is wrong and that transversely polarized gluons 
at leading order give rise to a transverse gauge link. Although their calculation showed 
how the transverse gauge link could be derived, it also contained a few points which 
deserved further clarification. Some of these points were given attention in Ref. 1 35|] 
but other points were unintentionally left open. All issues will be discussed here in 
considerable detail. 



We start with IMv/" , 



as given in Eg . 13 . 1 01 (represented in Fig. l3.3h ') 






]{k- 



dVi i-K) 



piy—rrp-+i£ 



T/'''']^ArT(P,Pi)y'Mk)87ati (Jt-lAT+'n)r] . (3.19) 



Performing the integral over p^ would lead to J drj^ A'^i^]) as was also encountered 
when calculating the longitudinal link. In gauges where the gluon fields vanish at in- 
finity this would form a perfect matrix element, but in those gauges where the fields do 
not vanish the matrix element is divergent and needs subtractions^. 

^A similar problem occurs also with the longitudinally polarized gluons. In order to keep the matrix 
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In order to consider the transverse fields at infinity more precisely one expresses the 
A-fields as 

A-(/7) = A^c^, 77\ 77r) + ^(t,) - A«(oo, jj\ r^j)] . (3.20) 

In this decomposition the ?; -dependence has vanished in the first term. When perform- 
ing the // -integral over the first term of Eg. 13.201 in Eg. 13. 191 (which is proportional 



to exp [ipiirj - ^)]) a delta-function in p^ emerges-'. This delta-function eliminates the 
large momenta in the denominator in Eg. 13. 191 making the contribution leading in M/Q. 
When performing the /;|-integral one finds that the second term in Eg. l3.20l contributes 
at order M^/Q^ and can be neglected. The choice for taking the subtracting point in 
Eq. 13.201 at infinity might seem arbitrary but actually it is not. Choosing the subtrac- 
tion point at 77 = - oo would produce, after an integration over p^, the divergent term 

J^ dj]^ [Aj(ri)-Aj(-oo^ ,ri^ ,riT)]. Hence it is the ie prescription in Eg. l3. 19l (or process) 
which forces us to take the subtraction point at infinity in order to keep the contributions 
from both terms in Eg. l3.20l finite. 

We continue with the first term and find 13511 



^""^^-Ar^^-f^'P^^'P'^ I 



(f^ d/7+ (fi]T 



,ip(+'PlT(ri-i;)T+'P, in-^y 



X ]im{P,S\i//iO)y'A(k)iyai 



Jt-i/ij+m 



pI^o (k - p\r - Iff- + ie 



X 



gdPjd^T ■AT(c^,r]\(T)HaP,S}c{l+0{M^/Q-)), (3.21) 



where C is some constant. To arrive at an equivalent expression the existence of a pure 
gauge at infinity was assumed in Eq. (38) of Belitsky, Ji, Yuan 126,1 . This assumption 
makes the treatment for QCD less general and is avoided here. 

Having the gluon coUinear means that the delta-function in p^ also pushes the trans- 
verse momentum down. In order to proceed we needed to interchange those limits. 
Therefore, when taking the limit /?| — » we keep pir finite. This exchange of limits, 
which is also present in Ref. I2(tI| and seems to be unavoidable, might lead to problems 
at higher orders in as . 



elements containing tlie longitudinal gauge link convergent one should require the longitudinally polarized 
fields to vanish at infinity. This requirement is supported in those gauges where the gluon-fields at infinity 
can be described with a scalar potential (see also Ref. t26ll ') 



f 



</>{ri) = J df ■ MO- 
c 

In order to have the potential finite at ;;" = 00 which seems to be physically reasonable, the A*-fields need to 
vanish at that point. This argument does not hold for transversely polarized gluons. 
^The boundary terms were also studied by Boer, Mulders, and Teryaev in Ref. l87ll . 
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Doing now a partial integration and applying the equations of motion one finds 
by performing first the 77 -integral, the /7|-integral, and then the rest of the integrals 
(A(k)[-iArM-l/iT+m] = [A(k)Ut-i/,T-m]+gJd''ki A-Jk,ki)ja]Ut-i/,T+m]) 

4 ^'' 

= J^e'^f <P, S |«A(0)/A(fc)/(-/g) Jd^r ■ Ar(^, f, (t)^{^)\P, S }, 

c 
x{i+0{m^/Q'-)) 

+ 8 rii^^^-^lime''V>'("-«<P,5|.A(0)/ fd%Al(k,ki)y, 

•It 

^-piT + 'n 



/ 



/(-/g) d(T ■At{^,i^\^t)^(S)\P,S\, (3.22) 



{k — piY — m? + ie 

c 

where the first term is exactly the first order expansion of the transverse gauge Unk. 
The second term was not considered in Ref. 12611 . while in Ref. 13511 it was thought to 
be suppressed in MjQ. In the case of jet-production in DIS the second term does not 
appear (the outgoing quark is assumed to be free), but when considering the fragmenta- 
tion of a quark into a hadron neither observations of Ref. 126113511 are in general valid. 
The second term is as leading as the first term and will be considered here. 

Similar terms were also encountered when deriving the longitudinal gauge Unk. 
Also here one finds that when including the diagrams in Fig. 13.3b and Fig. 13. 3t that 
the second term is exchanged with a term proportional to Aaa- Similarly as for the 
longitudinal gauge link, it will be assumed that this behavior can be generalized to all 
orders. When taking the Feynman gauge for the hard part, one finds that the diagram 
in Fig. l3.3t has an additional contribution at next-to-leading order in MjQ. In the next 
subsection where we will discuss the subleading order in MjQ, it will be pointed out 
where this term gets canceled. 

Generalizing the result above to all orders by inserting n transversely polarized 
gluons, one obtains for the hadronic tensor 

IMWt"^ = r dV d^k 6\p + q-k) 2Mv^^^^ + other diagrams, (3.23) 

2M< = fdVi . . .dVn '^'^ '^''^':, ' 'f"^" e^-P-P'^^e'^P^-P^^"' . , .e"^P'-'-P'->'''-'e"'"'^" 
nAj I (2n')^"'^ ^ 

X {-gTiP, 5 |^(0)/A(fc)r„„f,„ ^,^ ^^'^ , . . . . r., f/, ■ ' ' ^ 



' (k-pn) -in^+ie {k-piy--m?-+ie 

xfAf'^{i1„)...A"'(mm^)\PS\{l +0(M/Q)). (3.24) 
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Following exactly the same steps as done for the single transversely polarized gluon 
one obtains 



2Mv<;^=jdVi...dV 



4 j4 d'^^d'^T]i . . .d'*77„_i ,., 



'ip-PiX„'{p\-P2)n\ J{p„-2~P„-\)ri„-2Jp„-\ri„-i 



iiny 



An 



X {P, S |l/r(0)/'A(^)y„„_, J— ...Ja, Tj — 



{k-pn-\) -mr+ie (k—p\Y-nv-+ie 

In-lT 



x(-ig) J df„ ■ A(oo,^\(„T)(-gr-'A'''-'(Tj„-i) . . .A"'(,^^)i/,(^\P,S)c 

c 
x{l+OiM/Q)). (3-25) 

Using now an identity which holds for non-Abelian fields'* 113 511 



It in It 



ilT fn-l,r 

It It It 



5« Jd^ir ■A((^)jd^2T ■A((2)...JdUi.T -^(4+1), (3.26) 

identity for non-Abelian fields 



one can subsequently perform all the integrals and obtain 

C (i („-i 

xA"'{<x,,e,(iT) ■ ■ ■A''"(c^,e,i„T)4'(^)\P,S)Al +0{MIQ)). (3.27) 

This is the n*-order expansion of the transverse gauge link X°° ' ^ (C, ^t)- 

The gluon insertions in the conjugate part of the diagram can be evaluated as we did 
for the longitudinal gauge link. Taking the complex conjugate and interchanging /i and 
V one obtains i:~ '^'(0,0. 

The complete leading order result 

In the previous paragraphs we resummed gluons from the distribution correlator on one 
side of the cut being either longitudinally or transversely polarized. When taking the 
combination of longitudinally and transversely polarized gluons on one side of the cut, 
one finds in leading order only contributions from those diagrams in which the gluons 
couple in a specific order. This can be seen as follows. 

••in contrast, Eq. (47) of Ref. |26] does not hold in QCD. 
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Figure 3.4: Two contributions in the diagrammatic approach. The parton-lines entering 
from the bottom of the graphs belong to the distribution correlator. 



Consider a diagram containing a single A^-insertion. The sum over all possible 
longitudinally polarized gluon insertions between the inserted vertex of the transversely 
polarized gluon and the fragmentation correlator (see for instance Fig. 13 .41 a) yields a 
longitudinal gauge link between infinity and the AT-field which is at some coordinate 
77. The leading contribution of this diagram comes when this Aj--field is at infinity 
and therefore this longitudinal gauge link vanishes. Only the longitudinally polarized 
gluons, which couple between the vertex of the A^-gluon and the virtual photon (for 
instance Fig. 13. 41 b). contribute at leading order and provide the longitudinal gauge Unk. 
It is exactly this order which also appears in the full gauge link. 

Taking combinations on the left-hand-side and on the right-hand- side of the cut is 
relatively easy. The contributions on each side can be evaluated without using informa- 
tion from the other side of the cut. One simply obtains a product of both insertions, so 

Up to now we inserted gluons from a distribution correlator to an amplitude and 
looked at their final contribution in the cross sections. It turned out that the trans- 
versely polarized gluons at infinity contribute at leading order and yield the trans- 
verse gauge link. The same calculation can be done for transversely polarized glu- 
ons from the fragmentation correlator. The only leading contribution comes when 
those fields are at minus infinity, but those fields do not contribute in the chosen gauge 
{A^iji) - Ajirj^, -00, rjj) - 0). This means that the transversely polarized gluons from 
the fragmentation correlator always appear at subleading order in this gauge. In order 
to find the gauge invariant fragmentation correlator we assume that the expression for 
the cross section is gauge invariant. In that case, there is only one link operator possi- 
ble for the fragmentation correlator which vanishes in this particular gauge. That is a 
gauge link running via minus infinity. This gauge link is also found when one would 
have chosen to work in A^ = Aj-(oo, 0, ^t) - gauge. 

In this process we have found that inserting gluons on an outgoing quark yields a 
gauge link via infinity, while coupling gluons to an incoming quark gives a gauge link 
via minus infinity. In other processes we will encounter the same behavior Coupling to 
outgoing parton-lines leads to gauge links via plus infinity, while coupling to incoming 
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Figure 3.5: Two projected single transverse target-spin asymmetries measured by HER- 
MES Q. The following cuts were applied: Q^ > 1 GeV^ W^ > 10 GeV^ and 
0.1 < y < 0.85. According to Boer, Mulders 117 311 the asymmetry in figure (a) is pro- 
portional to /iii/j while the asymmetry in figure (b) is proportional f^^: ^Di (in the 
scaling limit). 



parton-Unes leads to gauge links via minus infinity. 

The hadronic tensor at leading order in MjQ can now be expressed as (the color 
factor 1/3 has been absorbed in the fragmentation correlator) 



IMW' 



-s 



(T-pT (fkr 6^{pT + qT- kj) Tr 



x{\+0{MIQ)) + 0{as). (3.28) 



This result includes an infinite amount of diagrams at leading order in MIQ and is al- 
most equivalent to the result of Boer, Mulders |78|. The only difference is that the 
correlators above are fully gauge invariant (without assuming Ar-fields at oo" to vanish 
inside matrix elements), forming a natural explanation for transverse momentum de- 
pendent T-odd distribution functions. By integrating the hadronic tensor over qj, these 
T-odd distribution functions disappear, because the resulting /77-integrated distribution 
correlator contains a gauge link of which its path is just a straight-line on the light-cone 
and does not produce T-odd effects. The integrated distribution functions ei{x), frix), 
and h{x), as introduced in Ref. 1781 . are therefore zero within this approach. 

From the hadronic tensor, like the above expression, several asymmetries were cal- 
culated in Ref. |78]. Two of these asymmetries were recently measured by HERMES. 
In the scaling limit (Q^ — » 00) the asymmetries given in Fig. l3.5h and Fig. 13.5b become 



Ml) 



proportional to the Sivers function, f^j ', and the Collins function H 



7-L(l) 
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3.2.2 Next-to-leading order in MIQ 

In the previous subsection we restricted ourselves to leading order in MjQ. Here the 
calculation will be continued to next-to-leading order We will start by analyzing in 
detail the contributions from a single gluon insertion at order in MIQ (Eg. 13.91 and 
Eg. 13. in and include the other contributing diagrams later on. 
Evaluating Eg. l3.9l (see also Fig. l3.6h ) one finds 



A(^) 






2Mm^; = ril^e'>^(P,5|^(0)/TrC 

CO 

xgJdTj- KA+(77-,^\^r)]^(^)|f,5),(l+O(M/0). (3.29) 

This contribution was discovered in Boer, Mulders llll4ll . Another contribution at sub- 
leading order comes from a transversely polarized gluon, see Eg. l3.11l and Fig. 13.6b . 
Using the relation p'^/{-p^+c+ie) - -I +{c+ie)/{-p^+c+ie) and a decomposition like 
Eg. 13 .201 that term gives 



^'-P^T J (2;r)4 ^ '^^ ^'^ 



A{k) 



-7a7 y 

■7 



2k- 
xigA-^i^) - gA^^,e ,^t))H^)\P,S)AI +0{MIQ)). (3.30) 



The contributions from Eg. l3.29l and Eg. l3.30l can be combined into 






/ll^'«- 



5|(A(0)/Tr'^ 

r 



A{k) 



-7aT7 V 

■7 



2k- 
xgjd,j-Gp'(T^-,e,^T)>P(0\P,S},(l +0{MIQ))+0{g\ (3.31) 



The above result will be generalized to all orders in g in the following paragraphs. 

The sum over all possible insertions on both sides of the cut yields in leading order 
the full gauge link. Some of the insertions, like 2Mv/^ , encountered in Eg. 13. 91 and 

2Mv/'^_ + _ in Eg. 13. Ill give a suppression in MjQ. In the sum over all insertions and 
considering the next-to-leading order in MIQ such terms appear only once, either on 
the left or on the right-hand-side of the cut. 

Consider the case that this term is on the left-hand-side and that the gluon from the 
distribution correlator, which is related to the suppression, is longitudinally polarized 
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Figure 3.6: Various diagrams contributing at next-to-leading order in M/Q. 



(similar to 2Mw^ , as in Eq. l3.9t . In that case the other gluon insertions in the dia- 

i -In- 
gram will contribute to the gauge link. Those gluons cannot have their vertices between 

the vertex of the A^ -gluon giving the suppression and the vertex of the virtual photon. 
The reason is that the suppressed term is proportional to y" which prevents a longitu- 
dinal link between it, and that the suppressed term vanishes when fields at light-cone 
infinity couple between the suppressed term and the virtual photon. The conclusion is 
that the suppressed term always couples directly next to the photon vertex and the other 
leading insertions are adjacent to the fragmentation correlator. This contributes to the 
hadronic tensor as 



2Mw'"' , (generalized) 



if 



X (P, 5 |^^(0)/A(A;)i:'^l(0, rfYiri- 



-l/i 



17- 



■2k p| + ie 



rgA^(rj)ilf(^)\P,S%. (3.32) 



By performing the integral over /7|, using the relation J di] id"li^^\0,ri )A^(ri) = 

J^d?]^ id^id"X.^^HO, If), and applying the identities of the previous chapter (£0. 12.701 
to work out the transverse derivative (note that (G^"(oo, rj^, rjj)) - 0), one finds 



2Mvi^'' , (generalized) 

r 

CO 

x(l+0(M/Q)). 



-£«- ^'(0-, «,-)£--. f^(Or, ^r)^A«(oo, ^+ ^j) 



rHo\p,s),y=e 

(3.33) 



If the suppression term comes from a transversely polarized gluon like 2Mw^^ + _ 
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in Eg. 13. Ill then one generally has (following the same arguments) 
2MW!^_ + _ (generalized) 



/^>' r-u 






-e 



•P^+ipi (1-0 



X {P, S |0'(O)/A(^)i:W(O, Ti-)irai 



-p{y 



-2k pt + ie 



rgA"j.(Tj)iJf(^)\P,S),, (3.34) 



which can be rewritten into 
2Mvv^''^+ _ (generalized) 



d'^ 



-p\y- 

J {2nf 
xg 



e'''^<P,5|iA(0)/Tr'^ 



A(^) 



2k- 



x(l+0(M/Q)). 
The sum of Ea. l333l and Ea. l3.35l gives 
2Mv/'' + 2Mwt''_ , _ (generalized) 

''~P\T ' ' 1 ' 

-yarJ' y 



(3.35) 



-f 



^e'"^{p,s\my"Tf 

(2;r)4 



A{k) 



2k- 



xg fd^ £(0,ij-)G^"(,j)£^-^\T^-,nHe\P,S%l.^^. (1 +0(M/Q)). (3.36) 

CO 

If the suppression is on the right-hand-side of the cut, one obtains the complex conjugate 
of this expression with /z and v interchanged. 

Integrated over p-, the matrix elements as above appear often in cross sections, 
leading to the following abbreviation for it 



«Jj;.(x,pr,P,5) 
X (P S 1^/0)^ Jd77-i:[^](0, 77-)G^"(77)£^- ^Un-,n4'i(aP,S\ 



(2;r)3 

r 



,.=^.=0. (3.37) 



Note that ^tiQ-- has only nonzero transverse components. Using the identities of the 
previous chapter. Eg. 12. 701 one can show that 



<l>f}J(x,pT,P,S) = ^^^^"^(x,pT,P,S)-p'^ <:>^^\x,pT,P,S), (3.38) 
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where 



<i>^^\.(x,pT,P,S) 



_rd^ 

J (27T 






e'PHP,SmOW^KO,niD"^^j(aP,S}cy=o . (3.39) 



p*=xP* 



At subleading order in MIQ one also encounters a contribution from the interaction 
of a transversely polarized gluon from the fragmentation correlator with the hard part 
(see Fig. l3.6t ). Such a contribution gives 



IMW" = r d^p d^k 6\p + q-k) 2M\/^ ^^^^^ + other diagrams, (3.40) 



<^r^l{kM)f 



li + h +m 



(p + k\)^ — m?- + ie 



■1 



8 d% Tr 



,x>,c 



Tr^ 



yf'Al^(k,k,)r^y, 



(3.41) 



This result can be generaUzed by considering gluons from the distribution correlator as 
well, yielding the gauge link in the distribution correlator. When applying the Feynman 
gauge for the elementary part, the gluon from the distribution correlator which couples 
to the gluon going to the fragmentation correlator contains an additional term. This 
term cancels exactly the term which arises when calculating the leading order transverse 
gauge link with the Feynman gauge for the hard part (see the previous subsection). The 
result in the applied light-cone gauge reads 



2^<I/™, (generalized) = ^Jd^^i Ti" ^W/ Tr^ (A«^(^,^i))/^r. 



(3.42) 



Because this result is obtained in the light-cone gauge with retarded boundary con- 
ditions, the fragmentation correlator appears to be non-gauge-invariant. Assuming 
gauge invariance of the expression for the cross section one finds that the correla- 
tor jdk^ d'^ki A1^(k,ki) corresponds to the correlator / ^A^g_\^" (z'\ kr, Pi,, S h)j y°, 
which is gauge invariant and defined as 

Af}J(z-\kT,Pi„Sh) 

x<Q|g /drx"^' ""(±00+ OGtXQ, r, m)-C'" HC, /7^<A('7)in;^x>ou,cL-=o (3.43) 

Also here one has the relation 

A^^\J{z-\kT,Ph,Sh) = Afl(z-\kT,Ph,Sh)-k"j A^^\z-\kT,P,„Sh), (3.44) 
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where 

^fl{z-\kT,Ph,Sh) 

3 ^ J (2;r)32£p^ J (27r)3 
X {D\£^'"*'^\OT,^T)-C^'-^\±^\e)iD-IH^)\Ph-,Px)ouui-=a ■ (3.45) 

Considering now all insertions on both sides of the cut one finds that the hadronic 
tensor including MIQ corrections now (finally) reads 13511 (for notational reasons the 
contribution from 7°[A^" (z"', fc/-, Ph,Sh)V"/' was put in (ju ^ v)*) 



IMW^"^ I A^pt d^kr d^ipr + qr - kj) 



Tr" 



^^^\xB,PT)rA^'\z-\kT)Y 



+ Tr^ 



.-^r<i>^-]"(xs,PT)r^'-\z-\kT) 

Q^/2 ^ 



- A" V',^r)/<l>W(xB,pr)r.-^/ + (j" ^ v)* 

"^ ^ eV2 

hadronic tensor for semi-inclusive DIS includingM/Q corrections 



This expression completes the descriptions of Mulders, Tangerman 111511 and Boer, 
Mulders ill4ll by including the transverse gauge link. Apart from diagrams which are 
connected to the equations of motion, all possible gluon-interactions between the corre- 
lators and the elementary part have been included. Although an infinite set of diagrams 
was calculated including M/Q corrections, the final result still looks remarkably sim- 
ple. It is basically expressed in two types of correlators, O and O^ic, and similar ones 
for fragmentation. For O one can simply plug in the parametrizations of the previous 
chapter. We could in principle parametrize 3)^-1 g as well but it will turn out that this 
correlator only appears in certain kind of traces. Using the equations of motion those 
traces can be rewritten in terms of the akeady defined functions of <i>. 

3.2.3 Some explicit cross sections and asymmetries 

In the previous subsection the hadronic tensor was derived including next-to-leading 
order corrections in M/Q. This tensor was expressed in correlators which have been 
parametrized in terms of distribution and fragmentation functions. Inserting these para- 
metrizations, the expressions for the longitudinal target-spin and beam-spin asymme- 
tries, as published in Ref. |36], will be given. Three other papers considered recently 
beam-spin asymmetries as well and motivated the publication of Ref. |36]. Using a 
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model calculation, Afanasev and Carlson estimated in Ref. 111611 the beam-spin asym- 
metry and compared their results to CLAS data. Yuan obtained in Ref. | llTl an ex- 
pression for the asymmetry and made an estimation by using a chiral quark model and 
a bag model. Metz and Schlegel claimed in Ref. illSll . wh ich a lso considers longitu- 
dinal target-spin asymmetries, that the calculation of Ref. fllw] is incomplete. They 
completed the model by including other diagrams which, as the authors point out them- 
selves, a re no t compatible with the parton model at order M/Q. In Ref. 1 36] the analysis 
of Yuan f 1 IT"] is completed by including quark-mass effects and the new function g-'-. 

Explicit leading order cross sections can easily be obtained by replacing the correla- 
tors (I)'^^ and A'"^ in Eg. 13 .461 with the explicit parametrizations of the previous chapter 
(see Eg. 12371 |Z651IT^IT^ . Note that our choice of light-like vectors, Eq.|Tll is 
easily related to the frame in which the correlators were defined, Eq. 12.471 12.751 At 
order M/ 2 the other correlators O'^tL and A't ] need to be included. They will be 
handled here by making a Fierz-decomposition which relies on the identity 

Tr^[AB] = aib\ + 02^2 + a"b^a + a1b4a + afbsaii^ 

where oi = ^ Ti-«[A], ^2 = ^ ^^[Ajsl a? = ^ Tr^[Ar^, 

fl" = \ Ti^[Aiy"ji], af = -^ Tr"[Acr"f^y^] and similarly for bt. 

^ 2 'V2 

Using the Fierz-decomposition one encounters O^-ic and Ag-ic only in particular 
combinations of traces which allow for a simplification by the use of the equations of 
motion. Showing only the argument connected to the polarization of the parent hadron, 
the correlators 0,5-1 ^ and Aq-\q are decomposed into 

<l.W/(5) = <l>W/(0) + <J/(5,) + ^"/(^r), (3.47) 

A[:j/(5,) = a[:1^"(0) + a[:1/(5„l) + h^^^JiSHT). (3.48) 

Using the relation [/0-m]i^ - 0, which straightforwardly gives [iiy+a^^Dy-m-)^]\l/ - 
[iyD^-iy^D'^+imo-''^+ie''^P'^'ya-y5iDp]i^ - 0, one finds for an unpolarized target (the 
functions depend on x and p^) 

1 P^ 

-Tr[oW^"(0)cr/] = iMx e - im fi+Mxh-—hf, (3.49) 

iTr[<l)W/(0)/cr/r5]=0, (3.50) 

iTr[oW/(0)r^]-iefTr[<l.W^^(0)/rV5]=/^?(x/-+/-^/7f+«g--/i), (3.51) 

and for unpolarized observed hadrons (where the functions depend on z and z^kj) 

i Tr[Al:]^"'(0)(r,-] = iM^E^-^ - imzD\^ + M,,H^-^ - ^^Hf'^ , (3-52) 
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iTr[A^:]^"(0)/cr/r5]=0, (3.53) 



\ T>-[Ai:l/(0)r-]+^4, Tr[AH/(0)/r-r5] 



m 



= ;t-(D^H+,- ™ ^^H_,-G^H_^£,H). (3.54) 



For longitudinal target-spin asymmetries one has the relations 

1 

2 



Tr[a)W^"(5L)cr/] = (3.55) 

1 p^ 

-Tr[<l)W/(5£)/(r/r5] = -m^L^iL + /M.Y5LeL + Mx^l/zl - i^^z./^^^, (3.56) 






lTr[cDW/(5,)y1 _ ^ «/^Tr[<J^ r5,)/r^r5] 



= -f^fpTpixSift - 'j^^Lhf^ + ixSigt - iSigiL)- (3.57) 

It turns out^ that these relations are sufficient to rewrite the appearing correlators O^ji^ 
and Afl-ip in terms of the distribution and fragmentation functions. 

The vectors in which the cross section is expressed need to be related to the vectors 
of the Cartesian basis as introduced in the previous chapter. The following relations can 
be found for the vectors P, q, and P/, {e± = (e, + e^)/ \1) 

Q Q 

q = —p^- F ^+' (3.58) 

V2 V2 

xbM^+OjM^IQ') Q + 0{M^IQ) ,,„. 

P = = e- + = e+, (3.59) 

eV2 XB^ 

z„Q + 0(Ml/Q) Ml-P,l+0{MllQ^) 
Ph = -p e_ + e+ + Pi,^. (3.60) 

V2 zh V2e 

Comparing these relations with Ea. l3.2l (and neglecting M^/Q^ corrections) one finds 

V2 
n+ = e+, n- ^e- + —z-Ph±, qri. = -Ph±/z. (3.61) 

ZhQ 

This gives the following relations for any two transverse vectors nir and nj 

niT-riT - mT±-nT±. (3.63) 



In the way the calculation is presented here, it is remarkable that the equations of motion contain suf- 
ficient information to rewrite the appearing correlators (t>" .„ and A" „ in terms of il>, il>" A, and A?. It 
has been pointed out by Qiu in Ref. 1 1 19] that one can avoid the quark-gluon-quark correlators by using the 
equations of motion at an earlier stage in the calculation. That could make the result less surprising. 
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So, by applying the Fierz-decomposition, using the equations of motion, expressing 
the result in the Cartesian basis, and using FORM 1 120], the hadronic tensor can be 
obtained. The tensor, also given in appendix B.AI reads explicitly 13611 



2MW""' = [2MW^" + IMW"/] [l+O [M'/Q-) + 0{as)) , 



(3.64) 



2MIV7 ~ 2z 



fdVd- 



% S\pt^ - ^i _ kr^) 



X< I 



2eYkr'l 

Q 



M 



--fxs^//^H + -f/,<-'--/iG"' 



j-[-] 



M,, 



Ml/ 



lefpTl [ m , , , M,, , , 



{'^'■■^•'*>'^-'-^'^''] 



+< 5z,^iz,£»' 



,[-] 



Pt±-S 



+/Si 



lel^e^^'h 



t t± "-Txp 



+iS 



Q 
Q 



M 

£,-LH 



—8itD\ 



8iL- 



lH. 



h^ nl-1 . 



M 



M 



ftiiD/ hf,D/' - Xb BlH, 



[-] 



M,, . , £l- 



xsgtD^' + -^hf, -hf,D' 



H 



2MW^s 



f «= Izjd'pr d^kr 6\pt^ '^' ^^^^ 



x-g^7/, /)';'+ 



-] , 8^±kTi_ ■ Ptx_ - kr^prl , _,_ 



MMi, 



hfH. 



-/iD-'+/i +xs—h^ 

z Mh z 



2eYk/l 

Q 

2eYpT"l 

Q 



-[-] 



U-] 



M,, . , /f [- 



XBf^D\-' + '-^hf'-^ + —I-hfH^^ 



2MMi, 



U' v}p 



+S 



+s 



kT±E± St±p + St±e± krj^p 
2Mi, 



'i^lL'^l 



h.rHf-' 



M 



Q 

le, e^ pt^p 

Q 



M 



— XBhLHf-^-—g,i,Hf^-^+g 



G^i 



[-] 



M, 



Ml, 



M '^ z MMi, "• ' ■'^ ' 



the hadronic tensor of semi-inclusive DIS 



where the approximation-signs indicate that transverse target polarization has only been 
taken into account at leading order and any polarizations in the final state have been 
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result of contraction with L„i, 









^p 



-LP 



ieYb'l 



ip 



4G' 

4g! 

4Gl 

4G: 

^^ 

4^ 

^\ 
46' 

y' 
46: 



(-l+y-//2) 

(1 - y)(a ■ Cxb ■ e^ — a ■ e,- b ■ e,) 

(- 1 -H }')(« ■ e^b ■ Cy + a ■ gy /? ■ e f) 

(2 - y) Vl -y a ■ e^ 

(y - 2) a/I -y « ■ e>' 

/leCy - y^/2)(a ■ e , b ■ e^ + a ■ e^ b ■ e j) 
^y Vl "^ « ■ ^>' 



Table 3.1; Various contractions given frame-independently (see also Ref. 111511 ') 



discarded. The distribution functions depend here on xb and p^ while the fragmentation 
functions have the arguments z and z^A:^. It should be noted that the hadronic tensor 
satisfies q^W^'' - qvW^^ - 0, expressing electromagnetic gauge invariance. 

In the above equation the hadronic tensor is expressed in terms of the various dis- 
tribution and fragmentation functions. In order to obtain the cross section, Eq. 12.231 
the hadronic tensor still needs to be contracted with the leptonic tensor L^y. Those 



contractions are facilitated by the use of table 13.11 After having made these con- 
tractions the result is not yet a simple product between the various distribution and 
fragmentation functions, but rather a convolution (J d^pj d^kj 6-^(pT ±-Pii±/z-kT ±) 
fiixB, Pj)Dj{z, z^k^)). The convolution can be made into a simple product by multiply- 
ing the cross sections with some factors Pi,j_ and subsequently integrate over d^P/,j^ . 
Note that the leptonic tensor is independent of Phj_. 

The unpolarized /"/u^-integrated cross section is given by (L|,y denotes the unpolar- 
ized part of L^v) 



fd^Ph^ 



EvuEy 



dV 



d^l' d^P/, 






sQ^ 



y^ 



zMxb)Di{z) 



x(l+0(M/Q)+0(as)). (3.65) 

unpolarized cross section for semi-inclusive DIS at leading order 



One of the single spin asymmetries we will study here is the beam-spin asymme- 
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0.04 - 



P,=0.038+/-0.005 




-0.02 

-0.04 I- ALu=Pi+P2Si"* 



Figure 3.7: Azimuthal beam-spin asymmetry measured by CLAS L121.1 . The applied 
cuts are e^ > 1 GeV^ W^ > 4 GeV^, and 0.5 < z < 0.8. 



try. This asymmetry has been measured by CLAS |121] and HERMES Il22r . The 
asymmetry will be defined as 



iLf/ 






(3.66) 



where Ae - I corresponds to a polarization pointing towards the target. For a general 
weighted azimuthal asymmetry the following definition is introduced 



Al';f=fd^Pi,^Pi,^-aA^,„y 



(3.67) 



Weighting the beam-spin asymmetry with Pi,j_ ■ e, and integrating over d^Ph± 
gives [3ff\ 





-2y^]\-y MM,, 


\"^ 7f, W-^Hd) ^h f,/--LH(l) j.„7.W-L[- 


1(1) 


- (l-y+y2/2)f,Di Q 


Im ■'' ' M ■'' XBZeH^ 




Mh ^ ^ Mh J 






x(\+0{MIQ) + 0(as)), (3.68) 




azimuthal beam-spin asymmetry for semi-inclusive DIS 



where the functions depend on xg or z. 

The asy mme try given a bove contains besides the contributions given in Mulders, 
Tangerman IllSll and Yuan llll7ll . two additional terms: the term proportional to hfDi 
and the term proportional to g''~Di. Presently, all six contributions are unknown, mak- 
ing this asymmetry unsuited for studying one function in particular. To extract in- 
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formation one should get a handle on some of the con tributions, either through phe- 
menological studies (e.g. see Efremov et al. 11231112411 ') or model calculations. For a 
review on the Collins functions, Hf, the reader is referred to Amrath, Bacchetta, and 
Metz il25ll and the references therein. Models of the Boer- Mulders fu nction, hf, have 
been studied by Gamberg, Goldstein, and Oganessyan 1 126L ll27lll28ll . Boer, Brodsky, 
and Hwang 1 129|, and Lu and Ma |130, 131,]. Hwang has attempted to co nstruct a 
model for g-'-, but encountered problems with factorization in the model I132ll . 

More easy to interpret is the asymmetry of the produced jet which can be obtained 
from Eg. I3.68l bv summing over all possible final-state hadrons and integrating over 
their phase space. Neglecting quark mass contributions the azimuthal asymmetry of the 
jet is directly proportional to g-'-. 



Mxb) 
beam-spin asymmetry for jet production in DIS 



•■' Q (l-3' + r/2) /i(xb) 



As explained in the previous chapter (see subsection l2.5.2> . the function g^ exists owing 
to the directional dependence of the gauge link which is also connected with the non- 
validity of the Lorentz invariance relations. This makes the asymmetry of increasing 
interest for the understanding of the theoretical description (are such functions really 
non-vanishing?) even though this function does not have a partonic interpretation (like 
all the other T-odd functions). One way to access the asymmetry for the jet-momentum 
is to make an extrapolation by first considering the leading hadron, then the sum of the 
leading and next-to-leading hadron, etc. 

Next we consider target-spin asymmetries which have 
been measured by HERMES [133, 134, 135] and COM- 
PASS 113611 . Considering these asymmetries one needs to / 
express the polarization in the lab-frame into the Cartesian q ■. _...■ 
basis, Eq. 12.31 In an experimental setup the target is usually 
polarized along or perpendicular to the electron beam. Both 

directions are combinations of e,, e,. and e^. If the target „. , „ .,, ^ ^. 

• J- „ , ■ J ■ ■ ■ ~ a A a Figure 3.8: Illustration 

IS longitudinally polarized then it has nonzero S l and St±- f^, i u f 

. . 11^^ or tne lab-trame. 

The c ontributio n of the latter was suggested by Oganessyan 

et al. f 13 /, 138] to access transversity via longitudinal target-spin asymmetries (see for 

related work Diehl, Sapeta |139| ). 

Defining to be the angle between q and 1 in the target rest frame (see also Fig. l3.8> . 
the following relation can be obtained 



^Vr^ + 0(MVe^). (3.70) 
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Polarizing the target along the beam leads for St± and 5/. to 



Sl^+1+0{m^/Q^), (3.71) 

|SrJ = ^Vr^ + 0(MVe^). (3.72) 



Defining the longitudinal target-spin asymmetry as 



C(2MW^r^,-2MW^r_,) 






where 5^ = 1 denotes a polarization against the beam-direction, one obtains ||3 



(3.73) 



"■UL' 
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+ (l-y)(x^zA,<-'">) 






-(1-..//2)(^x../-a)] 






x{l+0(MIQ) + 0(as)), (3.74) 






longitudinal target-spin azimuthal asymmetry in semi-inclusive DIS 



where the functions depend on xb or z- The functions f^ and G^*'', which were 
neglected in previous analyses, have been included. 

Compared to the beam-spin asymmetry there are more terms contributing to the 
asymmetry. The sizes of the contributions are at present unknown but the asymme- 
try is certainly not dominated by the contribution from transversity hi (see also Her- 
mes 1 140]). On the other hand, by using the different y-dependence one can in principle 
extract the contributions from the Sivers function, f^^: \ and transversity. For this ex- 
traction different collision energies are required (y « Q^Uxbs)). Experimental statistics 
could be improved by integrating over all other external variables. 



3.3 The Drell-Yan process 

The calculation of the hadronic tensor in Drell-Yan is very similar to the calculations 
we performed in the previous section. Only the most important steps will be listed. 
In the Drell-Yan process the set of light-like vectors ({«_,«+) with n_ ~ n+ and 
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«_ ■ n+ = 1) is chosen such that 



P^ = 



Pi^ 



^ Q 



XiM\ 

eV2 X1V2 

Q X2MI 
=n- + —n+, 

X2V2 eV2 



Q + OiM^IQ) Q+0{M^IQ) 

q — — n_ H — «+ + ^7. 



V2 



V2 



(3.75) 



Sudakov-decomposition for Drell- Yan 



In this process one can also resum the contributions of 
gluons coming from the distribution correlator and connect- 
ing with the hard part. The difference with respect to semi- 
inclusive DIS is that the quark-fragmentation correlator is in- 
terchanged with an antiquark-distribution correlator describ- 
ing an incoming parton. This results in a gauge link run- 
ning via minus infinity for the distribution correlator The 
first order calculation illustrates this difference in direction 
very clearly. Inserting a longitudinally polarized gluon from 
the distribution correlator, like Fig. 13.91 gives the following 
contribution to the hadronic tensor (arguments connected to 
parent hadrons will be suppressed in this section) 
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Figure 3.9: A 
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ing antiquark. 





2Mwt 



IMW" = r dV d^k 6^{p + k-q) [iMv/^.^ ^^ + other diagrams] , (3.76) 



<^Up,Pi)r<^(k){iy-ti)i 



-(Jt + i/i) + m 

n 

(k+pip--m?-+ie 



0.11) 



In order to perform the integral over p^ the denominator, {k + p-i)^ - irp- + ie, is again 
simplified into 2p\k^ + ie (eikonal approximation). The integral over pj takes now the 
form 



/ 



dp\ 



2/?|fe" -I- ie 



AliTl) = ^e(r - ^-)Atir]). 



(3.78) 



This illustrates that there is a contribution when 77 -^ -oa, but not when 77 ^ 00. The 
pole in p| is now on the other side of the real axis with respect to semi-inclusive DIS, 
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leading to a gauge link via minus infinity. One finds for 2Mv/^^ 1™ 

x{-ig) ^ dl^- A^{tj-,^t)HS)\P,,S,)A^+0{MIQ)) + 0{^a\ (3.79) 

—00 

which is exactly the first order expansion of a gauge link via minus infinity, and an 
additional term, ^a, which gets canceled similarly as in semi-inclusive DIS. 

We note here that via the coupling of gluons to an incoming parton we have derived 
a gauge link via minus infinity. This is similar to the gauge link for the fragmentation 
correlator which arose from coupling gluons to an incoming quark. 

The subleading corrections are also very similar*". The complete result reads 



IMW" = - I A^pt d^kr d^^ipr + h - ?7-)Tr° 



<l>"(^i,Pr)/'i'"te,^r)/ 



7a^r^Ul(xuPT)r^^-\x2,kT) 

x{l+0{M^/Q^)+0(as)), (3.80) 

where 't'J!_, was defined in Eg. 13. 371 and (the contraction over color indices has been 
made expUcit) 



<:>^'(x2,kT,P2,s2)^ I ;:_,: e 



-ik^ 



^ rd^|2_d£^ 
J (27r)3 

x{P2,s2\i/fa^)j:i'=J(o,e)hbmP2,s2)c\(-=o , (3.8I) 

O'* (X2, kr, P2,S2) = 3)M«(X2, kr, ^2,^2) + ^? OM(x2, kj, ^2, ^2), (3.82) 



(P' ..(JC2,KT,P2, ^2) = I ; — 



«P2,52|[/Z)|'«A/(^)]/it^(0,r)'Ay,i(0)l^2,52)c|r=0 • (3.83) 

jr=X2P2 



The hadronic tensor in EQ. I3.80l has been derived by employing the diagrammatic ap- 
proach. Effects from instantons are for instance not included. Such effects could be 

^The two terms at subleading order appear opposite in sign. This sign difference comes from having a 
quark or antiquark propagator. 
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relevant for the angular distribution of the lepton pair and have recently been studied by 
Boer, Brandenburg, Nachtmann, and Utermann 1 141']. 

In order to make a connection with the Cartesian basis the external momenta are 
expressed as (xiX2S - Q^ + 0(M^)) 



Q Q 

^ = — - e_ + — - e+, 

xiM\+0{M\lQ') Q + 0{M\IQ) 
P\ = = e- + = e+ + P\±, 



fiV2 



X1V2 



P2 = 



Q + OiMllQ) X2MI + 0(M\IQ^) 



xi 



V2 



eV2 



e^ + Pi^. 



(3.84) 
(3.85) 
(3.86) 



Comparing these relations with Eg. 13.751 one can derive that (neglecting M^ jQ^ cor- 
rections) 



qTi. 

fiV2 
n+ = e+ -, 

qT± - -2xiP\_i_ - -2x2P2±- 



(3.87) 



To express the hadronic tensor in the Cartesian basis the following identities can be 
used 



twist 2^ I'/'i/ ^ twist 2^ I y 1^1 J r^ 



1'!wL2(^i'm)-^1>LL2(xi,Pr.) 



'5wL2(^l'/'r^) 



2V2e 



A-A^, 



^,- ,0(^2, kr) = -^-^<^r' ,Jx2,kT) 

twist 2^ z» ./ / 2 twist 2^ ■^' -* -^ 



2V2e 



■0[m^IQ^), (3.88) 
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2 V2e "'""'' 



2V2e 



+ (9(MVe^). (3.89) 
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This yields for the hadronic tensor ||35| 



2MIV'"' = i fd^pr d-kr 6-(Pt± + h^ - grx)Tr° [cD"(xi,prx)/<I>'^(x2,fcrx)/] 

2 Vie 2 Vie 

+ ^r''1'l:l;(^i,Prx)/<I>'^(x2,fcrx) - ^/<l'[,:l"efe,fcrx)r''®'"'(x,,prx) 



+ (A<^vr] (i + c)(M2/e-) + o(Q'5)) 



(3.90) 



the hadronic tensor for Drell-Yan 
in the Cartesian basis including Ml Q corrections 



In the parametrizations of the correlators, which are «± -dependent, the ii-t are replaced 
by e± (that is where |„^^e^ stands for). The functions in the parametrizations depend on 
xi,pjj_ or X2,kjj_. Integrated and weighted hadronic tensors are given in appendix B.AI 
The above result for the hadronic tensor is similar to the leading order result of Ral- 
ston, Soper 1 44 1 . The only actual difference is the presence of gauge links in the parton 
distributions. The Drell-Yan process can be studied experimentally at RHIC 1142| and 
by the proposed PAX-experiment at GSI. This has recently generated several theoretical 
studies 1124. 141. 143. 144. 1451 . 



3.4 Semi-inclusive electron-positron annihilation 

For electron-positron annihilation the calculation is also very similar. By including the 
possible gluon insertions on the outgoing partons one obtains gauge links which are 
running via plus infinity. 

The calculation is set up by choosing the light-like vectors ({«+,«_) with «_ ~ n+ 
and ri- ■ n+ = 1) such that 



Pli2 = f"- + — F"+' 

Z2GV2 V2 

Ph, = — F"- + F"+' 

V2 zieV2 

q = n_ + n++qT. (3.91) 

V2 V2 

Sudakov-decomposition for electron-positron annihilation 
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The complete result for the hadronic tensor reads (suppressing arguments connected to 
parent hadrons) 



/ 



IMW" = 3 d>r d^/tr 8\pj + kj - qr) Tr' 



fi- 



A"(zr'./'r)/A"(z2''^r)/ 






+ (JJ^ v)* 



{l+0{M^/Q^)+0(as)), 



(3.92) 



where A' , ' was defined in Ea. l3.43l and 
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aM^(z-i, fer, Ph,Sh) = aM«(2-i, ^r, p,,Sh) + k"j AW(z-', fer, Ph,Sh), (3.94) 



A^^f^{z-\kT,Ph,Sh) 

34^J (2;r)32£pJ (27r)3 "^^' '^ 



(0 ,+oo )|P/„5/,)out,c 



X oui{Ph,Sh\£^'^'- "'(Or, /?r)X''^' "'(icx^-, /7-)/£)>,(?7)l^)cL-=o 






(3.95) 



The external momenta in the Cartesian basis read 

Q ^ Q 

Ml + O(MfjQ^) ziQ + OiMl/Q) 
Ph. = — ■ 7= e- + ;= e+. 



Z2GV2 

l2 



^ 



^ z,Q+0(Ml/Q) Ml+0(Ml/Q^) 

Ph, = ;= e_ + e+ + Ph.±, 



V2 



Zl 



eV2 



(3.96) 
(3.97) 
(3.98) 
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giving the relations (neglecting M^/Q^ corrections) 



Q 



n+ - e+ 



qr^ - -PhtJ./zi 



(3.99) 



Using the following identities to express the hadronic tensor in the Cartesian basis 

a[-j / -I \ A-A+ .[-] , -1 J^^f^- 

^wL,2(2i '/'r)=-^A;j^,2(^, ,pr)-^ 



-^w!st2(^^^/'rJ 






^wist2<-^l '^r±)-— |— ^„i,t2<^^l '^7-±) 



V2e 



+ <9(MVe^), (3.100) 



A;wL2(^2.^r) = ^A;„|^^,(z-,^,) 



= ^lwL2(^2'^fcrJ 



/^-/^4 



o(MVe'), 



(3.101) 



one finds 133 



2MW""' = 3 



+ [-^AW(zr',/^rx)/AW(Z2Urx)r'' - ^AW(zr',m)/AW(Z2Urx)r'' 
+ t|r''AW;(z7',Prx)/AW(z^',^rx) - f4^''^";(^r',^rx)r''AW(zr',p,J 

+ 0/<^v)*] (l+0(MVe')+C)(a5)) • (3.102) 

hadronic tensor for electron-positron annihilation 
including next-to-leading order in M/Q 



The integrated and weighted hadronic tensors are given in appendix B.AI 

Subleading order effects in electr on-po sitron annihilation have already been studied 
by Boer, Jakob, and Mulders in Ref. [1031 The only real diff'erence of the above result 
with Ref. 1 1031 is that the fragmentation functions contain here a fully closed gauge 
link. 



3.5 Fragmentation and universality 

In transverse momentum dependent correlators we encountered gauge links via plus or 
minus infinity. In subsection l2.6.2l we found that due to the possible interplay between 
the two mechanisms for T-odd effects (the gauge link and final-state interactions) , it is 
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Figure 3.10: Two typical diagrams which appear in model calculations and produce 
T-odd effects. Diagram (a) appears in electron-positron annihilation and diagram (b) 
appears in semi-inclusive DIS. The dashed-lines represent some kind of particle (for 
instance a scalar diquark). 



not possible to relate those correlators which have different link structures. However, 
if one of the mechanisms turns out to be absent or heavily suppressed, then relations 
between fragmentation functions having diffe rent gauge links can be derived. 

In 2002, an article written by Metz llOOll appeared in which universality of frag- 
mentation was argued. In that paper the effect of the gauge link is modeled through a 
set of diagrams in semi-inclusive DIS and electron-positron annihilation, see Fig. l3.10l 
Within the model it is found that the Collins function appears with the same sign in 
both processes. In order to achieve this result the momentum in the loop in Fig. l3.10l is 
not integrated over the full integration domain. This restriction is justified to make the 
model consistent with the parton model but also makes it unclear whether the univer- 
sality result still holds if this restriction is lifted. 

Th e res ult of Ref. Ik 100 1 was formalized in a subsequent paper by Collins and 
Metz llOlll in which it is claimed that the gauge link is taken into account. Having 
similar diagrams as in Fig. 13.101 the authors eikonalize the quark-propagator next to the 
exchanged gluon, meaning that its transverse momentum is assumed to be small. This 
procedure is similar to the restriction made in Ref. 1 100] because it effectively limits the 
integration domain over the momenta. When comparing the two processes, it was found 
in Ref. 1 101| that the pole of the eikonal propagator does not contribute^ in the cross 
section and consequently the difference between the processes vanishes'*. According to 
the authors this shows that fragmentation functions are universal in QCD. 

A first comment on both papers is that restricting the integration domain of the 
momenta in the loop makes the interesting result less rigorous. It is not clear whether 



'Gamberg, Goldstein, and Oganessyan find in Ref. [ 146] in a slighdy dift'erent formuladon that the eikonal 
pole does contribute in T-odd efi'ects, giving a sign dependent fragmentation function if only the gauge link 
is considered (no final-state interactions). This result has been challenged by Amrath, Bacchetta, and Metz 
in Ref. 1 125|. However, if the result of Ref. 1 146] holds then, together with the final-state interactions, this 
model would illustrate the universality problems with fragmentation. 

After a single momentum integration all poles fall on one side of the real axis. This allows one to close 
the contour on the other side of this axis and one finds that the difference between the processes vanishes. 
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the obtained result finds its origin in the analytical properties of Feynman diagrams or 
whether it results from the applied restriction. Secondly, it is not clear whether the 
presence of vertices, through which the production of hadrons from quarks is modeled 
(like the quark-hadron-diquark vertex), introduce additional analytical structures which 
might be the origin for the result. 

When comparing the results of Ref. llOlll with the analyses in this thesis one ob- 
serves several diff'erences. For instance, the eikonal pole turns out not to contribute in 
the analysis of Ref. flOlT while we found in subsection l3.2. ll that this pole leads to the 
gauge link. In addition, the diagrams in Fig 13. 101 which should give the eff'ect of the 
gauge link (although not proven in Ref. llOlll '). give the same sign for T-odd functions 
in both processes. This is in contrast with the analysis in subsection l2.6.2l where it was 
found that there should be a sign-difference if final-state interactions were to be absent. 
From that point of view the effect of the diagrams in Fig B.lOl should be considered as 
a final-state interaction and not as an eff'ect of the gauge link. 

The idea of not having a gauge link in the fragmentation functions is by itself in- 
teresting. Fragmentation functions are defined through matrix elements involving the 
vacuum which is invariant under translations. It might be possible to use that property 
to show that fragmentation functions are already gauge invariant without the presence 
of a gauge link. If a gauge link is not needed to establish gauge invariant fragmentation 
correlators then the derived gauge link in fragmentation functions for semi-inclusive 
DIS and electron-positron annihilation might not influence the expectation value of 
the matrix elements. If this turns out be true then this would make the fragmentation 
functions sign-independent universal which would also explain the results obtained in 
Ref. 1 101 1 . In addition, in that case the Lorentz-invariance relations for fragmentation 
functions might be valid as well (see als o the discussion in subsection l2.5.2> . 

Summarizing, the arguments in Ref. flOl] lead to an interesting result although sev- 
eral issues need clarification. Another point we discussed was that gauge links might 
not be relevant for defining gauge invariant fragmentation functions. This could allow 
for a nonperturbative proof of universality of fragmentation. Experimental and theoret- 
ical studies on universality of fragmentation are recommended. 



3.6 Deeply virtual Compton scattering 

Up to now we derived in this chapter the gauge links in the unintegrated correlators for 
various high-energy scattering processes. In all these processes interactions between the 
correlators and the hard diagram were included in the form of gluon-lines. These gluons 
coupled to external quarks (or fields) on which the equations of motion were applied to 
derive the gauge links. Studying a different kind of process in which gluons are coupled 
to internal instead of external partons is therefore interesting. The process studied in 
this section is deeply virtual Compton scattering (DVCS). Although the gauge link in 
this process has been derived in several ways and is perfectly known, the method which 
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P-Pi+Q 



p\\ P-q' \\p + A P-Pi\\ Pi^P+'l \\p + A 



(c) 



Figure 3.11: Various contributions in DVCS. 



was applied in previous sections seems not to be presented for DVCS in the literature. 
After some introductory remarks that calculation will be presented here to illustrate the 
consistency of the applied approach. Note that in this section the considered diagrams 
are amplitude diagrams. 

In the DVCS-process an electron scatters, via the exchange of a virtual photon, off 
a hadron, giving, besides the production of a physical photon, the hadron a momen- 
tum change A (with A^ ~ -M^). This process gaine d sig r iificant po pularity after Ji 
showed its connection with quark angular momentum 11471 Il48ill49ll . Together with 
the intrinsic spin of the quarks, this provides indirect access to the quark orbital angu- 
lar momentum as well. Besides the contributions of Ji considerable progress has been 
made by many others. For an ov erview and introduction the reader is referred to Goeke, 
Polyakov, Vanderhaeghen JTsoll . Diehl lHJlll . and Belitsky, Radyushkin 1 152]. 

In order to describe the DVCS-process the diagrammatic approach can be em- 
ployed. We will slightly deviate from the more common notations of DVCS in order to 
follow the notations of this thesis. Restricting to the lowest order in g, as, and ae.m., 
there are two contributions to the amplitude which are represented in Fig. 13. lib and 
Fig. 13. lib . They are given by 



-i(-ie)(iQfef 
iM = u(k',A')yyu(k,A)— '^/^-^ (q'Xiyri'ef.iq'), 



q 



J {p + qr + It 



7 



+ <D(p, P, P')fi / ^'^^. /] + 0{g\ 
ip - q'Y + le 

where <!>(/?, P, P') is now an off-forward correlator defined as 



(3.103) 



(3.104) 



<i>ijip,P,P 



J (27T 






e'''^out<i",5'|lA>(0)^,(^)|^,5)i„,c 



(3.105) 



In the above definition the spin-labels of the hadrons and color indices of the quark- 
fields have been suppressed. 
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Introducing a Sudakov-decomposition to facilitate the expansion in M/Q 

P = =n_ H =«+, q = — =«- =n+, (3.106) 

Q^/2 jcV2 V2 V2 

g^" = ^'^^ - <«: - «>i , e^" H e'^-'^''«+,,«_,„ (3. 107) 

where x - xg + O (m^ / Q^j, and restricting us to A^ ~ Q and A^^ M^, one obtains 

in leading order in M/Q (applying also a Fierz-decomposition) 



T"' - i'^Z. mn^ + ^F^ edJ [1 + ^(^/2)] + 0(g), (3.108) 

^F^;. » = j^P' ^^T^ (^"^(^^'^ ^ (O.r^C/'^)) , (3.109) 

Tl. inn. = / d/'^ ^._^.1a^_,., (^T^(^^) - (iK'Pip')) ^ (3.110) 



where 



^(P^) = ^ J ^e'"'^"out<i^'i«/'(0)rXr)mn,c, (3.111) 

^(y^^) = ^ J ^^e'P*^\u^{P'\my^y54'(n\P}in,c. (3.112) 

amplitude for DVCS in leading order in Ml Q 



Since all the correlators are on the light-cone the only other possible leading contribu- 
tions come from the resummation of A^-gluons. The above result is therefore complete 
in the light-cone gauge (to all o rders in g) at leading order in MIQ and corresponds to 
the result obtained by Ji in Ref. lll48ll . Assuming the expression for the cross section to 
be gauge invariant leads to a straight gauge link in a covariant gauge between the two 
quark-fields in the correlators F and F above. These correlators can be parametrized 
into functions (generalized parton distribution functions) which can be measured in 
experiments or be predicted by theory. 

Instead of assuming the cross section to be gauge invariant we will consider now 
a covariant gauge and derive the gauge link by including the longitudinally polarized 
gluons explicitly. The insertion of a single gluon, connecting the correlator and the hard 
part (see Fig. 13. lib ), contributes to the amplitude as 

"'''' 15^^ + other diagrams, 

I^+A + m 



rf" = Tpf^ ^jj^^ + other diagrams, (3.113) 



^-^^.inD^ = (o/dV dVi Tr^' [<i>Xip,PuP,P')ri-^ 



(p + qY — m? + ie 



Migy-t,)i ^ //-//. 1^^- ^ v]^ (3.114) 
(p - Pi + q) — m^ + le ^ 
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where 

_|_^e'"f e'>'^''-f'out<P'l 0-/O) Al(jj) if,i(0 |P)i„,,. (3.115) 

Analyzing the /?| -dependence of Eg. 13.1141 and making the parton model assumptions 
one finds a pole at p^ ^ p^ + q^ i^ 0. Calculating the integral over p^ by taking the 
residue one finds 



^) + ie 



Tl^ __ = (/) rd«+ ^g-'Vr (p'|j,(o)y(/) 'iy__ 

r 

X (-0 Jd77- A+(;7-)e'(^'+''*"'»(r)|P)i„,c (1 + 0(M/Q)) . (3.116) 

oo 

We see in the above expression that the first order expansion of the gauge link runs 
between oo and ^" instead of between and ^ . When evaluating the ;?^ -integral one 
obtains the Heaviside function ^(-t/ ). From the two conditions, oo > ?]> ^^ and 
> 7/ , it turns out that of the gauge link only the part between ^~ and remains. This 
result can be rewritten as follows 



X |^l^ut<i''l'A(0)r^(-0 fd//- A+(;7-).A(r)l^>in,c 



2n 

r 



r 



+ (/)67'out<P'l«A(0)/r5(-0 JdT?- A+(77-),A(r)l^>m,c 



x(l+0(M/Q)). (3.117) 

The equation above is a copy of Ea. 13. 109l including the first order gauge link expansion. 
Although the momentum of the gluon-line did not vanish (no pole at p^ = 0), the 
procedure of including longitudinally polarized gluons from the correlator as illustrated 
here still provides the gauge link. The presented calculation can be straightforwardly 
generalized to all orders in g. In leading order in M/Q one obtains Eg. B.llOl and the 
matrix elements in Eq. B.llll and Eq. l3.112l containing the path-ordered exponentials. 



86 3 Electromagnetic scattering processes at leading order in a^ 

3.7 Summary and conclusions 

The diagrammatic approach was ap plied to va rious scattering processes. Following the 
ideas of Boer, Mulders, Tangerman |114, 115], and Belitsky, Ji, Yuan |26], we assumed 
factorization for the correlators and produced the leading order in as cross section for 
semi-inclusive DIS, electron-positron annihilation, and Drell-Yan. In Ref. IllSll cross 
sections were obtained by choosing effectively the A^ + A'' - gauge, leaving only 
the A^-fields from the correlators to consider. That result, which includes subleading 
order corrections in M/Q, was studied in Ref. I114ll in more detail by considering the 
A = gauge. In that gauge longitudinally polarized gluons between the correlators 
and the elementary scattering diagram need to be considered as well. In Ref. 111411 
they were explicitly calculated to order g^ of which the result was generalized to all 
orders in g by using Ward identities. Afterwards it was pointed out in Ref. [26] that 
there is an additional contribution which was not identified before (the transverse gauge 
link). The theoretical description was here extended by including all leading order as 
corrections and including M/Q corrections. The obtained result is similar to the results 
of Ref. Ill 5.1 except that T-odd distribution functions are now also included and the 
approach is fully color gauge invariant. 

The hadronic tensors were expressed in gauge invariant correlators containing gauge 
links. These gauge links arise from including leading-order-ffs -interactions which con- 
sist of gluons interacting between correlators and the elementary scattering diagram. 
The paths of these gauge links run via plus infinity if the gluons are coupled to an out- 
going parton, and via minus infinity if the gluons are coupled to an incoming parton. 
Following the same approach we also considered the DVCS-process in which gluons 
couple to internal parton-lines. The path of the obtained gauge link for this process is 
just a straight line between the two quark-fields as also encountered in DIS. 

Using the expressions for the hadronic tensor in semi-inclusive DIS, some explicit 
asymmetries were derived. One of the more interesting asymmetries is the beam-spin 
asymmetry for jet-production in DIS, which is proportional to the distribution function 
g-'-. This function was discussed in chapter |2] and originates from the directional de- 
pendence of the gauge link. Whether or not such functions are really non-vanishing is 
at present an open question. Its measurement would contribute to the understanding of 
the theoretical description. 

In the second part of this chapter, it was pointed out that the interesting arguments 
for universality of fragmentation as given in Collins, Metz 1 1011 contain some issues 
which need clarification. Therefore, the difference in paths of the gauge links in the 
fragmentation functions appearing in semi-inclusive DIS and electron-positron anni- 
hilation could lead to different fragmentation functions for the two processes. This 
might create a problem for extracting transversity via the Collins effect. Alternative 
ways to access transversity, besides Drell-Yan, are for instance //(/' production in pp- 
scattering 11531 . and via interference fragmentation functions 1154. 155. 156. 157il (in 
which the nonlocality of the operators is light-like). 
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3.A Hadronic tensors 

Based on Ref. |!351] some hadronic tensors will be given. Contributions from antiquaries 
can be included by interchanging q with -q and /i with v. The following definitions 
were used (suppressing arguments of parent hadrons): 

^f"{x) = ^A^pT p'i ^^''\x,pt), ^fix) = - ^A^pT p^T OW(x,/^r), 

%(x) = J (fpT <i>^^^-'(x, pt), ^lix) = J d^pj ^f"{x, pt), 

AUz'') = [d^kr A^^^"(z-\kT), A«(z-') = f dhr A^^^"(z-\kT). (3.118) 

Semi-inclusive DIS 

The translation of Eq. 13.461 into the Cartesian basis is similar as was done for semi- 
inclusive electron-positron annihilation or Drell-Yan (see Ref. b35.] for details). The 
unintegrated hadronic tensor reads 

IMW' = fd^PT d^kr 5\pT^ - PhJz-kT^)Tr^{[<^{xB,PTA.)rA(z\\kT^)Y) 

+ ['^{xB,PTA.)r^A{z\\kT^)Y - <:>(xB,PTjr^f^A(z-\kTjy 

+ (ju^y)*)](l+6)(MVe') + <3(«5))L^^^. (3.119) 

The integrated hadronic tensor reads (employing Eq. l3.38l and Eq. l3.44> 
fd^Ph^ IMW ^z^Tx° [<t(xB)/A(z-')/ 



+ 



(fl2^/0«(xB)/A(z-') + fl2^/A«(z-')/<t(x5) 

+ {fi^vr)\[\+0[M^IQ'-)+0{as))[^^^^^. (3.120) 

The single weighted hadronic tensor reads 

Jd^P,,^ Pl^ IMW' = z' Tr« [ - (5W"(xb)/A(z-1)/ + 0(xb)/A[^1''(z-')/] 

x(l+0(M/Q) + 0(as))\ . (3.121) 
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The Drell-Yan process 

The integrated hadronic tensor reads (employing Ea. l3.38l and Ea. l3.82> 



[(fqTi. IMW" = - Tr" [0(jci)/€)(x2)/ 



^ eV2 eV2 

-— -=r o', (xo/ofe) - -r -— ^ 

2gV2 2 gV2 



+ -^^/«',, (^i)/i'fe) - ^/^^^^-'"(xO/Ofe) 



2eV2 '^ 2 gV2 

+ Ou ^ y)*)] (l + O(MVe') + C)(a5)) L^^^^. (3.122) 

The single weighted hadronic tensor reads 

x{\+0{MIQ) + 0(as))\ . (3.123) 

Semi-inclusive electron-positron annihilation 

The integrated hadronic tensor reads (employing Eg. 13. 44! and Ea. l3.94> 



f(fqT± IMW" = 3 Tr" [A(z2')/A(z7')/ 



+ (^yA«(z2 ')/A(zr') - -e^/A-(zr')/A(z2^) 
^eV2 eV2 



%/AW''(zr')/A(zr') - /22i^AW«(zr')/A(z2') 
e V2 2 V2 

(fi^vy)]{l +0{M^/Q^)+0(as))l^^^^. (3.124) 



The single weighted hadronic tensor reads 



Jd^qr^ q"r^ 2MW''' = 3Tr« [A(z2')/AW'^(zr')r'' - A[*(z2')/A(zr')/1 



x{l+0{M/Q) + 0{as))\ . (3.125) 
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Color gauge invariance 
in hard scattering 

processes 




In the previous chapter several processes were studied at tree-level (leading order as). 
The hard scale was set by an electromagnetic interaction involving two hadrons and, as- 
suming factorization, cross sections were expressed in nonlocal scale-dependent corre- 
lators. These correlators contain gauge links of which the path depends on the process. 
As first noted by Collins |25], this produces a sign-flip for T-odd distribution functions 
when comparing Drell-Yan with semi-inclusive DIS. 

Following the same approach in this chapter, scattering processes will be analyzed 
in which the hard scale is set by an QCD-interaction (besides the participating hadrons). 
We will find that the gauge link does not only depend on the process, but even within 
a process different gauge links appear; the gauge links will depend on the hard part or 
subprocess. How awkward this at first may seem, the procedure appears to be consis- 
tent. A prescription for deducing gauge links will be given together with some results 
for gluon-gluon correlators. 

The fact that the gauge Unk depends on the diagram and not only on the process 
immediately raises questions on factorization and universality. This particular topic 
will be discussed in the last part of this chapter. 
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4.1 Gauge links in tree-level diagrams 

In this section gauge links in QCD-scattering processes will be considered. Having 
an as -interaction in the hard part at lowest order in g, the number of ways in which 
the gluons can be inserted is richer than for the previously discussed electromagnetic 
processes. Not only will we find more complex gauge links than in the electromagnetic 
processes, but we will also observe that the path of the gauge link depends on the 
elementary scattering diagram within the process. 

In the next subsection a simple QCD-scattering process (quark-quark scattering) 
will be worked out. Using the results of the previous chapter it will be possible to 
derive the gauge links very quickly. The calculation will elaborate upon the presence 
of complex link structures in QCD. 

In the second and third subsection we will consider processes in which gluons as 
external partons are present in the form of jets. Although the same approach is followed, 
the calculation in the second subsection is technically more involved. We will obtain 
some results for Drell-Yan and semi-inclusive DIS in which an additional gluon-jet is 
being produced. 

The structure of the obtained results suggests a general prescription for deriving 
gauge links appearing in diagrams. This prescription will be given together with some 
examples in the last subsection. 

As a final remark, most of the calculations will be done at the amplitude level. We 
will often anticipate by only keeping terms which contribute to the cross section. The 
intermediate results will be the amplitude diagrams in which path-ordered exponentials 
appear. However, in order to absorb these path-ordered exponentials into the correlator, 
one needs to consider the calculation at the cross section level. It is then also essential 
to sum over the colors of the other external partons. 

4.1.1 Gauge links in quark-quark scattering 

The quark-quark scattering subprocess is relevant for hadron-hadron collisions and will 
be studied in this subsection. The gauge links are derived for the correlators which 
are attached to the elementary scattering subprocess (squared amplitude diagram). All 
external partons are assumed to be separated by large momentum differences. We shall 
start by considering gluon insertions from a single distribution correlator and derive the 
gauge link. Results for the other correlators will be given at the end of this subsection. 
Let us begin by considering the subprocess in Fig. l4.1h (the resulting gauge invari- 
ant correlator will be given in Ea. l4.6t . Including a longitudinally polarized gluon (A^ 
with momentum /?! ~ n+) from a distribution correlator (belonging to the quark enter- 
ing the graph from the bottom), that gluon can be coupled in four different places; to 
the incoming quark, to the outgoing quarks, and to the exchanged virtual gluon. The 
inserted gluon introduces an extra integral over pi which needs to be performed. An- 
alyzing the p|-dependence one observes two classes of poles in p^. The first class 
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P-Pi 



(a) 




P-Pi 




Figure 4. 1 : Figure (a) represents a squared Feynman diagram contributing to quark- 
quark scattering. In figure (b) and figure (c) a gluon is inserted from the correlator 
into the scattering diagram. In figure (b) the dash indicates a pole at p'^ ^ 0, while in 
figure (c) the dash indicates a pole at p| ^ 0. 



consists of poles at p^ x Q which arise from gluon insertions on external partons in 
which the introduced propagator (by the inserted gluon) goes on shell (see for instance 
Fig. 14. lb ). The second class contains the poles at p^ + corresponding to other internal 
parton-lines going on shell (see for instance Fig. 14. lb ). 

In the previous chapter we have seen that the poles at p| a; give a contribution 
to the gauge link but the meaning of the poles at p^ ^0 is unclear Fortunately, in 
the sum over all possible gluon insertions it turns out that the poles of the second class 
cancel each other. For instance, in the case of a single gluon insertion one can show the 
following cancellation by doing the calculation explicitly 






= 0, 



(4.1) 



where the dashes indicate the taken residue. This result can be generaUzed to all order 
insertions' and was noted in Ref. p37l|. 

The cancellation as above is observed in various processes suggesting it to hold for 
amplitudes in general. Unfortunately, such a proof does not seem to exist (not even for 
QED) and is not simple to establish. Therefore, two general arguments in favor of this 
cancellation will be presented, but it should be stressed that via an explicit calculation 
of the cases considered in this thesis the same results can be obtained. 



By choosing the momenta in a convenient way (use <I>(p, pp 
is not much more complicated to show than the equation above. 



Pl,...,Pn-\- Pn, Pn)), the general proof 
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The first argument is connected with a Ward identity, which states that a longi- 
tudinally polarized gluon (on shell) with nonzero momentum does not couple to an 
amplitude which has only physical external partons. Since the inserted gluon (with 
pt + 0) is coupled to all possible places, except for one particular interaction which 
is already absorbed in the correlator definition (see for example Fig. 13.1b ). such glu- 
ons should not be able to contribute in the sum. The second argument in favor is that 
in the present case (see for instance Fig. I4.lt ) the pole at pt ^i corresponds to the 
internal gluon-line being on shell. However, this configuration is physically not possi- 
ble, because the coupling between two on shell quarks (with unequal momenta) and an 
on shell gluon is forbidden by momentum conservation. This is also connected to the 
fact that all tree-level Feynman amplitudes are finite as long as the external partons are 
well separated. According to this argument the cancellation in the above equation takes 
place in a nonphysical area. This argument for a single gluon insertion can be extended 
to an arbitrary amount of insertions. 

Having only the poles at p^+ k. 0, the task to perform is to evaluate those poles which 
arise from inserting gluons on the external parton legs. The calculation of inserting lon- 
gitudinally polarized gluons to a single external parton leg is similar to the calculations 
performed in the previous chapter An explicit example will be shown below. By in- 
serting the gluons on a single quark-line and following the same steps as in the previous 
chapter, one can derive the gauge link to all orders. There is only one important dif- 
ference: the color matrices of the inserted vertices (of the inserted gluons) cannot be 
simply pulled into the considered correlator because of the presence of the virtual gluon 
inside the graph. Therefore, those color matrices are "standing" on the quark-line on 
which the gluons were inserted. Taking combinations of insertions on several external 
legs is relatively straightforward because the manipulations of the inserted gluons on a 
certain external leg can be performed without using information from the other external 
legs. As a result, the combination of insertions gives a product of gauge links of which 
their color matrices remain on the quark-lines on which the gluons were inserted. 

As an example, let us reconsider the subprocess given in Fig. 14. Ih (also given in 
Fig. l4.2h ). Summing over the colors of the incoming and outgoing partons, its contri- 
bution to the cross section is proportional to (where O was defined in Ea. l2.37> 



cTpig. ^a = K Tr"^ [tbta^(p)] Tr"^ [f^f J 

= — Tr*' [<l>(p)] (where K is some constant). (4.2) 

Using the machinery of the previous chapter^, the inserted longitudinally polarized glu- 
ons from the correlator <i> produce in the amplitude diagram (left-hand-side of the cut) 
the following gauge links: X.^^- ^ (007^") on the outgoing quarks and X.^^' ^(^^-co") on 
the incoming quarks (the coordinate ^ is related to the field i/r(^) in the correlator 0(/7)). 
Insertions on the right-hand-side of the cut yield similar gauge links. The result of the 



Issues related to the equations of motion will be discarded in this section. 
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vL(^T-cx!-) L(-cx!TO 



(a) 




Figure 4.2: Illustration of gauge links in a squared Feynman diagram contributing to 
quark-quark scattering. Figure (a) represents an ordinary squared Feynman diagram. In 
figure (b) the result of the gluon insertions is presented. It is shown on which quark- 
lines the color matrices of the indicated gauge links (L) are acting. The superscripts ^^, 
and ^T or O7, have been omitted. 



all order insertions is graphically illustrated in Fig. 14.2b in which it is shown that the 
color matrices of the gauge links are at this moment on the external parton-lines. The 
result of the insertions reads (the notation will be a bit sloppy; the gauge links appear 
under the ^ integral which is present in the definition of O) 






Using now the identity (A^ = 3) 



if)ij(f)ki - -6ii5kj - —SijSki, 



(4.4) 



one finds 



0"Fig. lO fc 



2^/-l 



■ Tr^ [(D(p)XW(0,r)X'°HO,r)] 



12 



TrC [O(p)i:t+l(0, D] Tr^ [-C^^KO, DJ I (4.5) 



where the transverse gauge links were added for completeness. This addition can be 
performed uniquely by assuming consistency (gauge invariance). Using Ea. l4.2l for the 
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Figure 4.3: Various quark-(anti)quark scattering subprocesses. 



normalization, the gauge invariant quark-quark correlator is obtained 1371 l4i 



^Fis.W7\.ii(x,PT,P,S) 



-f 



(27r)3 



■e-PUP^Sl^rjiO) 



: [^_£^^Ho,n^r^ -c^°Ho,n- j-c^^Ho,n-c^°\o,n]Map,s) 






(4.6) 



In the correlator for this diagram a combination of path-ordered exponentials appear. 
Such gauge links were not discussed before and although this new result looks rather 
awkward at first sight, it is actually invariant under color gauge transformations. It 
should also be noted that the appearance of these complex gauge link structures is not 
just a peculiarity of QCD, but also appears in the QED-version of Fig. 14.21 In that case 
the gauge link is i;t+l(0,r)-£'°'(0,r) which was obtained in Ref. L37.I . 

Using the above procedure one can derive the gauge links which appear in other 
quark-quark and quark-antiquark scattering diagrams, see Fig. 14.31 There is only one 
issue which still needs to be addressed. When gluons are inserted from several correla- 
tors at the same time, interactions between the different insertions need to be included. 
In the previous chapter we discussed processes in which two correlators were present. 
The interactions between the insertions from the two different correlators were simpli- 
fied by choosing a suitable gauge. In the present case where four correlators are present, 
this trick cannot be applied. However, the distribution gauge links as given here are ex- 
act when two quark-jet production is considered, giving confidence in the results in 
which the fragmentation process is included. The results for quark-quark scattering 
were presented in Ref. |40] and are given in table l4.1l Note that in the limit of ^7- -^ 0, 
all gauge links in table l4.1l reduce to a straight gauge link, XP^- ^ (0^,^^), between the 
two quark-fields, see also Fig. 12. lib . 



4. 1 Gauge links in tree-level diagrams 



95 



\ correlator 


quark 


quark 


diagram 


distribution 


fragmentation 


Fig.|4ji 


A^WTr'^XM _ L£M£l°] 


±£l-] Tr^ X[°]t _ iXt-lXt^lt 


Fig.|431i 


2_£[+] Xr*- i]l°l - -i][+]X'°^ 


3 £[-] Xr*^ £[□]! _ 5£[-]£[n]t 


Fig.|4.5b 


jL^WTrC^Mt + I^H 


^i:HTrC£M + |XW 


Fig.lO; 


|£WTrCxMt-ixH 


IXHTr^XM-iXW 


Fig.lOli 


IXWTrCxMt-iXH 


|£HTrC£M-ii:W 



Table 4.1: The gauge links in the correlators which connect the external parton legs 
appearing in Fig. 14.21 1431 The gauge links in the antiquark correlators in the figures are 
this case the Hermitean conjugates of the gauge links in the quark correlators. 



4.1.2 Gauge links in semi-inclusive DIS and Drell-Yan with an ad- 
ditional gluon in the final state 

Introduction 

We will consider semi-inclusive DIS and Drell-Yan in which an additional gluon is 
radiated with some transverse momentum (at least more than the hadronic scale). The 
additional gluon is assumed to be observed and therefore its momentum is not integrated 
over. Possible contributions to the processes are given in Fig. 14.41 

In order to obtain the gauge link in the distribution correlator all insertions must 
be taken. The poles al pt i^ Q cancel each other for all possible insertions which is 
compatible with the arguments presented in the previous subsection. Only the poles 
pf X 0, coming from insertions on the external quark and outgoing gluon, will be 
considered here in detail. 

Inserting gluons only on the outgoing quark in semi-inclusive DIS makes the cross 
section proportional to (see also previous subsection and Fig. 14.4b ) 



cr~Tf[th£^^\0,nta^(p)]S. 



ah. 



(4.7) 



where the gauge link appears under the ^ integral in <l>(p) (sloppy notation). The above 
expression is by itself not gauge invariant, the possibility of inserting gluons on the 
radiated gluon simply has to be included. From this point of view there can be no reason 
to neglect such contributions. The interesting point here is that one can simply argue 
what the gauge link should be. If the procedure is to be consistent then the only possible 
gauge link which can appear in the correlator is a gauge link via plus infinity. Gauge 
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Figure 4.4: Figure (a) and figure (b) represent contributions to respectively semi- 
inclusive DIS and Drell-Yan in which an additional jet is being produced. 



links via minus infinity or loops, as encountered in the previous subsection, cannot arise 
because the gluons are inserted on partons which are all outgoing (incoming partons are 
needed to produce such gauge links). When taking insertions on the radiated gluon into 
account, there is only one possible way in how those insertions can contribute to obtain 
the proper gauge link. Inserting gluons on the radiated gluon should lead us to the 
following replacements 



th -» 



i:"^-f*(0",co- 



)tbZ°''^\^-,Q-). 



(4.8) 



Together with the insertions on the outgoing quark this provides us with a gauge link 
via plus infinity in the quark-quark correlator lISTIl . Ea. l4.7l becomes 



= ^TrC[XW(0,r)1>]- 



(4.9) 



With simple arguments, the effect of inserting gluons on a radiated gluon in semi- 
inclusive DIS was derived. Using the same replacement rules, we can now also study 
the gauge link structure of Drell-Yan plus an additional outgoing gluon (see Fig. 14 .4b '). 
Using the same rules, one finds in that case for the gauge link 13711 



|£W(o,r)Tr^ [i:t°""(o,r)] - ^x'"ko,r)- 



(4.10) 



The result is here a combination of gauge links via plus and minus infinity. Although 
the result is not very appealing, it is gauge invariant. The appearance of gauge links via 
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plus and minus infinity appear because insertions were taken on incoming and outgoing 
partons. Based on the previous subsection, such a result was to be expected. 

The replacements in Ea. l4.8l were derived by assuming the procedure to be consis- 
tent with color gauge invariance. In the following it will be argued that the replacement 
rules are in fact the proper rules. 

The explicit calculation 

The derivation of gauge links will be presented which underlies the results of Ref . 13711 . 
The amplitude in semi-inclusive DIS and Drell-Yan in which an additional gluon is 
produced will be considered at order as , and the gauge link will be derived by includ- 
ing gluon insertions. Before we begin with the technical derivation we note that if the 
polarization vector of the outgoing gluon is replaced by its momentum, then the ampli- 
tude vanishes which is a result of a Ward identity (at this order in as there are no ghost 
contributions). 

For the coupling of an inserted gluon on the external gluon-line some work needs 
to be performed. Consider the outgoing gluon with momentum / for a given diagram 
in a particular amplitude (for instance the left-hand-side of the cut in Fig. 14.4b '). In 
order to indicate the presence of the outgoing gluon, the polarization vector eJ (I) will 
be introduced, where c is the "color charge" of the gluon (c e 1, . . . , 8), a is the color 
index (a e 1, . . . , 8), and/? is the usual Lorentz-index. Although eL (/) is proportional 
to 6ac (as e" ~ 6"), this property will not be used because it is preferred here to show 
explicitly the presence of this external outgoing gluon. Using these definitions, the 
considered diagram can be expressed as (/, j, k,l e 1, 2, 3, and the Dirac indices are not 
expUcitly shown) 

tij8je^:iJH0Mf;,,(p,l)ci>„(p), (4.11) 

in which M denotes the remainder of the expression for the diagram and f^' is a 3 x 3 
color matrix. Although the calculation will be performed at the amplitude level, the full 
correlator O is present in the above expression for economics of notation. We will not 
rely on the part of <1> which is in the conjugate amplitude. 

Inserting a gluon with momentum pi ~ n+ and polarization «+ on the radiated gluon 
replaces O by O^ and changes the above expression into (in the Feynman gauge) 



/ 



j4„ ,i ( Vga'a j-bsa 

■'(I- piY + ie 



X [g"-{l - 2p,f + gApy + IT - 21-8"^) ef^\l) Mf;„(p-pul) K.m^P'P^^ 

-Jd'P^'ij j~']T.: f 



X 



{-2l-g"P + (l- p,rn1 + nifl) e^»(0 M^p-puD <i>i,,i(p, Pi). (4.12) 



When referring to this equation, we will refer to the right-hand-side. The first term of 
Eg. 14. 12! does not change the original Lorentz structure of the diagram and will yield 



I 



■qSm ~ 
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the first order expansion of the gauge link. The second and third term have a less 
clear meaning. In the processes studied in this subsection the second term vanishes in 
the sum over all diagrams in the amplitude, and the third term does not couple to the 
Hermitean conjugate amplitude and vanishes therefore in the cross section. Both terms 
do therefore not contribute which is a result of a Ward identity which can be applied 
since I - pi ^0. 

In other processes like two-gluon production the sec- 
ond and third term do not vanish. As will be shown later, 
it turns out that these terms together with the ghost contri- 
butions exactly cancel the nonphysical polarizations of the P^?^Pi' _ 

radiated gluon in the Feynman gauge. This cancellation is p-pit i^ IP 

far from being trivial and it should be considered as a firm _^^ "^^^ 

consistency check. 

When continuing with the above expression, Eq. 14.121 
one should realize that there could be an additional p^- Figure 4.5: Example 
dependence in the rest of the diagram (in M), see for in- 
stance Fig. 14.51 Before closing the contours and evaluating the integrals by taking their 
residues, the powers of p^ in the denominator and in the numerator should be counted. 
Using that ft ~ // it is possible to show that the denominator is one order higher than 
the numerator which justifies the closing of the contour at infinity. Using this result the 
p|-dependence in the possible present internal quark propagator (in M) can be analyzed 
more closely. The p^ in the numerator can be discarded because that term removes the 
pole at pI - leaving the pole at p^ ^0 behind. Such poles do not contribute which 
can be checked explicitly. The /:>|-dependence in the denominator can be discarded as 
well because only the poles at p^ = will be considered. 

Continuing with the non-vanishing terms of Eq. 14.121 and performing the integra- 
tions leads to (using that if'"^"ti, - [fj, f^]) 

first term of Eq. j4.12> 

= 5/^;""(Ofca7M;;;,(/,,o 

oo 

X J (H? e"'^<^' S Wmi-ig) J d77-A:(;7-, 0, ^T)U^m S } 

r 

+ 8/^al^'(0(tats)ijMf;,,(p,l) 

X J ^e'PHP,SmO)(-ig) j d,]-A:(,]-,0,^T)M^mS}, (4.13) 



which is exactly the first order expansion of replacement stated in Eq. 14.81 The above 
equation expresses that two link operators are produced for insertions on outgoing glu- 
ons. In the following we will generalize this result to all possible insertions. 
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-pi 



Tot™" 

A 



=^ 



(b) 



p-pif 



As»g.i^S»^ 



^ 



(c) 



Figure 4.6: Various insertions in the squared amplitude. The discussed interference 
terms in figure (a) and figure (b) cancel the contribution of figure (c). For the sum of 
these three diagrams one remains with the non-interference contributions (the part of 
the inserted vertices consisting of metric tensors) of figure (a) and figure (b). 



When attaching more than a single gluon we have to consider the Lorentz structure 
of the inserted vertices a bit more closely (see the three terms in Ea. l4.12> . The inserted 
vertex which is the closest to the electromagnetic interaction will be called the first 
vertex (see for instance the left triple-gluon vertex in Fig. 14. 6h ). and so on. When 
inserting a number of gluons the second term of the first inserted vertex still vanishes 
by a Ward identity in the sum over all amplitude diagrams. The combination of the 
first term of the first vertex with the second term of the second vertex also vanishes 
because of a Ward identity, and so on. The same holds for the third term of the last 
vertex which does not couple to the conjugate amplitude. The conclusion is that the 
only non-vanishing contributions come from the part of the inserted vertices which are 
products of metric tensors and interference terms. These interference terms consist of 
contractions between terms like the second and third in Ea. l4.12l belonging to different 
insertions. 

Consider two adjacent insertions with momenta pi - p2 and p2 on a gluon-line 
with momentum I - p\, giving the additional gluon propagators: (I - pif' and {I - pif' 
(see Fig. l4.6h '). The interference term between those two insertions produce a term in 
the numerator proportional to (/ - pa)^ canceling one of the gluon propagators. This 
interference term is canceled if one includes the contributions from the insertion in 
which the considered vertices interchanged (see Fig. 14.6b ) and the insertion in which 
the two gluons couple into a four-gluon vertex-^ (see Fig. 14.6b ). When an interference 
term does not come from adjacent insertions, then this cancellation occurs in a slightly 
more complicated manner It will be assumed that similar cancellations also appear 
when more than two of those terms interfere. 

Attaching A^ gluons to the outgoing gluon and considering now only the metric 
tensors of the insertions (like the first term of Ea. l4.12> modifies in the Feynman gauge 



^ Up to now we have not considered insertions in which the gluons from the correlator cross. When 
allowing for N such insertions one should actually add for all diagrams a symmetry factor 1 jN. The diagram 
which has the four-gluon vertex counts in that case double and should be multiplied by two. 
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Eg. lTTTI into 



/(I 



^ (-2/-)(-o ^„„^ ^ j-2r)(-0_^.,,,e3 



(-2/ )(-0 ,^„ 



-f 



'^(l-piY + ie (l-p2)^ + ie '" (1-Pn)^ + ie 

X 5/ ^^"'"(0 m;};,(/7-p,, 1.-;:^,^_,,(/,, p,-p2, ..., PN-i-PN, Pn\ (4.14) 



where we already anticipated on leaving out terms which do not contribute to the cross 
section. Considering the color structure explicitly. Eg. 14. 141 reads (see Eg. 14. 161 for 
clarification) 



Eg. flm 



/d>, 



■ • ^^Pn 



1 



1 



-P{ + 'e 



-Pn + 'f 



A sum over all possible combinations of ta and A^ color matrices ti. with the 
index / on the right of ta going down (from left to right) and on the left of tg 

going up. The sign of each combination is (_ l)number of terms right of r„ 

^^P^KiK;^.a,,m<^P^^---^Pn). (4.15) 

The factor in the middle contains the following terms in the case of four gluon inser- 
tions, 

tl,tlJatiJij - ti.tjljl^ti^ + ti^tiJathU2 + ■ ■ ■ (4.16) 

Performing the integrations over p-, leads to a path-ordered product. Using the fol- 
lowing identity for non-commuting a's and b's^ which can be proven by using induction 



I d?7i I d772... I d77i^ 



b{T]\) . . . b(ri„,)a(riM) ■ ■ ■ a(?7m+i) 

+ all possible coordinate permutations such that the b's are 

path-ordered and the a's are anti-path-ordered 
c c c 

= I d77i I d772 . • . I dj]^ b(rii) . . . b(rjm) 

c c c 

X I d77„,+i I dri,„+2--- I dT]Na(riN)...a(ri^+i), 



(4.17) 



identity to obtain a product of ordered exponentials 



So [a(Xi),aix2)] # and [a{xi),hix2)] i= 0. 
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one obtains for the part of Eq. 14.151 which consists of m f-matrices left of ?„ and n t- 
matrices right of ta (with n + m - N), the following expression 

[ti, . . . t,JJi, . . . ?,,„„],.. gj e^'''(0 My,p, I) 

r r r 

xjdj]-^Al{T]N)JdT]-„_^AUr]N-i) . . J d%^^^AlJT],„+l)M^)\P,S}. 

(4.18) 

This expression represents the m* order expansion of X.^^'^ (^" , °° ) multiplied with the 
n* order expansion of -C^'^'^ (oo ,^ ). 

Including now all terms of Eq. 14.151 by summing the above expression over all 
possible values of n and m, and summing over all A^, one obtains 

^S''(o^/M«;,(p,/) 

X J^e'>^<P,5|0-,(O)[x^-f'(r,'x,-)f„Xfr-^'(<x,-,^-)j..^,(^)|P,5), (4.19) 

which is exactly the replacement rule as stated in Ea. l4.8l 

By taking all possible insertions on a radiated gluon the derivation of the gauge link 
was shown in considerable detail in the Feynman gauge. The color matrices of these 
gauge links are convoluted with the rest of the diagram and are not directly acting on 
the fields in the considered correlator (see the next subsection for how this expression 
gets handled further). The calculation for an incoming gluon is very similar and its 
result will be stated in the next subsection. 

4.1.3 Prescription for deducing gauge links in tree-level diagrams 

The previous subsections suggest a general procedure for evaluating gauge links ap- 
pearing in tree-level squared amplitude diagrams. Although not rigorously proven^, 
a prescription will be given together with some examples. We note that the derived 
gauge links appear under the f integral of the considered correlator (the notation is a bit 
sloppy). 

In order to determine the gauge link for a correlator (its parent hadron is mov- 
ing mainly in the n+-direction), which is connected to a specific elementary tree-level 
squared amplitude diagram having physical external partons, the following set of rules 
can be applied: 

'in all considered cases the prescription yielded the coiTect result. The structure of the calculations in the 
previous subsections indicates that a general proof by induction should be feasible. 
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1. Write down the contribution of this diagram to the cross section keeping only 
the color factors. Rewrite this expression into a product of the amplitude dia- 
gram (left-hand-side of the squared diagram) and the conjugate amplitude dia- 
gram (right-hand-side of the squared diagram). Color charges of external partons 
must be made explicit. For the considered correlator we simply use <i>ab, for the 
other incoming and outgoing quarks we apply the vectors (m .'^ and m ) indicat- 
ing their color charge / and component j and similarly for antiquarks (v and 
V . ). For the gluons we introduce f^"*'"'""' where a is the color charge and b is 
the color component. Any structure constant, f'"^ has to be expressed in terms 
of the basic color matrices, r'^ via the relation /"'"' = (-2/) Tr'^[tc[ta, tb]]- 

2. Taking all the insertions from a quark distribution or antiquark fragmentation 
correlator on the external partons in the amplitude diagram (except for the par- 
ton connected to the considered correlator), can be translated into the following 
replacements (outgoing refers to crossing the cut): 



J ji ^ 

an outgoing antiquark: v. — > X ^' (^", oo^) v ■ ' 



an outgoing quark: m. ' — > m . ' X - (oo ,^ ) 



. an outgoing gluon: ^>-™'' fj;> ^ Xf^' ^\^-, ex,-) ^'^"'"tjff xfj- ^\^-,r) 

• an incoming quark: u — > Xf.^' (^^, -oo") u. 

• an incoming antiquark: v — > v X- (-oo^,^") 

1 (*3),in^(/j) r^T, £^ / J — -\ (a),in.(b) fBr, £^ / - j — \ 

• an mcommg gluon: e^' t]-' ^ J:]^- (^ ,-oo )e^'' f^/X^^ (-oo ,^ ) 

The coordinate ^ is the argument of the parton-field in the quark distribution or 
antiquark fragmentation correlator If the considered correlator is quark frag- 
mentation or antiquark distribution correlator, the same rules apply but with ^ 
replaced by 0. See also the correlators given in section fT/^ and section IT?I for 
further clarification. 

3. For the insertions in the conjugate amplitude diagram one makes the following 
replacements when considering a quark distribution or antiquark fragmentation 
correlator (outgoing refers to crossing the cut): 

• an outgoing quark: u.'^X.-^ (0", oo^) u ■ ' 

• an outgoing antiquark: v.'' — > v X -' (o°^, 0") 

. an outgoing gluon: 4^->-™V(;' ^ X^^ ^'(0", «,-) 4^^-°X/ 4" ^'(-"' »") 

1 (*)t (*:)t rOr, f * / - n\ 

• an mcommg quark: M -^ u . \L. (-co ,0) 

• an incoming antiquark: v — > X-^ (0", -cx>^)v 
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• an incoming gluon: e]^ ' , ^ -^ik (0,-oo )ej^ 1]^^ Zj (-co,0 ) 

and replaces by ^ when considering a quark fragmentation or antiquark distri- 
bution correlator. 

4. Using that u]^ - 6rb, e], = Srh, etc., simplify the obtained expression. At this 
point it is essential to sum over the colors of the other external partons. In the 
final expression the transverse gauge link pieces are included which can be done 
uniquely. 

5. Divide the expression by the normalization which can be obtained by replacing 
allX'sby 1 and (D by 1/3. 

6. The result is the gauge link of the considered correlator. 

To illuminate this set of rules three examples with explicit results will be given. 

Example: Drell-Yan 

We will reconsider the gauge link in the distribution correlator in Drell-Yan. The parent 
hadron is assumed to be moving mainly in the n+-direction. The result of each step of 
the prescription is as follows: 

1 : (T = /TTr*' [<l)[v*'^^(A:)v''^"(A:)]], where /Tis some constant. In analogy with Dirac 

spinors one has in this trace [v*'^''(A;)v'-'^^' (A;)],) = v {k)v {k). 

3 : cr = KTx^ [OX"'^' ^"(0", -oo-)[v('->(fe)vW"'"(fe)]Xfr' ^\-oo-,^)\. 

where in the last step the transverse gauge link was included. 

5 : normalization is K. 

6 : the gauge link is X'"^(0, ^ ). 

Example: Drell-Yan and gluon-jet 

In this example the gauge link is calculated for a distribution correlator (with parent 
hadron along n+) in Drell-Yan with an extra gluon-jet in the final state, see for instance 
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Figure 4.7: A tree-level diagram contributing to three jet-production. 



Fig. 14.4b . The results read: 



1 : (T = KTr^ [Of,[vW(^)yW'(fe)]f,]ef ""'e^^'^ 



3 : o- = /TTr'^ [OX"^- ^ (0", oo-)thX:'^' ^ (00", 0")r'"' ^ (0 



■Or, re 



r"r, f ^ /-A- 



") 



X [yW(fe)vW'(^)]i:^- ^'(-cx,-,^-)£f- ^\^, <x,-)t„j:^r. r (oo-,^-)]eW™'e(^ 



.(s)out (.s)outt 



4 : cr = (ii:/2)Tr^[<DXW(0,r)] Tr'^[i:'°^"'' (0,^)1 - -^ Tr'^[<l)X"(0,r)]- 

6 

4 

5 : normalization is - K. 

3 

6 : the gauge link is -£^^\0,^)Tr'^[£^°^HO,^')] - -^"(O,^)- 



Example: Three jet-production in hadron-Iiadron collisions 

In this example we consider the gauge link for the distribution correlator, with parent 
hadron along n+, attached to the lower quark-lines in the diagram in Fig. 14.71 This 
diagram contributes to three jet-production in hadron-hadron collisions. The results 
read: 



1 : cr = //^TrC[(I)f„f,f,M^Mfrff,]/^-^^^ej''^™'"'"ef™'"l'e^-''"V(''>™'e«°"'4-''^'"^ 
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2 : cr = KTl^[^tatbtc[u^u]X.^'''^io<y--o<,-)tdZ^''-^\-oo-,oo-)[tf,tg]X.^^'^\c<,-,^-)] 

(/i)outt J/)outt i;)in (/!)out M)out (;)int 
^ ^/ ^g ^rf ^a ^c f i 

(/i)outt J /)outt i;)in (/?)out Ji)out ij)™"!" 
^ V « rf « <■ * 

4 : o- = ^ Tr*^ [$£"(0,^)] + j Tr^ [<l)i:W(0,r)-£^°'(0,r)] Tr*^ [X'^^'XO,^)] ■ 

5 : normalization is K/3. 

6 : gauge link is ^£M(^q^^-^ + ^i:W(0,r)X'°l(0,r)Tr^ [^[^^"'"(O,^)! ■ 

4.2 Relating correlators with different gauge links 

The gauge invariant correlators, appearing in the diagrams, will here be related to the 
correlators introduced in chapter|2 We will find that the T-odd functions, appearing in 
different diagrams, can in general differ by more than just a sign. 

The gauge link structure consists in general of several terms of which each term 
consists of a link (-C'^"^' or £,^'^) multiplied by some traces of gauge loops. After an inte- 
gration over the transverse momentum of a correlator these gauge links collapse on the 
light-cone, one remains with a simple straight gauge link. The transverse momentum 
integrated distribution and fragmentation functions are therefore universal^. 

To treat the first transverse moment of T-even distribution functions, one can use 
time-reversal to project out the T-even structure of the unintegrated correlator This 
enables one to write the T-even part of the correlator as the average of itself and the 
correlator having the time-reversed gauge link (oo <-> -oo, as done in Eq. l2.62> . Using 
the identities in Eg. 12.701 which also gives the relation (based on Tr*' f,, = 0) 

5" TrC[£M(0,r)l| =0, (4.20) 

one can show that the first transverse moment of the unintegrated correlator is link- 
independent. The first transverse moment of T-even distribution functions is therefore 
universal. 

The first transverse moment of T-odd distribution functions having some gauge link 
is treated analogously. The T-odd part of the correlator is projected out by taking the 

In a fragmentation correlator the gauge links appear in two different matrix elements. To show the 
universality of integrated fragmentation correlators with an arbitrary gauge link, one can choose the light- 
cone gauge. 
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difference of itself and the time-reversed one (see for example Eq. I2.63> . Using the 
identities in Eg. 12.701 and Ea. l4.20l one can show that the first transverse moment of a 
T-odd distribution correlator equals the T-odd correlator as defined in Ea. l2.63l up to 
a constant. Therefore one can relate the first transverse moment of T-odd distribution 
functions to the ones appearing in semi-inclusive DIS. 

For example, let us consider the first transverse moment of the T-odd part of a 
distribution correlator containing the gauge link: £! - |X'^'(0,^")Tr''[X'°^'(0,^")] - 
jXy^Cd,^^). Using time reversal the first transverse moment reads as follows 13711 

x(<P, S |^(0)[|Xt^l(0, D TrC[X'°^"'"(0, r)] - \t^'\^, nM>{o\p. s >c 

- <p,5i0'(o)[|i:H(o,r)Tr^[^'°^(o,r)] - lx.^^Mo,n]H^)\p,s},)Lo 

= Ji'r"'" (4.21) 

The above expression is equivalent to | times the first transverse moment of the T-odd 
correlator as defined in Eq. 12.631 The first transverse moment of the T-odd functions 



(like the Sivers function f^j ) is therefore |-times the T-odd function which is mea- 
sured in semi-inclusive DIS. 

For fragmentation functions the interplay of the two possible mechanisms for T-odd 
effects created the problem for comparing the fragmentation functions in semi-inclusive 
DIS with the functions in electron-positron annihilation. The first transverse moment 
of those fragmentation functions read 



^W'^C^-i)^ iy f_lZ£ f ^e'/'''^V<- vp, p„i,wm rOr,rm+ -^^+^ 



i^Tiz-') = ^ > : ( 7^^^^:^ J^e''^''^^'^ut<n,i'zl«A(0)i:"^-^"(0^±-+)|Q), 



3^J (2;r)32£pJ 2n 

r 



+ j:°^'^(+oo+,^^)iD" 



ilr(0\Ph,Px}out.c\„-=r=o- (4.22) 

i]T=i;T=0 



By making comparisons in the light-cone gauge, the first transverse moment of a frag- 
mentation correlator with some gauge link can be expressed in the transverse moments 
of the fragmentation functions appearing in semi-inclusive DIS and electron-positron 
annihilation. 

As an example, let us consider the correlator with the gauge link structure X.' = 
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Table 4.2: The comparison of functions having complex gauge links with the functions 
measured in semi-inclusive DIS and electron-positron annihilation. The factors in the 
table are obtained by joggling with color matrices. The factors are thus effectively a 
function of A^c. 

a: The factors for T-odd distribution functions. The first moment of a T-odd function 
appearing Fig. 14. 21 and Fig. l4.3l is the T-odd function of semi-inclusive DIS multiplied 
by the given factor. So <^p^°^'^]l = factor x OlT-^dd]^ 

b: The factors for T-even and T-odd fragmentation functions, the first transverse mo- 
ment of the functions in Fig. 14.21 and Fig. l4.3l can be written as a linear combination 
of the functions appearing in electron-positron annihilation and semi-inclusive DIS. So 
AF,g.(...):; = (At+l^ + AM^)/2 + factor x (Al+J^ - AH^)/2. 



-C'°^(0,^^)-C'^^(0,^^). When taking its first transverse moment one obtains 



""^ ^' ^ 3 4^ J (2;r)32£pJ 



In 



-e'/'ftfVz 



X outiPh, Px\4'a(0)-Ct ^"(0^ ^^)4} ^"(-'^^ '^^)l")c 

— oo 



x(Q| 






rOr, f 



+ C'('-"'-'-")^7. (-",r)'^" 



'/« 



^g{mPh,PxUu 



|'r=r=o- (4.23) 



By comparing this expression to A^"^' in the light-cone gauge one finds 

Af^(z-') = 2A[*(z-')-A[-]'^(z-i). 



(4.24) 



The above expression states that the functions with this particular gauge link are equiv- 
alent to twice the functions measured in electron-positron annihilation minus the func- 
tions which are measured in semi-inclusive DIS. That relation holds for T-even and 
T-odd fragmentation functions. 
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Using the method as given above, the functions appearing in quark-quark scatter- 
ing can be compared with the functions appearing in semi-inclusive DIS and electron- 
positron annihilation. The results are given in table l4.2l 

In this section it was shown how the integrated functions and the first transverse mo- 
ments can be compared with the functions in semi-inclusive DIS and electron-positron 
annihilation. For the higher transverse moments this comparison is not always possible. 
In those higher moments new matrix elements are encountered which do not appear in 
semi-inclusive DIS or electron-positron annihilation. 



4.3 Unitarity in two-gluon production 

The gauge link does not only depend on the process but also depends on the diagram, 
forming a potential danger for unitarity. If a theory is unitary then the cross section for 
the production of nonphysical boson polarizations is canceled by the cross section for 
the production of ghosts an d ant ighosts. The fact that QCD is an unitary theory is not 
trivially to see. In Ref. » 17lll58ll 't Hooft and Veltman showed that in the sum over all 
diagrams unitarity is maintained to all orders in QCD. Having different gauge links (and 
thus different functions) for each squared amplitude diagram, it is not clear whether the 
approach still obeys unitarity. In this section it will be argued at lowest order in as that 
including a gauge link in the correlators does not produce nonphysical polarizations in 
two-gluon production in quark-antiquark scattering. This may be considered as a firm 
consistency check. Before starting the gauge link calculation let us review, without the 
gauge link, the unitarity proof for this process to some extent. 

Consider the gluon-production amplitude of the process 





(4.25) 



When squaring this amplitude there are nonphysical polarizations contributing in the 
Feynman gauge. In the complete square those polarizations have either e(/i) being 
forwardly polarized (e"(0 ~ /") and eih) being backwardly polarized (e''^(/) ~ l", the 
bar denotes reversal of spatial components), or vice versa. Together they contribute to 
the cross section as 



|Eq. (|425}f = 2v{k)ktdu{p) 



f 



•dab 



(h + hf 



v{k)htd'u{p) 



f 



■d'ba it 



ih + hf 



(4.26) 
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Figure 4.8: Two elementary scattering squared amplitude diagrams and the gauge links 
corresponding to the correlators attached to the quarks below (parent hadron is assumed 
to be moving along the n+ -direction). The gauge links are derived by applying the 
prescription of the previous section. As shown in this section, the prescription gives the 
correct result when the gluons crossing the cut have a physical polarization. 



The ghost-production amplitude is 





(4.27) 



Taking the square of this amplitude exactly yields Ea. l4.26l with an opposite sign^. The 
cross section for the production of nonphysical polarizations is therefore canceled by 
the ghost and antighost contributions. One remains with purely transversely polarized 
outgoing gluons. This explicitly shows unitarity at this order in as . 

Including the gluon insertions for the gauge link one finds in general different gauge 
links for different diagrams (see for instance Fig. l4.8> . Having different gauge links for 
the diagrams and thereby also having different functions, it is a priory not clear whether 
such a cancellation between ghosts and nonphysical gluon polarizations still occurs. 
In the remainder of this section the rather technical argument will be briefly outlined 
that this is still the case. Important intermediate results are illustrated in Fig. 14.91 and 
Fig. 14.101 The reader who is not interested in the derivation of these results can skip the 
next paragraphs and can read the last paragraph of this section. 

When inserting a gluon with momentum /:>, on an outgoing gluon-line with mo- 
mentum / and carrying the vector index a, the inserted vertex has the following form 

'The opposite sign comes from the cut fermion loop. 
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-21- g"'" + (/ - pifn"^ + «f r, (4.28) 

where a' is contracted with the internal part of the amplitude. When considering multi- 
ple insertions, the terms of the inserted vertices, which are like the first term of Ea. lT^ 
give back the same elementary amplitude diagram with gauge links according to the 
replacement rule as given in subsection 14.1.31 step 2 (this has been shown in subsec- 
tion l4.1.2T i. Arising from multiple insertions, we also discussed in subsection l4. 1 .2l in- 
terference contributions between terms like the second and third term of Ea. l4~28l These 
interference terms are expected to cancel against insertions containing four-gluon ver- 
tices. Of all the insertions, this leaves us to discuss the vertices in which combinations 
between terms like the first term with the second or third term appear. Those combina- 
tions were canceled in semi-inclusive DIS and Drell-Yan in which an additional gluon 
was radiated, but they do not cancel in the present case. 

Consider A^ gluon insertions in which the vertices contain A^ - 1 terms are like the 
first term in Eq. 14.281 and one term like the third term in Eq. 14.281 This particular 
combination contributes as having the outgoing gluon being backwardly polarized (it 
becomes proportional to /i ■ e(/i)). Coupling this term to the conjugate amplitude leads 
to a forwardly polarized gluon on I2 in the amplitude (see the discussion at the begin- 
ning of this section). In the sum over the amplitude diagrams only one non- vanishing 
contribution remains. This contribution comes from the diagram with the two triple- 
gluon vertices (like the third term of Eq. I4.25t and is proportional to /"' . It gets con- 
tracted with the inserted vertex (n'^l"), making this particular contribution proportional 
to {-i)r/{-2rp^ + ie) when the inserted gluon propagator (corresponding to this ver- 
tex) is included. When comparing this with the first term of Ea. l4T2l one finds that the 
contribution here is opposite in sign and is divided by a factor of 2. Being opposite in 
sign this contribution is as if we would have coupled the gluon to an antighost (with a 
ghost on I2) giving 1 /2 times the A^* order expansion of the gauge link. 

When considering A^ gluon insertions with A^ - 2 terms like the first term of Eg. 14. 281 
and 2 terms like the third term of Eq. 14.281 one finds that this contribution can also 
be interpreted as coupling to an antighost multiplied with -(1/2)^ times the A^* or- 
der expansion of the gauge link. Considering now A^ insertions and summing over 
all possible third terms (number m) and first terms (total number of combinations is 
N\/{Nl - ml)/ml), one finds that this contributes as if we would have inserted gluons 
on an antighost but with the following factor in front of the A^* order expansion of the 
gauge link 

'^ (_l)m+l Nl (l\^ 

- 1 - - . (4.29) 



z 



_"' {N-m)\m\ \2, 

m=l ^ ' 

Summarizing, when inserting gluons on the external gluon-line carrying momentum li 
one encounters besides the discussed replacement rule several other terms. Some of 
these cancel when insertions with the four-gluon vertices are included. The others can 
be interpreted as insertions on an antighost-line (in the amplitude) with 1 
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Figure 4.9: Illustration of the calculation. The * indicates that of all the insertions on 
h and I2 only certain parts are included. Each of those parts consists of several (zero 
or more) terms like the first term of Eq. 14.281 and several terms (one or more) hke the 
second or third term of Eg. 14.281 



the A^* order expansion of the gauge link. The replacement rule is similar as for the 
gluons. 

Inserting A^ gluons on external ghost-lines and antighost-lines is relatively simple. 
One obtains the original amplitude diagram times the A^*-order expansion of the gauge 
link multiplied by ( 1 /2)^. The replacement rule is similar to the replacement rule as for 
insertions on external gluon-lines, but the factor (1 72)^^ at each order of the expansion 
prevents one forming a gauge link which would make the correlator gauge invariant. 
However, the insertions on an antighost-line together with the contributions of inser- 
tions on an external gluon-line, in which terms appear like the third of Eq. 14.281 can 
be recast into a ghost-antighost production diagram with a full gauge link ((1/2)'^H- 
Eq.E21= !)■ 

The contributions, arising from insertions with terms like the second term of 
Eq. 14.281 work out in very similar way. Those terms can be interpreted as coupling 
to a ghost-line. Together with the insertions on a ghost-line the result can be written as 
the production of ghost-antighost multiplied with a full gauge link. The result of the 
above calculation has been represented in Fig. 14.91 

Summarizing, when inserting gluons in the two-gluon production amplitude, one 
obtains the gauge link for each amplitude diagram plus some additional terms. Part of 
these additional terms (interferences between the second and the third term of Ea. R3^ 
should be canceled by four-gluon vertices as explained in subsection l4.1.2l The other 
terms together with the insertions on ghost-lines combine into full gauge links with an 
elementary ghost-antighost production-diagram. We remain now with the same set of 
elementary scattering diagrams in which gauge links according to the replacement rules 
appear Those gauge links are at this moment not absorbed in the correlator The result 
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Figure 4.10: The result of the insertions on /i and I2 on the amplitude level. In obtain- 
ing this result properties of the conjugate amplitude were used. In the diagrams it is 
indicated on which quark- lines the color matrices of the gauge links are standing. 

is illustrated in Fig. 14. 101 

Insertions in the conjugate part of the diagram yield the same result. Since the 
insertions do not modify the property that if one external gluon is forwardly polarized 
the other has to be backwardly polarized, one can even take combinations of insertions 
on the left-hand-side and on the right-hand-side of the cut and obtain the same result. 

To show unitarity is now relatively straightforward. On the amplitude level, one can 
show that if one gluon has a nonphysical polarization, then the other gluon must have 
a nonphysical polarization as well in order to have a non-vanishing amplitude. This 
contribution, which comes from the diagram with the triple-gluon vertex, is canceled 
in the cross section by the ghost-antighost contributions, because the gauge link for 
ghost-antighost production is similar to the gauge link of that contribution. This shows 
for this example that the procedure of obtaining gauge Unks is consistent. 



4.4 Gauge links in gluon-gluon correlators 

In the description of evolution and in the treatment of hadron-hadron collisions, one en- 
counters quark-quark correlators and gluon-gluon correlators. In this section the gauge 
links in gluon-gluon correlators will be studied. Defining transverse momentum depen- 
dent gluon distribution and fragmentation functions is a difficu lt task. In several arti- 
cles this topic has been studied (see for instance Collins, Soper IIIILI112I1 . Rodrigues, 
Mulders 1 159], Burkardt |90], and Ji, Ma, Yuan 1 160]), but a definition, including hard- 
part-dependent gauge links as obtained for quark distribution functions, has not yet 
been derived on the basis of the diagrammatic expansion. Using the same techniques 
as for the quark distribution correlator and assuming the theory to be consistent, it is 
possible to "derive" gluon distribution correlators. Although one would rather like to 
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show the consistency of the theory, the approach to be followed here provides at least 
some insight. 

Similarly as for the quark-quark correlator the photon- 
photon correlator can be derived using the LSZ-reduction 
formalism*^. In the Coulomb gauge a photon-photon corre- 



\ / 



lator without gauge link (appearing for instance in photon / \^ 

photon-distribution scattering, see Fig. 14. lit reads 

Figure 4.11: A con- 

<^°'^ (x n P S)= f^^-^-^e'f^ tribution to photon 

'^ ' ^' ' J (27r)3 photon-distribution 

x{P,S\A"(0)A^aP,S)\e-o- (4.30) scattering. The corre- 
'p*=yP* lator is connected to 

Although this correlator is not gauge invariant, the expres- *^ photons entermg 

sion for the cross section is actually already gauge invariant 

which can be seen as follows. Making a gauge transforma- ° P ' 

tion modifies Ea. l4.30l into 

<l)^(p, p,S)J^^-^^'PRP, S |(A"(0)+5''H(0))(A^(f)+c^H(f))|P, S )Lo ■ (4-31) 

J (^^) p+=xP+ 

When performing a partial integration the derivatives become proportional to the mo- 
mentum p. These terms, however, do not show up in the cross section because they do 
not couple in the sum over all diagrams (Ward identity). The expression for the cross 
section is thus shown to be gauge invariant. 

Being gauge invariant, the expression for the cross section can be compared to 
any other gauge invariant object. If they are the same in a certain gauge (for instance 
h(p)-A=Q), then they are the same in any gauge. This justifies the use of the correlator 

Eq. (E3U) -> r "^^ / f e"'^hs(p)h^(p){P, S IF'^mpyl^i^m S }Lo , (4-32) 

with n{p)^ = 0, n{p) ~ p, hip) ■ p - I. So by starting with a certain correlator 
connected to the cross section we were able to derive its gauge invariant form. 

In general processes other kind of correlators appear, also containing gauge links. 
In principle these gauge links can be obtained by following the same procedure as 
outlined in the subsection 14. 1 .31 for the insertions. The result one obtains is the A- 
fields together with the presence of gauge links. Now to show to which gauge invariant 
expressions the obtained correlators correspond is more difficult than the previously 
considered case. The main problem is that one picks up a contribution of the gauge 
link when doing a partial integration as performed in Eq. 14.311 Such contributions are 
similar to the gluonic pole matrix elements which we also encountered for the quark- 
quark correlators. Since there is no reason to neglect those terms, this forms at present 
the main obstacle for obtaining gauge invariant gluon-gluon correlators. 

Identifying A-fields witli part ons is a complicated procedure in the LSZ-formalism in an arbitrary gauge 
(see for instance Itzykson, Zuber II6III I However, in the Coulomb gauge it is relatively simple, motivating 
the choice for this gauge. 
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Figure 4.12: Two processes with gluon distribution functions attached to the gluons 
entering the graphs from the bottom. Figure (a) contributes to Higgs-production and 
figure (b) contributes to two jet-production in hadron-hadron scattering. 



So how can one proceed? Well, one could argue that although it is not completely 
clear how to obtain the correlator in terms of the F-fields, the longitudinal gauge link 
as obtained by the presented rules should still be correct. If the approach is consistent 
then these longitudinal links should be closed at infinity by adding transverse gauge 
links, which can be done uniquely. In addition, if the cross section is gauge invariant 
(as it should be), then one can compare this expression in the light-cone gauge with any 
other gauge invariant expression. Comparing with the same expression, but with the A- 
fields replaced with the F-fields provides then a gauge invariant definition. It provides 
a solution for the moment, but more research is definitely needed here. 

Using this approach one is able to obtain gauge invariant correlators. As an exam- 
ple, results will be given here for the two subprocesses illustrated in Fig. 14.121 Using 
the prescription of subsection 14. 1 .31 and factorizing the color matrices with which the 
gluons couple into the hard part, the gauge Unk for the diagram contributing to Higgs- 
production can be found to be^ (Fig. P^-.lJh t 



(D: 



a/} 

Fig. ITT^ fl 



drd^^,X 



^J (27r)3 

x{P,S\ Tr"^ F^f^(0)£}-\0, f)F^''{^) [6'\0, f)^ 



l^,5>cL=o ■ (4.33) 



All matrices in the above expression are 3x3 matrices. This expression is equivalent to 
the result of Ji, Ma, and Yuan who study in Ref. [160] this particular process as a first 
step to factorization of effects from intrinsic transverse momentum in hadron-hadron 
collisions. The equivalence can be seen by the following relation 2 Tr*" {ta£}^\Q,^^)ti,£}^'^'HO,^^)'\ - 
-Cj,^ (0, ^") in which the A/-fields on right-hand-side are contracted with the structure 



In this diagram there is a virtual loop present which we have not considered before. In the next section 
we will study gauge links and loop corrections and we will argue that poles in the loop from insertions can 
be discarded. 
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constants (if""') instead of the color matrices. In a sense this relation was obtained in 
subsection l4.1.2l 

Another example is the gauge link of a subprocess which contributes to two jet- 
production in hadron-hadron scattering (see Fig. 14.12b ). The gauge invariant gluon- 
gluon correlator is found to be 



^Fig.EHIi-^J (27r)3 
[3 



x{P,S 



Ti-c [F+"(^)[xW(o,^-)]"'"f +'^(o)xt^(o,r)] Tr*^ [xf°'(o,r)] 



■ Tr*^ [f^''{^)£^'^^\0,^-)F^'^(0)£}^HO,^-)] 



|P,5>cU=o • (4.34) 



In this new result various gauge links via plus and minus infinity appear It is not 
possible to rewrite this result in terms of gauge links in which only structure constants 
are used (as in the previous example). 

In Ref. I9Q1 Burkardt studied gluon distribution functions and suggested a sum rule 
for the Sivers quark and gluon distribution functions. It may be good to point out that 
those gluon distribution functions, containing gauge links via plus infinity, appear in 
semi-inclusive lepton-hadron scattering, like two-jet production. A contribution to that 
cross section would be Fig. 14.12b in which the incoming gluons from the top of the 
graph are replaced by incoming virtual photons. 

In this section a method for obtaining gauge invariant gluon-gluon correlators has 
been suggested although a significant amount of work remains to be done. Note that 
one still needs to look for observables sensitive to the path of the gauge link. Those 
observables should be sensitive to the intrinsic transverse momenta of the gluons. One 
particular observable will be discussed for hadron-hadron scattering in the next chapter. 



4.5 Factorization and universality 

Up to now processes have been described by using the diagrammatic approach in which 
correlators were attached to an infinite number of hard scattering diagrams. These cor- 
relators, like O, <i>A, and A^^i's as defined in section l24l contained a certain renormal- 
ization scale and it has been assumed that these correlators could be factorized, enabling 
one to calculate in principle its scale dependence process-independently. We found at 
leading order in as that in several semi-inclusive processes, this infinite set of hard 
scattering diagrams and correlators combined into a finite set of diagrams convoluted 
with a finite set of gauge invariant correlators (containing a gauge link). Although the 
starting point was purely process-independent, it turned out that the hard scattering part 
determined the path of the gauge links in the final result. In this section the validity of 
the applied factorized approach will be discussed. 
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Whether processes allow for a factorized description has been studied for several 
decades. Ellis, Georgi, Machacek, Politzer, and Ross studied the issue of factorization 
for semi- inclusive PIS and Drell-Yan 1 162, 163]. By applying methods of Libby and 
Sterman il64ll65ll . Collins and Sterman were able to show in Ref. 1 166] to all orders 
in as that semi-inclusive electron-positron annihilation is free of infra-red divergences 
(divergences appearing in separate diagrams when momenta of virtual or radiated par- 
tons vanish, I — > 0). Subsequently, Collins and Soper introduced in Ref. I IIILI112I1 
infra-red free factorization formulas for electron-positron annihilation with two almost 
back-to-back hadrons being observed. Those factorization formulas were constructed 
at high qj- (order -qj ~ q^ = Q^ » M^) and for fully integrated over qj-, where qj is 
the transverse momentum of the virtual photon with respect to the observed hadrons. A 

factorization theorem for low qr {q\ M^) was also proposed. In that theorem the 

infra-red part, coming from the vertex correction in which all momenta become soft, 
was factorized from the jets into a soft factor In the same paper also an attempt was 
made to describe Drell-Yan, but a factorization theorem could not be obtained due to 
the interplay of final-state and initial-state interactions. Possible problems due to this 
interplay were also noticed by Doria, Frenkel, and Taylor L167II . 

The study on Drell-Yan was con tinued in several papers amo ng which pap ers o f 
Collins, Qiu, Soper, and Sterman 156, 104, 168, 169, ncllllIlTl, and Bodwin ] l?!. 
In the end it was believed that the problems of initial and final-state interactions were 
under control, because the interactions between spectators were expected to be on a 
longer time scale and should therefore vanish by unitarity'*'. This yielded fairly well es- 
tablished factorization theorems for (77-integrated Drell-Yan and small (77--unintegrated 
Drell-Yan 111311 . In 1992 Collins included straightforwardly polarizations of participat- 
ing hadrons into the factorization theorems Il74ll and suggested in Ref. l23ll a factor- 
ization theorem for semi-inclusive DIS at small ^7. 

In 2002, Brodsky, Hwang, and Schmidt showed that final-state interactions between 
spectators lead to single spin asymmetries in semi-inclusive DIS ]24]. Based on this 
surprising result, Collins concluded that the studied interactions are actually on a shorter 
time-scale than it at first sight seems. He pointed out in Ref. 12511 that such interactions 
could give problems for transverse momentum dependent factorization theorems (at 
small qj) in Drell-Yan in which both initial and final-state interactions are present. 

The discussion on factorization and universality for transverse momentum depen- 
dent cross sections was rece ntly c ontinued in several papers. In 2004, Ji, Ma, and Yuan 
argued and claimed in Ref. 11751 1 176.1 to have shown all order factorization theorems 
at small qj for Drell-Yan and semi-inclusive DIS. A one loop calculation illustrated the 
factorization theorems, which was generalized to all orders by using power counting. 
Complications from combining final and initial-state interactions, which appear at two 
loops or higher (see for instance Fig. l4.13L were not explicitly addressed. Further, it 
was not discussed whether or not the fragmentation functions are process-independent. 

By showing explicit results, it was subsequently pointed out in Ref. L37- .38- .3911 
that there could be a problem for the transverse momentum dependent factorization 

'"Being on a longer time scale it was not expected that those interactions could influence the short time 
scale production of the virtual photon. 
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theorems at small qj for semi-inclusive DIS and Drell-Yan. The problem appears due to 
the presence of the gauge link which has a non-expected behavior in processe s in which 
there are initial and final QCD-states present (in the language of Ref. | lOll ll75L Il76ll 
this issue appears at two loops and higher). This will be elaborated upon in the next 
subsections. 

The study on un iversa lity and factorization was continued 
by Collins and Metz llOlll . Inspired by earlier work of Metz llOflll , 
the conclusion was drawn that the fragmentation process is uni- 
versal for which an additional argument was provided (which 
was discussed in section l33t . Universality of transverse mo- 
mentum dependent distribution and fragmentation functions 
was consequently claimed for semi-inclusive DIS, electron- 
positron annihilation, and Drell-Yan, but the issue of final and 
initial-state interactions in Drell-Yan was not further illumi- 

nated. P^S^'"^ ^.13: A 

combined initial and 

The mentioned references use various starting points for flnal-state interaction 
the discussion of factorization, which prevents straightforward -^^ Drell-Yan 
comparisons of their results. However, one important point, 

namely the coupling of almost collinear and longitudinally polarized gluons to various 
places in elementary diagrams yielding the gauge link, is obtained in essentially the 
same manner in all the approaches. Being similar, the results on this point should be 
equivalent. Another point in the discussion of factorization are the virtual corrections 
which have not been studied so far in this thesis. Therefore, we shall consider the vertex 
and self energy corrections in the next subsection. Together with the results on gauge 
links appearing in tree-level diagrams, the validity of factorization theorems for various 
processes will be discussed in the second subsection. 



4.5.1 Virtual corrections 

In this subsection gauge links will be discussed for diagrams in which virtual correc- 
tions appear. Intuitively one can already guess the outcome. Since virtual corrections 
do not modify the nature of the external particles (incoming or outgoing), one does 
not expect to find different gauge links for diagrams in which virtual corrections are in- 
cluded. In this subsection technical arguments for some specific cases in semi-inclusive 
DIS will be given to show that this idea indeed holds. It will be argued in the Feynman 
gauge that the first order expansion of the gauge link is the same as for the diagram in 
which the virtual correction is absent L38.1 . The conclusion drawn for semi-inclusive 
DIS also holds for Drell-Yan and semi-inclusive electron-positron annihilation. The 
first order expansion of the gauge link is sufficient to illustrate the complications re- 
lated to factorization. 
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Vertex correction 

We will begin by studying the problem in QED and then extend the result to QCD. 
There are three possible ways to insert a single photon with momentum /?i to the vertex 
correction diagram (showing the integral over pi explicitly). 



/ 



dVi 




Ip-Pi pi 



(4.35) 



where the lines at the bottom of the graph are connected to the relevant part of the 
correlator ^\{p, p\)- 

The momentum dependence p\ can be routed back to the correlator via the photon 
propagator or via the electron propagators. We will choose here the latter. Having this 
/7 1 -dependence, the propagators contain poles in p\ which can be evaluated by taking 
their residues. By introducing dashes for the taken residues, this pole calculation is 
made more explicit 




(4.36) 



where the possibility that two propagators going simultaneously on shell has been dis- 
carded. It is fairly straightforward to show that the bracketed terms vanish by doing 
the calculation explicitly. The last diagram contains the pole on the external parton and 
yields the first order expansion of the gauge link. 

This result will now be extended to QCD. The virtual photon and the inserted photon 
are replaced by gluons, and the electrons are replaced by quarks. The calculation is in 
QCD slightly different, because in the QCD-version of Ea. l4.36l the order of the color 
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matrices in the bracketed terms is different. This can be solved by using the freedom on 
how to route the /:> i -momentum through the loop, giving us the following identity 




(4.37) 



Substituting this identity in the QCD-version of Eg. 14361 all bracketed terms (present in 
Eq. 14.361 14.37> are canceled. We remain with the last diagram of the QCD-version of 
Ea. l4.36l (giving the first order of the gauge link), the last diagram of Eg. 14.371 and the 
diagram in which the gluon is inserted on the virtual gluon. The latter diagram contains 
the inserted vertex consisting of three terms: -21^ g"' " + if^ l" + {I - pi)"' 11% (similarly to 
Eq. l4.28> . The first term cancels the last diagram of Eq. R371 while the second and third 
terms are canceled by terms appearing in a similar way when treating the self energy 
corrections (see further below). The conclusion is that the first order of the gauge Unk 
for the vertex corrected diagram is similar to the uncorrected diagram. 



Self energy correction 



The calculation of the gauge link in a diagram in which the 
parton connected to the considered correlator has a self energy 
correction (see Fig. 14.141 1 is similar to the calculation of the 
vertex correction. It can be straightforwardly shown that the 
gauge link is the same as for the uncorrected diagram. The 
calculation for the diagram, in which the self energy correction 
is on the other external parton, is conceptually more difficult. 
That calculation will be presented here in more detail. We will 
begin with QED and then extend the result to QCD. 

WTien taking a self energy diagram which is not directly 
connected to the considered correlator there are three possible 
places to insert a photon for the gauge link. This gives the 



Figure 4.14: Self 
energy correction on 
the parton belonging 
to the correlator. 
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following in terms of their residues 



j^r. 




(4.38) 



The term between brackets on the right-hand-side again cancels. The first of the three 
remaining terms can be calculated explicitly and yields the first order gauge link mul- 
tiplied with i-)6Z2 {6Z2 is of order g-). The second term is more difficult to calculate, 
because it contains a double pole. To circumvent this double pole problem we wiU 
exponentiate this diagram together with the last diagram. 

Consider A^ (including N = 0) consecutive self energy loops on an electron-line. 
Inserting a single photon either before or after each loop and summing over A^, the 
result can be written as product between two geometrical series of which one depends 
on pi and yields the first order gauge link expansion multiplied by Z2, and the other 
just gives Z2. Together, this gives the first order gauge link expansion multiplied with 
Zj ~ 1 -t- 26Z2. Note that in this sum also the first order link diagram without self energy 
correction is included. The total result of all insertions to order g^ is the first order link 
expansion (order g) times 6Z2 (order g^). 

In QCD one has to route the momentum partly through the virtual gluon as done 
when treating the vertex correction. Similarly as for the vertex correction, one obtains 
the first order gauge link expansion for the sum over the insertions, but some terms 
remain which come from the insertion on the virtual gluon of the self energy correction. 
Similar results are achieved when the self energy correction is on the other parton leg 
(see Fig. l4.14l i. Together, the remaining terms of the self energy corrections cancel the 
remaining term produced by the vertex correction. 

The conclusion is that the first order of the gauge link remains unchanged when in- 
cluding virtual corrections in semi-inclusive DIS. The same conclusions can be reached 
for semi-inclusive electron-positron annihilation and Drell-Yan. 



4.5.2 Evolution, factorization, and universality 

Scaling violations arise in a natural way when including higher order corrections in as . 
In this subsection as corrections in combination with gauge links will be discussed. 
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Figure 4.15: Various elementary scattering diagrams convoluted with gauge invariant 
correlators contributing to electron-positron annihilation. All corrections have correla- 
tors with gauge links via plus infinity. The derivation of those gauge links in similar to 
semi-inclusive DIS and Drell-Yan. 



One of the corrections, which was not discussed, is the self energy correction of the 
inserted gluons (or gauge link). Since virtual corrections are not expected to modify the 
gauge link structure of the graph, it will be assumed that its result can be rewri tten as a 
charge renormalization (as, see for related work Collins, Soper, Sterman lIlQSll '). 

Semi-inclusive electron-positron annihilation 

To compare the process at two different scales Q^, corrections in as should be included. 
At order as there are various contributions of which some typical examples are given 
in Fig. 14.151 In the figures the gauge links are indicated which result from all gluon 
insertions in the elementary scattering diagram. The gauge links are in this case all 
the same. When constructing the correlators, gluon radiation was absorbed up to some 
scale we will call ju^ ~ M^. When increasing the scale of the process, additional gluons 
will be radiated having a transverse momentum of at least fi^. Since we are considering 
small qr their transverse momenta cannot be to high (they should still be in the order 
of the hadronic scale). Although the gauge link was determined in subsection 14. 1.21 
for gluons having a large transverse momentum, the technical derivation for the gauge 
link is expected to be the same in the limit of -Ij — » ju^. That conclusion can also be 
reached by following the general arguments as given in subsection l4.1.2l 

When integrating over qr and over the unobserved gluon (with momentum I) one 
finds at this order that the infra-red divergences'^ are canceled i n the cross section, 
which is consistent with the result of Collins and Soper r ill[lll2il . When considering 
finite small qr and integrating the unobserved gluon over a restricted region Ij < A^, 
where A^ is of some hadronic scale, the cancellation of infra-red divergences can also 
be found. It may be good to point out that if the radiated gluon is integrated over a 
restricted part of its full phase-space, the bilocal operators in the correlators are still off 



' The infra-red divergences appear wlien / — ► in separate (liard) scattering diagrams which are convoluted 
with correlators. Following the general arguments of subsection l4.1.2l the gauge links depend on the nature 
of the external particles, giving every correlator a gauge link via plus infinity. 
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the light-cone and transverse momentum dependent effects (like T-odd effects) are still 
included. The factorized approach appears to be consistent. When constructing factor- 
ized correlators appearing in cross sections for small qj, no problems are encountered 
in perturbation theory at this order in as ■ Since soft gluon radiation does not modify 
the original gauge links, ladder diagrams (like the last diagram in Fig. 14.151 1 can be 
included when constructing the correlators, giving them a scale-dependence. 

Drell-Yan 

The situation is quite different for Drell-Yan. Considering as corrections in a similar 
way (see Fig. 14. 16^ . one encounters various gauge links in the correlators for the correc- 
tions. Gluon insertions on the radiated gluon were included to obtain gauge invariant 
correlators. If the approach is to be consistent (gauge invariant) then those contributions 
cannot be neglected or circumvented unless there are additional analytical properties of 
the correlator'^. Note that when the transverse momentum of the radiated gluon is of 
a hadronic size, the interactions appear at a similar time-scale as the insertions on the 
incoming antiquark in Drell-Yan or as the insertions on the outgoing quark in ordinary 
DIS. 

The above results point to difficulties when calculating the scale-dependence of 
the overall process. When considering as corrections to compare different scales, the 
behavior of the correlators is different. Since brehmsstralung diagrams have different 
correlators then the tree-level diagram, the effect of this kind of radiation cannot be 
absorbed in a correlator or other constructed objects (like for instance a soft factor) in 
the approach followed here. The scale-dependence is thus significantly more difficult to 
calculate than in electron-positron annihilation. It should be noted that when simplify- 
ing this process to QED, the problem of the gauge links in combination with radiation 
does not appear Since the photon does not carry any charge, the gauge links do not 
change when photon radiation is included. QCD has a different behavior here. 

In the previous chapter, section 13.21 the equations of motion were applied in the 
calculation of gauge links. Considering the diagrams, it seems that when increasing 
the scale of the process, which produces more gluon radiation, the gauge links in the 
correlators get modified. If radiation effects can be absorbed into the correlators ap- 
pearing in the factorized form <l> i8> O i8> //, then it seems as if the equations of motion 
were applied in the wrong way, or that the equations of motion are not invariant un- 
der scale transformations'^. Note that this problem does not appear when discussing 
electron-positron annihilation. In the following it will be indicated how this problem 
can be circumvented. 



For instance if T-odd distribution functions are zero or if gauge links are not an intrinsic property of tlie 
nucleon, but a vacuum effect instead. 

The same problem also appears in semi-inclusive DIS (see Fig. 13.21 . When applying the equations of 
motion we assumed analytical properties of the gluon connecting the fragmentation correlator (it was assumed 
to be outgoing). This assumption could be invalid and may need improvement. 
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Figure 4.16: Various elementary scattering diagrams convoluted with gauge invariant 
correlators contributing in Drell-Yan. At this order in as , all virtual corrections have 
correlators with gauge links via minus infinity, while all real corrections have correla- 
tors with the gauge link £' = l^l+l Tr*^ £^°^' - ^£^-\ 



We will try to construct the factorized correlator connected to the external par- 
tons which enter the process from above in perturbation theory, although this is strictly 
speaking not allowed. Constructing the correlator, the incoming hadron and its remnant 
is expanded in free parton states. Note that only the connected part of the diagrams 
needs to be considered. Since those partons are then essentially free, the gauge link in 
the lower correlator can be determined by using the results of the previous sections. For 
the subprocess in Drell-Yan in which partons, which cannot be absorbed in the lower 
correlator, do not cross the cut, one finds the gauge link for the lower correlator to run 
via minus infinity. However, if one of the partons, connected to the upper hadron, ra- 
diates a gluon to the final state (note that its energy depends of the process), then the 
gauge link in the lower correlator gets modified. So for each component of the Fock- 
state expansion of the upper correlator, one can calculate the gauge link for the lower 
correlator However, each Fock-state component will have in general different gauge 
links, making the procedure unsuitable for constructing the upper correlator factorized 
from the rest of the process. This problem only appears for cross sections which are 
sensitive to intrinsic transverse momentum. After an integration over qr, one finds at 
leading order in M/Q that all gauge links in the correlators are on the light-cone and 
run along straight paths between the two quark-fields (see Fig. 12. lib ). In that case the 
problems with gauge links and factorization disappear 

The above arguments illustrate that even when applying perturbation theory and 
a simple Fock-state expansion one encounters difficulties with factorization theorems 
sensitive to the intrinsic transverse momentum. This problem originates from gluons 
which are inserted on both initial and final-s tate partons. Since this issue has not been 
explicitly discussed in Ref. lll0lMl75Ml76ll . factorization for azimuthal asymmetries 
at small qj remains an open question for Drell-Yan. These problems are not encoun- 
tered when considering fully ^/--integrated factorization theorems at leading twist. In 
that case one finds that the infra-red divergences are canceled at this order in as . At 
subleading twist, the ^/--integrated factorization theorem contains correlators off the 
light-cone (see the previous chapter, appendix 13. A> and therefore it suffers from the 
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Figure 4.17: Various elementary scattering diagrams convoluted with gauge invariant 
correlators contributing in semi-inclusive DIS. At this order in as , the distribution cor- 
relator contains a gauge link via plus infinity. The virtual corrections have a fragmen- 
tation correlator with a gauge link via minus infinity, while all real corrections have a 
fragmentation correlator with the gauge link X.' = |X'^^-C'°^' - ^-C'^ '• 



same problems as for the transverse momentum dependent factorization theorem. 

Semi-inclusive DIS 

The situation of semi-inclusive DIS is similar to Drell-Yan (see also Fig. I4.17> . The 
real as corrections modify the original link structure of the fragmentation correlator. 
This forms a problem for factorization theorems which include eff'ects from intrinsic 
transverse momentum. The integrated cross section is still infra-red free at leading 
twist, but at subleading twist problems with factorization occur 

Concerning the construction of factorized correlators via a Fock-state expansion, 
when expanding the fragmentation correlator in free parton states one finds for each 
component the same gauge link in the distribution correlator This points out that for 
(multiple) jet-production in DIS factorization theorems sensitive to intrinsic transverse 
momentum eff'ects might be feasible. When expanding the distribution correlator in 
parton states, several partons will cross the cut which cannot be absorbed in the frag- 
mentation correlator. These partons influence the gauge link structure which forms a 
problem for constructing the distribution correlator in semi-inclusive DIS. 

It may be good to point out that if gluonic-pole matrix elements for fragmentation 
vanish, then the fragmentation process is univers al (univers ality of fragmentation func- 
tions has been advocated by Collins and Metz 110(1 llQlB '). In that case the problem 
with gauge links in the fragmentation correlator does not exist and factorization theo- 
rems could be feasible. ^_^ ^_^ 

As a final remark, in Ref. 111751 117611 factorization has been discussed including 
M/Q effects. However, since contributions from gluons at M/Q have been discarded 
(the Og-ic and Ag-i^ terms of chapter|3}, it is not clear whether the factorization theo- 
rems are complete as they stand. 
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4.6 Summary and conclusions 

In this chapter gauge links appearing in more complicated diagrams were considered. 
When deriving those gauge links by considering gluon interactions between the correla- 
tors and the hard part, two classes of poles were encountered. The first class contains the 
poles pt X Q which correspond to the limit in which the gluon momenta vanish. This 
class yields contributions to the gauge link. The second class contains the other poles 
and in the examples shown, those poles canceled. Supplementing general arguments 
in favor of this cancellation have been given and a set of rules has been conjectured 
that allows one to deduce gauge links for arbitrary diagrams. Applying these rules for 
gluon-gluon correlators, gauge invariant transverse momentum dependent gluon distri- 
bution functions were obtained. In a sense, these rules factorize interactions between a 
correlator and the elementary scattering part. 

It was found that gauge links do in general not only depend on the overall process 
but also depend on the elementary scattering diagram. The encountered gauge links also 
had more complicated paths than the ones appearing in semi-inclusive DIS or Drell- 
Yan. To compare the correlators with these different gauge links, we considered the 
first transverse moments of distribution and fragmentation functions. It was found that 
these transverse moments can be related to the transverse moments appearing in semi- 
inclusive DIS and electron-positron annihilation. For the higher transverse moments no 
relations were obtained. 

The fact that the gauge link in the correlator depends on the subprocess instead 
of the overall process could form a potential danger for the unitarity of the approach. 
However, in an explicit nontrivial example (two-gluon production) it was found that the 
factorized description in terms of gauge invariant correlators still respects unitarity. 

We also discussed the subject of factorization and universality for semi-inclusive 
DIS, electron-positron annihilation, and Drell-Yan. In semi-inclusive electron-positron 
annihilation we did not encounter any inconsistencies because QCD-states app eared 
only in the final state. For semi-inclusive DIS and Drell-Yan there was in Ref. l llOll 
Il75l ll761 recently important progress reported on small (77-unintegrated''* factoriza- 
tion theorems, but one particular issue of dealing with both final and initial-state in- 
teractions has not satisfactory been treated. This issue is related to the observation 
that gauge links in transverse momentum dependent correlators can get modified when 
considering radiation effects. This could create a problem for factorizing effects from 
intrinsic transverse momentum appearing in cross sections for semi-inclusive DIS and 
Drell-Yan. Since fully ^^-integrated factorization theorems contain at subleading twist 
effects from intrinsic transverse momentum (see appendix 13. A> . it forms at that order 
a potential problem as well. Factorization of intrinsic transverse momentum depen- 
dent effects in Drell-Yan and semi-inclusive DIS remains therefore at present an open 
issue P?*]. Similar problems also appear for semi-inclusive hadron-hadron scattering. 

The non-Abelian character plays an essential role in the discussion of final and 



q is the momentum of the virtual photon and —q- ~ M 
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initial-state interactions. For instance, in QED where photons do not carry charge, pho- 
ton radiation does not modify the gauge link structure. Although more complicated 
gauge links than the ones appearing in semi-inclusive DIS and Drell-Yan are still en- 
countered, factorization of effects from intrinsic transverse momentum in hard QED 
scattering processes should be feasible. 

It may be good to point out that in QCD the problems for semi-inclusive DIS orig- 
inate from the gauge link in the fragmentation correlator in combination with gluon 
radiation. If this gauge link turns out to be irrelevant (which means that gluonic-pole 
matrix elements for fragmentation functions would be suppressed), then the problems 
with factorizing effects from intrinsic transverse momentum disappear. In that scenario, 
which is favored by the in sectio n |331 differentlv given interpretation of the results ob- 
tained by Collins and Metz llOlll . factorization theorems including effects from intrinsic 
transverse momentum might be obtainable. Such factorization theorems can also be ex- 
pected to hold for (multiple) jet-production in DIS, because in those cases the problem 
with fragmentation functions is not present. 

As a last remark, in the diagrammatic approach used in this thesis, no assumptions 
were made for the correlators other than the assumptions similar to the parton model 
(see chapter|2l. Many issues might be unknown for that reason, one of them being the 
evolution equations for transverse momentum dependent functions (for a detailed study 
of this problem see Henneman |66]). It is possible that additional assumptions, which 
should be physically justified, might solve some of the issues. One possibility could be 
that gauge links in general do not influence the expectation values of matrix elements. 
In that case T-odd distribution functions would be zero and proofs on factorization and 
universality would be simplified significantly. Another possibility could be that gauge 
links correspond to certain interactions in the hard part in a particular kinematical limit 
(for instance the vertex correction in the limit that the gluon is collinear with one of the 
external quarks). In that case the gauge link is not an intrinsic property of the nucleon, 
allowing (part of) its contribution to be possibly absorbed in other matrix elements or 
factors. 
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Results for single spin 

asymmetries in 

hadronic scattering 



In the previous chapter the necessary tools were developed to determine the Wilson 
lines which appear in correlators in hard scattering processes. As discussed, it is not 
clear whether effects from intrinsic transverse momentum allow for a factorized de- 
scription in hadron-hadron scattering. We will assume it does and obtain results for 
single spin asymmetries. In these asymmetries the effect of the gauge link appears to be 
more than just a sign; it also determines the sizes of the asymmetries. It will be shown 
how the intrinsic transverse momenta of quarks can be accessed by observing unpolar- 
ized hadrons in opposite jets. For the ease of the calculation, only contributions from 
quark distribution and fragmentation functions will be considered. Gluon distribution 
and fragmentation functions are neglected, but can be straightforwardly incorporated 
by following the same procedure. 
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5.1 Introduction 

We will study the transverse momentum dependent distribution and fragmentation func- 
tions via single spin asymmetries in hadron-hadron scattering. Within the applied di- 
agrammatic approach, an odd number of T-odd functions is needed to produce cross 
sections for single spin asymmetries at leading order in as . At leading twist, T-odd 
functions are transverse momentum dependent. Integrated T-odd distribution functions 
are assumed to be zero and unpolarized integrated T-odd fragmentation functions ap- 
pear only at subleading twist |99]. Single spin asymmetries at leading order in inverse 
powers of the hard scale must therefore arise from intrinsic transverse momenta. 

The validity of the applied theoretical description depends on two issues. The first 
issue is connected to factorization. This was discussed in the previous chapter and in 
the present study it will be assumed to hold. The second issue is related to asymptotic 
freedom. It is expected that only those processes can be described perturbatively in 
which observed external hadrons are well separated in momentum space. The observed 
outgoing hadrons must have a large perpendicular momentum with respect to the beam- 
axis. For observed hadrons close to the beam-axis, one does not only have a problem 
with perturbation theory, but one will also encounter interference effects between the 
hard scattering process and the remnant of the incoming hadrons. In that case fracture 
functions need to be included. 

One of the most studied semi-inclusive cross section in hadron-hadron collisions 
is single hadron production. It is also this process in which the first single spin asym- 
metries in inelastic collisions were observed [20]. Since then, single spin asymmetries 
have been measure d in several processes (see for instance Ref. LI 77- .178- .179- ■18(a 
ll 8 iLll 82II1 83U1 H^lt. and extensive theoretical studi es have been made (see for instance 
Ref.TsT, e1 EiJEll EH EM EM EM EMl)- The main theoretical challenge 
with single hadron production at large transverse momentum is that there is no ob- 
servable which is directly connected to the intrinsic transverse momenta of the partons. 
This makes the extraction of the transverse momentum dependent functions complex. It 
could be that the only possible manner of extracting information is to consider a specific 
form for the transverse momentum dependence of distribution or fragmentation func- 
tions (for example exp[-A:^/M^]). Such a form can lead to problems with the gauge 
link because it generally contains higher order transverse moments. At present only 
the first transverse moment of functions can be easily related to the transverse moment 
of functions appearing in semi-inclusive DIS and electron-positron annihilation. In the 
higher transverse moments new matrix elements are involved which are more difficult 
to relate (see also section l4~2t . Besides, for the higher transverse moments convergence 
becomes an issue as well. 

The two hadron production process, in which the two hadrons belong to different 
jets and are approximately back-to-back in the perpendicular plane, does offer an ob- 
servable directly sensitive to the intrinsic transverse momenta. The fact that the two 
hadrons are not completely back-to-back can at leading order in as be interpreted as an 
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effect from intrinsic transverse momentum. This process will be studied in this chapter. 
Having obtained transverse target-spin asymmetries, a simple extension will be made 
to jet-jet production by just summing over the observed hadrons and observing the jet 
instead. This latter process has been studied by Boer and Vogelsang in Ref. 1 192] (see 
also Vogelsang, Yuan L193il for related work). The study here should be considered as 
an extension of Ref. lll92ll to hadron-hadron production including a full treatment of 
gauge links. 

To find out what observables are present, let us reconsider the Drell-Yan process. 
The extraction of the intrinsic transverse momenta is connected with momentum con- 
servation in the hard scattering cross section, which in the case of Drell-Yan is ex- 
pressed by a four dimensional delta-function in Eg. 13.761 The presence of a hard scale, 
originating from the electromagnetic interaction involving two hadrons, allows for a 
Sudakov-decomposition. This enables one to eliminate two delta-functions, leading to 
the fixed light-cone momentum fractions xi and X2- The remaining two-dimensional 
delta-function, d^ipT+kr-qr) in Eg. 13.801 is directly connected to the intrinsic trans- 
verse momenta of quarks. Note that the momenta pr and kr are transverse with respect 
to their parent hadron. These momenta can be accessed by considering azimuthal asym- 
metries (see Drell-Yan in appendix B.At . 

In the hadron-hadron production process, there are four correlators and two hard 
scales present. If we use the large momentum difference of the initial hadrons to fix the 
perpendicular plane (±) and the light-cone momentum fractions of the initial quarks, 
then the remaining two dimensional delta-function reads 6^(pi±+p2±-ki±-k2±). The 
momentum pi± is already transverse with respect to its parent hadron, so pi±=p\T±- 
This is not the case for A;,^ because the momenta of the outgoing hadrons, Ki, have large 
perpendicular components. It is convenient to make the decomposition kj-Ki I Zi+kn , 
where kjj is defined to be transverse with respect to Kj and is of a hadronic scale. The 
large momentum difference between K^^ and K2±, which is still present in the two- 
dimensional delta-function, can then be used to fix one of the remaining light-cone 
momentum fractions, z,-. This leaves one delta-function behind which contains the non- 
back-to-backness and forms a natural observable as we will see. In the case of two-jet 
production there is no other light-cone momenta to fix. In that case the sum of the two 
jet momenta in the perpendicular plane is akeady proportional to the intrinsic transverse 
momenta of the partons, providing one a two-dimensional vector variable. 



5.2 Calculating cross sections for hadronic scattering 

The calculation for the cross sections will be outlined. After considering kinematics, 
observables will be defined which are sensitive to transverse momentum dependent 
functions. The cross sections will be expressed in an elementary hard scattering cross 
section, fragmentation functions, and distribution functions. Those functions are de- 
fined through bilocal matrix elements and contain a gauge link which depends on the 
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squared amplitude diagram (which we also call subprocess). A simplification of the 
cross section will be achieved by introducing gluonic-pole cross sections. This enables 
one to express the cross section in simple hard scattering cross sections convoluted with 
the in chapter|2|defined distribution and fragmentation functions. 



Kinematics 

For the ease of the calculation we will work in the center of mass frame of the incoming 
hadrons. The incoming hadron with momentum P\ fixes the direction of the z-axis and 
the spatial momenta which are perpendicular to this axis will carry the subscript ±. The 
hadron with momentum f 2 enters the process from the opposite direction. A hard scale 
is set by s = {P\ + Pzf'- The pseudo-rapidity is defined as //,■ = - ln(tan(0,72)), where 
9i is the polar angle of an outgoing hadron with respect to the beam-axis. We introduce 
a scaHng variable x-,^ defined as x,j_ = 2|K,i|/ ^fs, where Ki is the momentum of an 
observed outgoing hadron. These definitions yield the following relations 

Pi-Ki = \sxi^e-'" + 0{M^), Pi-Ki^ jsxi^e'^' + 0{M^), 

Pi-K2^ \sx2^e-'^' + 0(M\ Pi-Ki^ \sx2^e'^' + O(M^). (5.1) 



The Mandelstam variables for the partons are defined as 

s = ipi+p2f, i=ipi-kif, u = (pi-k2f, (5.2) 

and fulfill s + i + u - p\ + p\+ k^^ + k\. The variable y, which is observable and the 
analogue of the y variable used in semi-inclusive DIS, is defined to be 

y=^- , ^ , A l+0{M^ls)). (5.3) 

s exp(77i - 772) + 1 ^ ' 



Defining observables 

Using the diagrammatic expansion (see also chapter|2j the cross section for two hadron 
production reads (0,'s are the azimuthal angles of the observed hadrons) 

1 d'K, d^K2 2 

dcr = J[ 

2s (InyiEK, (27t)^2Ek, 

xij_X2i_s d(pi d02 1 I ^IM^W 

- -Axi^dx2^dmdj]2^^JlUl+0{ — , (5.4) 



128(27r)4 '-" ^-^ " '^ 27T 27T \ \ s 
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where ^[^ = j;^ j'^^l^^(2;r)45(Pi+P2-^i-^2-^x) |MP which is generically ex- 
pressed as 

^ = fdVi dVa d'^/ti d'*fe2 (27r)'* 5'*(/?i+/?2-fei-fe2) 

X Tr^'^ { (5(/?i) (8> 0(/?2) ® A(yti) ® A(/t2) 

®H{pi,p2,kuk2)®H*{px,p2,kx,k2)], (5.5) 

suppressing the momenta and spins of parent hadrons. The symbol ® represent a con- 
volution in color and Dirac indices, and the hard elementary scattering amplitudes are 
denoted by H. The gauge links of the correlators can be derived by applying tech- 
niques developed in the previous chapter. Those gauge links depend in general on the 
subprocess. 

Each of the external partons and its parent hadron have both a nearly light-like 
momentum in a more or less common direction (parton model assumption). The intro- 
duction of the light-like vectors, np,, with np^ ■ hp^ - 1 (bar denotes reversal of spatial 
components) such that np. is proportional to P, in the asymptotic limit, allows us to 
classify at which order components appear The components pi ■ Ftp. appear at order 
yfs (and similarly for fragmentation), the spatial transverse components appear all at 
a hadronic scale. The light-like components which are left, pi ■ np., appear at order 
M^/ -\/s. Parton momentum fractions are defined as usual 

Pi ■ np, Ki ■ fiK, 

Xi = _ , Zi = - — ^. (5.6) 

Pi ■ np. ki ■ nK, 

The light-cone momentum components pi-np. can be simply integrated because they are 
suppressed in the hard parts. The hard scale s can be used to fix the incoming light-cone 
momentum fractions, xi and X2, in terms of zi and Z2- This gives 



/ 



dVi dVa d'^ki d^k2 6^{pi + p2 - ki - ^2) 



^1 



d(pi ■ npj d(p2 ■ np^) d(fei ■ n/f,) A{k2 ■ nK^) 6{xi - \r ■ P2) 6(x2 - -r ■ Pi) 



X I d^piT d^P2T d^kiT Sk2T 6^{qT^ - rj_) 



X fdipi ■ np,) d(p2 ■ np,) d(ki ■ nK,) d{k2 ■ nK.} (l + 0{M^ I s)) , (5.7) 

with qp = PIT + PIT - k\T - k2T, r = 1 . (5.8) 

Zl Z2 

The transverse parton momenta ({piT, P2t, ^ir> ^2r)) are four-vectors and defined to 
be transverse with respect to their parent hadron while the symbol ± means perpendic- 
ular with respect to Pi and P2- The vector qj is thus of a hadronic scale and all its 
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components are in general nonzero. In the case of jet-jet production the integrals over 
ki and k2 do not appear in the expression above, qr = Pit + Pit, and Kt/zi is replaced 
by 1^'^ . In that case the vector r± gives access to the intrinsic transverse momenta of the 
initial partons (^j-) which is similar to Drell-Yan. 

Another hard scale is formed by the scalar product of Ki and K2. This scale, which 
is present in S^iqTj^ - r^), can be used to express one of the light-cone momentum 
fractions in terms of the others, so ziixi,X2,Z2) or Z2ixi,X2,zi)- One remains with an 
integral over one momentum fraction. Since the choice for this fraction is arbitrary, a 
more symmetrical expression can be obtained by introducing an additional integral over 
x±, giving (for details see Ref. 14(JR ) 



/ 



d>i (fp2 (Tki (fk2 d'^ipi +p2-k\- ^2) 
4 



-■- I dx_L 

S-^Xi_lX2_l J 



d(pi ■ np,)d(p2 ■ np,) d(ki ■ rixt) d(k2 ■ rix,) 



6[x 



xb jci 



2-1 j_ 



(e'" +e'^'-)) 6{x2- ^x^{e-'" +e-'^'-)) 



xS 



-1 Xj_ 



— 6[z-2- — 

Xl± \ X2± 



/2 2 2 2 / ^T± ' ^IN 1 
d Pit d p2T d k\T d ^27 ^ p 2^1. sine 

: ^ d{pi ■ np,) d{p2 ■ np^) d{k, ■ "if,) d{k2 ■ hk.} (l + (9(M/ V^)) , (5.9) 



e2± 




e2N 



where e,^ = TiTu/IKul, 4^ = {-2 1 s)e^''P'' P ,^ :^ 

xPiv^i^p and 64> is defined in Fig. 15. II and is eiiv 

of order 0{M/ s/s). The above expression il- 
lustrates that the non-back-to-backness, 6(p, pro- eu 
vides access to the transverse momenta of the 
quarks (via qT±)- Weighting the cross section 
with sin 6(p produces a projection of qT which 
leads to the first transverse moments of distribu- 
tion and fragmentation functions. 

The following cross sections ai'e now defined Figure 5.1: Plane perpendicular to 

the incoming hadron momenta. 

r da- 
{d(r}= d<f,2—- 
J d(p2 

d.xij^ dx2j_ df/i d772 d0i 



32ns 



2n 



rdxj^ 

J X_L 



^(x^). 



(5.10) 
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{^sm(6cf>)dcr)= fd^z {sm(6(f>)^ 

dxij_ dx2j. drii di]2 d0i f dxj. ^ , . 

= 3W7^ ^J ^-i--^^^--)' 

(5.11) 

where 

Y(xJ = (Prhp,)(P2-np,)(Ki-nK,)(K2-hK,) 

<S,A(z2\k2T)<S>H{pi,p2,ki,k2)<S>H'(pup2,ki,k2), (5.12) 

2::^(^J = (Pr«P,)(P2-«F,)(^rnjf,)(^2-%,) 

X I A^piT d}p2T A^kiT A^k2T [qj] ^(xupir) <Si(S>(x2,P2t) <Si A(Zi\kiT) 

®A{z2\k2T)®H{pi,p2,ki,k2)®H*{pi,p2,ki,k2). (5.13) 

In the expressions the momentum fractions, x\, X2, Zi, and Z2, are a function of x^ via 
the arguments of the first four delta-functions in Eg. 15 .91 

Calculating cross sections 

In the case studied here, gluon distribution and fragmentation functions are neglected. 
Taking only the quark and antiquark correlators into account, the forms of the expres- 
sions for^, 2 and Tg will be presented. These were obtained in Ref. ll40ll by using 
FORM (U^. 

The following sum of diagrams contribute to ^^, E and E^j 

DfJ + 4';] + D^ + Dlf + Dt|l + D^^J + (K, - K2) + (q - q), (5.14) 

where the D stands for the diagrams as displayed in Fig. 15.21 For ^^ the following 
convolution can be obtained 

y[2 ~ I dx^ d^piT d^p2T d^kiT d%7- ^ p^^ _f'^ - ix_L sin 50 1 



dcr' 



[D] 



J] ^ fr(x„pl,) ff\x2,plr)^Df\zuz\k\,) Df\z2,zlk\,X (5.15) 



df 



where the /|.'^' and D[^' represent some distribution and fragmentation functions having 
a diagram-dependent gauge link, the light-cone momentum fractions are fixed by the 
first four delta-functions in Eq. 15.91 and dcr. .^j /df is an elementary parton scattering 
subprocess (or diagram) convoluting the functions. 
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Figure 5.2: Contributions in hadron-hadron production. In the figures only the hard part 
is shown, omitting the correlators and parent hadrons. 



In the unweighted cross section, ( dcr ), the intrinsic transverse momenta can be 
neglected at leading order in Mj yfs. In that case all the gauge links in the correlators 
are the same and on the light-cone, see Fig. 12. lib . The sum over diagrams can be 
performed, yielding a simple parton scattering cross section in E 

Also here the light-cone momentum fractions are a function of x± and dcr,yi; /df is the 
elementary parton scattering cross section. 

In a weighted single spin asymmetry each contribution consists of three integrated 
correlators (T-even) containing gauge links on the light-cone and one correlator (T-odd) 
of which a transverse moment is taken (as discussed in the introduction). The first 
transverse moment of a function can be related to the functions defined in chapter |2l 
(see also the previous chapter, section l4!2t . The first transverse moment of a fragmen- 
tation function is in general a combination of the ones appearing in semi-inclusive DIS 
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and electron positron annihilation. Containing an arbitrary gauge link, a fragmentation 
function can be written as 

dJ^'^*'* = aX') 5f ' + D<'\ (5.17) 

where C is some calculable constant and where definitions of Ref. ll35ll were applied 

Df\z) = 1 (z)W<'Vz) + dH<''(z)) , Df\z) ^ \ (dW('>(z) - dY'''\z)) . (5.18) 

The functions D \z) are independent of the gauge link direction, while the functions 
5. '(z) exist by the presence of the gauge link. The tilde functions therefore vanish 
if the transverse momentum dependent fragmentation functions appearing in electron- 
positron annihilation are the same as in semi-inclusive DIS (for a discussion on univer- 
sality, see subsection l3.5> . For single spin asymmetries this leads for S" to 



^ ^i- D i,j,k,l 



dcr'^] 

L 

df 



- I ^ Z ^■'''(-'> /^<-2) ZC(i) -ir- ^^^^'^ ^'fe^ + ■ ■ ■ ' ^^■i'^) 

^ -^-L i,j,k,l D 

gluonic-pole cross section 

where terms related to the first transverse moment of the T-odd function /j \x\) 
have been explicitly shown. The combination YjDC^c\,\Acr^l^li /df is given the name 
gluonic-pole scattering. The color factors, C\j , are obtained by comparing the diagram 
dependent functions, / , with the in chapter|21defined functions /. . For the frag- 
mentation functions D. , which are similarly defined as T-odd distribution functions 
(compare Eg. 15. 181 with Eq. l2.63> . one also finds color factors. For the functions D ' 
all the color factors are simply 1 because those functions are independent of the hnk 
direction. The functions D are therefore convoluted with ordinary parton scattering 
cross sections. 

It may be good to point out that since the first transverse moment of a T-odd distri- 
bution function is of order g and higher, the calculation of a single gluon insertion in the 
diagrammatic approach (for the gauge link) is already sufficient to obtain the color fac- 
tors C\-. . The fact that the all order insertions provide fuUy gauge invariant correlators 
can be seen as a consistency check. 



5.3 Results for cross sections and asymmetries 

In this section an explicit cross section and several single transverse target-spin asym- 
metries will be given. The unpolarized observed hadrons, or jets, are assumed to be 
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approximately opposite in the perpendicular plane. The results are based on Ref. ^40 
in which explicit expressions as a function of 3' can also be found. The quark-quark 
subprocesses together with the gluonic-pole subprocesses are given appendix IS.AI 



Unpolarized cross section for hadron-hadron production 



The unpolarized unweighted cross section for hadron-hadron production is found to 
be a simple convolution between the integrated functions and an elementary parton 
scattering cross section, reading 



( dcr ) — djcij, dx2± drji dri2 



dip I r dxx 

In J Xi 



2 /rui)/rfe)^^^2if^^z)f(zi)z)rfo) 



X 

?1'?2?3?4 



xh+Ol — \+0(as)\, (5.20) 
integrated cross section for hadron-hadron production 



m 



where the summation is over all quark and antiquark flavors. The arguments of the 
functions are defined through xi - Xx(exp(77i)H-exp(772)), X2 - Xi(exp(-77i)-i-exp(-/72)), 
and Zi - x,j./xj.. Furthermore, the variable s depends on x± and y through the relation 
s^xls/(4y(l-y)). 



Single transverse target-spin asymmetries 



In the sin 6(p weighted asymmetry for hadron-hadron production, gluonic-pole scatter- 
ing subprocesses are encountered. The asymmetry reads 
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( i sin((J0)dcr > 

d0i f dxj^ 

= dxij. dx2j. drii dri2 -^— cos(0i - tps) 

2n J x^ 



{ 7^ X /^■^}^^(xi)/f (X2)^ ^"^''^g"^-^- gf (zODffe) 



qiqiqjqA 



-^ > qxqiq^qA 

Xj_ ^js ■^—' 2 dt 

Xj. Vs ■'^-' 2 df 

x|l+<9(^j + 0(«s)j, (5.21) 

weighted asymmetry for hadron-hadron production 



where the summation is over all quark and antiquark flavors. The arguments of the func- 
tions are defined through xi - Xj^(exp77i + exp?;2), X2 - Xi(exp-77i + exp-?72), and 
Zi - Xi^/xj_. Further, s depends on x± and y through s = x^^s/{4y{l - y)). 



The asymmetry for hadron-jet production where the hadron belongs to a dififerent 
jet is obtained by replacing Di(z2) by ^(1 - Z2)^J2q which fixes x^ to be X2j. (the sub- 
script 2 refers to the measured jet), and setting all other functions to zero. This yields 
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< i sin(50)da- ) = dx^ dx2i drji Arji -— !- cos(0i - cps) 



(^ X K^\-^)fr^-2)' {'''''2r'' Dfizo 



X2± yfs 

Mo V-l ... , >.„. irn, > ^ 



dAo- t- 



Ts^T 



-^2^ "vAs ■'^ 2 df 

^-^ ' qiqiqyqt 

X2± V* 2 df 

^-L V qiqiqiqt 



. „ dA(T , 






:|l+(9l^j + <9(ffs)j. (5.22) 

weighted asymmetry for hadron-jet production 



The sin 50 weighted asymmetry for jet-jet production, where the jets are approxi- 
mately back-to-back in the perpendicular plane, reads 



27T 

M] \— 1 „„. inu , „,7, , , * dcr^j^2->iji3ij4 



< 2 sin((J0)dcr ) = dxij_ dx2x drji dT]2 -^^ cos(0i - ^5) 5(xij.-X2i) 



X 



{^ Z r^^"(-.)/rfe)| 



Mo vn un, , .? 'l'^'^<?Ii?;^93 94 



' qiqiqiqt 

ll+ol^] + 0{as)\. (5.23) 

weighted asymmetry for jet-jet production 



As discussed in the previous section, in jet-jet production effects from intrinsic trans- 
verse momentum can also be studied through the vector r". The asymmetry given here 
is just one projection of this vector. An r^-weighted asymmetry could give a better 
separation between the two contributions. 
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5.4 Summary and conclusions 

Effects from intrinsic transverse momentum were studied through single spin asymme- 
tries in hadron-hadron collisions. Assuming factorization, the diagrammatic expansion 
was applied to derive the tree-level expressions and to illustrate effects of the gauge 
link. That this effect can be more than just a sign was pointed out in chapter0]and is a 
generalization of the earlier work in chapter|3 Contributions to the cross sections were 
considered from quark and antiquark distribution and fragmentation functions. For the 
contributions from gluons a similar procedure can be followed. 

A challenge was how to extract the transverse momentum dependent functions from 
the cross section. It was found that by weighting with a specific angle (the non-back- 
to-backness), the asymmetry becomes proportional to three integrated distribution and 
fragmentation functions and one distribution or fragmentation function of which the 
first transverse moment is taken. This method has the advantage that no momentum 
dependence of functions needs to be assumed, offering a model-independent way of 
studying effects from intrinsic transverse momentum. Whether such a method can also 
be applied in the case of single hadron production remains to be seen. 

The distribution and fragmentation functions contain gauge links which were cal- 
culated by using the prescription of the previous chapter Those gauge links become 
relevant when considering transverse moments of the transverse momentum dependent 
functions. The first transverse moments of the correlators in hadron-hadron scattering 
also appear (with factors) in semi-inclusive DIS and electron-positron annihilation. For 
the higher moments such relations have not been achieved (see chapter|4j. 

The single spin asymmetries were obtained in terms of functions containing a subprocess- 
dependent (diagram-dependent) gauge link. Consequently, the asymmetries were rewrit- 
ten as a folding of the functions defined in chapter |2l with newly defined gluonic- 
pole cross sections. Those gluonic-pole cross sections are just partonic scattering dia- 
grams where each diagram is weighted with an additional factor (compare for instance 
Eg. 15. 241 with Ea. l5.25t . This factor comes from taking the first transverse moment of 



the T-odd functions. Except for the fragmentation functions D. ^ where such factors do 



not appear, the T-odd functions are of order g and higher (see for instance Eq. I2.71> . 
Therefore, these factors can already be calculated by doing a single gluon insertion. 
The all order insertions, as performed in the previous chapter to obtain gauge invariant 
results, are not necessary to produce these factors but of course provide confidence on 
the approach we have followed. Since single gluon insertions are sufficient to produce 
these factors, comparisons with the results of Qiu and Sterman |57] are possible. 
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5.A Partonic cross sections 

In this appendix the elementary and gluonic-pole scattering subp roces ses will be listed. 
The elementary expressions are taken from Bacchetta, Radici |156] in which minus 
signs from interchanging fermions are made explicit. Interchanging ki with fe is equiv- 
alent to interchanging u with f or y <-> (1 - y). For constructing the cross sections a 
useful relation is 4ji:iX23'(1 - y) = x]_[l + 0{M^/s)] which can be employed underneath 
the x^ integral. 

Quark- quark scattering 

The unpolarized quark-quark scattering subprocesses are given by 



d<L,,' 4nals^ + u2 




d&r^ , 

qq'^q'q 

dt 


Anaj s^ + f 


di 9.?2 p ' 


9*2 fi2 


d.^!,:;'l,, 4nal f 








di 21 s^ iu 




giving for the total cross section 








d&gq-.qq _ ^^^qq' ^qq' ^ ^^m'^q'q 


-2 


^" qq^qq 





dt dt dt dt 



(5.24) 



For unpolarized gluonic-pole-quark scattering (gluonic pole is associated with first 
quark) the above expression is modified into 

dCTonn^rin _ ^[tt,qq] qq' ^qq' ^[uil.qq] qq'^q'q _ r.^[tu,qq] '^'^qq^qq ,, ^,. 



dt ^-/c' di ^'/(i* di G,/(i) ^^ 

cm '■i^^'^cjm 

The relevant polarized quark-quark subprocesses are 

q^'^qTq' l^nO^ US ^^"^q^^q^q f^^s ^ 



Consulting table 14. 21 one has the following factors: C[l'f'|^ = 1/2, C[i'"f|'', = 1/2, and 

^[tu.qq] _ _T ,r. 
C,/(.l) ~ ^'^ 

' e relevant p 

dAa-[';l, , , 8^Q,2 .. dAo-l"" , 2,na^ 

q'q'—>q'q' _ OJiU^ US q'q—>q'q _ OJiU^ 

df "~ 9^2 F"' ^f ~~27s2 f 

giving for the cross sections 

dAo-T T dAo-";!, ,, dAa-['"l 



'qTq^qTq qTq'^q^' q^^q^ 



dt dt dt 

dA&qTq^^qq _'^^^q1qT^qq _ ^StTuI 

df " df ~ ~ 27 ^2 ' 



(5.26) 
(5.27) 



5. A Partonic cross sections 141 



For the polarized gluonic-pole scattering cross sections, the above expressions are mod- 
ified into 

dAo-T -T dAo-f"', ,, dAo-f'"' , 

""" I'l^gl'l _ p{tt.qq] q'q'^q'q' _ ^[ni,qq] q'q^q'q ^t- ^(,s 

di ~^G,D(I) ^f ^G,D(1) ^f ' ^^■^'^> 

If"] 



qTgq'-iqq _ „[m,qq] q^q^^qq ^^ ^g. 

df ~ '^•■^(2) di ' 

TableOgives the relations: Cf'/,f^^ = -1/2, Cl"'^^^^ = 3/2, C^J^f^ = -3/2. 



Quark-antiquark scattering 

The unpolarized quark-antiquark subprocesses are given by 

do-"'! -, 47rQ'2?2 + ,o2 do-'"' ,-, Ana^ f + f? 

qq'^qq -r/iix^ A -t- W ^ qq^q'q' ^niz^. t -f- u 

dF "912 P~' df "912 ^2 ' 

d^g^^47ra|«2 

df 27^2 ts' 

giving for the cross section 

H" - - dM"^ dfx'^'' dix'"' 

aO'qq^qq _ qq'^qq' qq^q'q' _ q W^W /^ oq-, 

df df df df 

For the unpolarized gluonic-pole scattering cross section, the above expression is mod- 
ified into 

H " - do-'"' do-'*'' do-'"' 

OO'gqq-tqq _ „[n,qq] qq'^qq' ^[ss.qq] qq^q'q' _ r.^[ts,qq] qq-'qq rs "^ n 

df ~^c,f(i) df +^c,/(i) df G,/(i) ^f ' ^^■^'> 

where 4%', = -3/4, C^^ = 5/4, c[:';;^f) = 5/4 (using tablegS. 
The polarized quark-antiquark scattering subprocesses are 

dAo-'"' ,_, ^„a^ tf. dAo-'7' _ ^„a^ .-> 

q'q'^q'q' OJlu.^ m q' q' —>qq OJlu^ u 

df ~~ 9^2 i2' Jf "~27S2 I' 

df ~ 9*2 F"' ^f - 27*2 7' 
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giving for the cross sections 

dAa-[-'/\ , , dAo-tfl, 

(5.32) 

(5.33) 

(5.34) 
df df 27 .?2 

For the gluonic-pole scattering cross sections the above expression is modified into 

df G,/(2) df G,/(2) df ' ' 

dAo-T- -T- dAo-''"_, ,., dAo-"/^ ,_ 

""" I'g^gg'g _ ^[",??] q'q'-tq'q' _ „lts,qq] q'tj^q'q ,, ~ ,, 
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df 
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^qTq 


q^q'^q^' 


q^^q^q 


df 




df 


At 


dA&^T ij_ 


^cfq 




8nal 



df G,D(1) ^f -G,D(1) jj- 

dAo- t- --, dAa-'"l _, 

df ~ "^'^('^ df ' 

where r^'"'*^^ - 5/4 r'"''''^ - 5/4 r^"'*''^ - 3/4 r^"'''^^ - -5/4 r'"'*^^ - -5/4 
(using tablef 
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Summary and 
conclusions 




Effects from intrinsic transverse momentum of partons were studied in several hard 
scattering processes with an emphasis on color gauge invariance. In order to describe 
the processes, the diagrammatic expansion was em ploye d which is a field-theoretical 
approach. Extending the work of Boer, Mulders fll4T and Belitsky, Ji, Yuan 12611 . 
factorization of effects from intrinsic transverse momentum was assumed, and by con- 
sidering an infinite number of diagrams the tree-level expressions including MjQ cor- 
rections were evaluated in chapter|3lfor semi-inclusive DIS, Drell-Yan, and electron- 
positron annihilation. In those processes transverse momentum dependent distribution 
and fragmentation functions were encountered which are defined through matrix el- 
ements in which bilocal operators are folded with a gauge link. From a theoretical 
point of view, the presence of this gauge link (also called Wilson line) is pleasant be- 
cause it makes the bilocal operator invariant under color gauge transformations. These 
gauge links are of increasing interest because by several papers, among which Brodsky, 
Hwang, Schmidt 124], Collins [25], and Belitsky, Ji, Yuan |26], it has been shown that 
the gauge links are not the same in every process and could lead to observable efi'ects. 
Distribution and fragmentation functions containing a gauge link were discussed in 
chapter 12 These functions describe the way in which quarks are distributed in a nu- 
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cleon or how a quark decays into a jet and a particular hadron. They form a bridge 
between theoretical predictions and experimental observations, and are vital for our 
understanding of the nucleon's substructure. In the diagrammatic approach these func- 
tions naturally appear when studying hard scattering processes. The exact form of these 
functions, including the path of the gauge link, is thus not a starting point but rather de- 
rived. 

In the transverse momentum integrated distribution functions the presence of the 
gauge link did not produce new effects in contrast to the transverse momentum depen- 
dent functions. In the latter the gauge link does not run along a straight path between 
the two quark-fields and therefore allows for the existence of T-odd distribution func- 
tions. The existence of such functions was conjectured by Sivers 12 iL 12211 in order to 
explain the observation of single spin asymmetries. As pointed out by Collins |25], 
those functions have the interesting property that they appear with different signs when 
comparing Drell-Yan with semi-inclusive DIS. This interesting prediction should of 
course be experimentally verified. 

Another possible source for T-odd effects was uncovered by Qiu and Sterman 1 84l 
18511 . They suggested that the presence of gluonic pole matrix elements could produce 
single spin asymmetries in hadron-hadron collisions. It was shown in chapter |2] that 
those matrix elements have the same form as matrix elements from which the T-odd 
distribution functions are defined. The presence of the gauge link and the gluonic pole 
matrix elements are therefore in essence the same mechanism. 

The path of the gauge link in the transverse momentum dependent functions runs 
in a particular direction via the light-cone boundary (see Fig. 12. lit . In chapter|3it was 
discussed that this direction is a potential source for new effects or functions (see for re- 
lated work Goeke, Metz, Schlegel 17911 1. While calculating the longitudinal target-spin 
and beam-spin asymmetries for semi-inclusive DIS in chapter|3] it was found that one 
of these new functions, g-'-, could produce a nonzero azimuthal single spin asymmetry 
for jet-production in lepton-hadron scattering. Given the reason for the existence of this 
function, this prediction deserves experimental verification. 

The presence of gauge links in transverse momentum dependent fragmentation 
functions confronted us in chapter|3with issues related to universality. It was found that 
transverse momentum dependent fragmentation functions could appear to be different 
when comparing semi-inclusive DIS with electron-positron annihilation. The reason 
for this difference is that for fragmentation functions there are two possible sources for 
T-odd effects: the gauge link and final-state interactions. If one of the two mechanisms 
is suppressed, relations between the two processes can be drawn. It should be pointed 
out that Collins and Metz obtained in Ref. 1 100, 101 1 in a quite general treatment frag- 
mentation functions which appear with the same sign in the two different processes. 
This interesting result was discussed in chapter|3]and deserves further attention. Also 
from the experimental side this universality issue can be addressed. For instance, by 
comparing properties, like the z-dependence, of transverse momentum dependent frag- 
mentation functions between different processes. 
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In chapter 0] it was discovered that gauge links can arise which are much more 
complex than the gauge links in the discussed electromagnetic processes. Besides be- 
ing more complex, the path of the gauge link turned out to not only depend on the 
process, but also on the subprocess (or squared amplitude diagram). These complex 
structures predict that T-odd distribution functions can in general differ by more than 
just a sign when comparing different processes. In fact only their first transverse mo- 
ment can be straightforwardly compared between different processes. A prescription 
was given in order to deduce the gauge link in squared amplitude diagrams, and it was 
shown in a two gluon-production process that nonphysical polarizations of the gluons 
are canceled among diagrams although separate diagrams are convoluted with functions 
having diagram-dependent gauge links. This is a firm consistency test of the applied 
approach. 

An illustration of these more complex gauge links was given in hadron-hadron pro- 
duction in hadron-hadron scattering. In hadron-hadron scattering the challenge is to 
extract the transverse momentum dependent functions. In chapter|5]such an observable 
was presented. This observable has the advantage that it does not require input on the 
explicit form for the momentum dependence of the functions. By using this observable, 
it was shown that effects from the gauge link yield more than just a sign in single spin 
asymmetries. The expression for the single spin asymmetry was found to be a set of 
elementary scattering subprocesses convoluted with universal integrated functions and 
a function of which a transverse moment was taken. In this latter function the gauge 
link depends on the subprocess. Since the first transverse moment of functions can be 
related to a set of "standard" functions, which also appear in semi-inclusive DIS and 
electron-positron annihilation (just a factor), the asymmetries were rewritten in terms 
of these "standard" functions folded with the newly defined gluonic-pole cross sec- 
tions. Gluonic-pole cross sections are just elementary parton scattering cross sections 
in which the various subprocesses (squared amplitude diagrams) are weighted with a 
particular factor. Contributions from gluon distribution and fragmentation functions 
were discarded for the ease of the calculation, but should be included in future studies 
to make realistic estimates of single spin asymmetries. The definition of transverse mo- 
mentum dependent gluon distribution functions was addressed in chapter|3 but needs 
further improvement. Just like experiments could verify in the Drell-Yan process the 
sign change of T-odd distribution functions, experiments should also be able to check 
the more involved appearance of T-odd distribution functions in hadron-hadron produc- 
tion. Other processes which also contain such more complicated effects are photon-jet 
production, and two jet-production. 

Besides leading to new effects, the gauge link also poses theoretical challenges 
among which the issue of factorization of intrinsic transverse momentum dependent 
effects. This issue underlies most treatments and results in this field incl uding this 
thesis. Recently, significant progress was made by Ji, Ma, and Yuan in Ref. I175[ll76ll 
and by Collins and Metz in Ref. 1 101]. They considered semi-inclusive DIS and Drell- 
Yan at small measured qj {q is the momentum of the virtual photon and q^ M^). 
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One particular issue which was not explicitly addressed is the path of the gauge link in 
distribution and fragmentation functions connected to higher order diagrams (two loops 
or higher). At those orders it was argued in chapterl^that the effect of gauge links might 
endanger factorization of intrinsic transverse momentum dependent effects in semi- 
inclusive DIS and Drell-Yan. This could form a problem for azimuthal asymmetries 
at small qj at leading order in MIQ and at subleading order for fully ^r-integrated 
cross sections. This issue awaits further clarification. Model calculations could be 
very useful here. Just like a model calculation uncovered the possible existence of T- 
odd distribution functions, a two loop model calculation could illustrate some of the 
important aspects related to factorization. 

In this thesis effects from intrinsic transverse momentum were studied in hard scat- 
tering processes. The existing theoretical formalism was further developed and several 
issues have been clarified. Various cross sections and asymmetries were obtained. Their 
measurement can contribute to our understanding of the nucleon's substructure within 
the framework of QCD and very likely will guide physicists in answering several re- 
maining questions. 
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Samenvatting 

Enkelvoudige spin asymmetrieen en ijkinvariantie 
in harde verstrooiingsprocessen 

De materie om ons heen is opgebouwd uit kleinere bouwsteentjes. De kleinste bouw- 
steentjes die wij kennen, worden de elementaire deeltjes genoemd en zijn op dit moment 
dus niet op te delen in nog kleinere deeltjes. Voorbeelden van deze elementaire deeltjes 
zijn het elektron en het foton. Op dit moment is de grootte van deze deeltjes niet bek- 
end, ze worden beschouwd oneindig klein te zijn oftewel puntachtig klein. In de studie 
van elementaire deeltjes speelt een bepaalde familie van deeltjes, genaamd hadronen, 
een belangrijke rol. In tegenstelling tot de elementaire deeltjes hebben hadronen wel 
degelijk een zekere grootte. Een voorbeeld van een hadron is bijvoorbeeld het pro- 
ton met een afmeting van ongeveer 0, 0000000000001 cm. Hadronen zijn echter geen 
elementaire deeltjes, ze zijn opgebouwd uit kleinere puntdeeltjes die quarks en glu- 
onen heten. Op dit moment weten we al aardig wat over de opbouw van hadronen 
uit deze quarks en gluonen, en daar zouden we graag nog veel meer over willen weten. 
Niet alleen vanwege de structuur, maar ook omdat de daaraan gerelateerde natuurkunde 
uitdagend en verrassend is. In dit proefschrift zullen dan ook de drie pijlers van de mod- 
erne natuurkunde: de relativiteitstheorie, de quantumtheorie en het concept symmetrie 
(in het bijzonder ijkinvariantie), worden aangewend om een beschrijving te realiseren. 

De natuurkunde die nodig is voor de beschrijving van elementaire deeltjes is in 
veel opzichten totaal anders dan waar wij in ons dagelijks leven aan gewend zijn. Zo 
spelen op die hele kleine schaal quantumeffecten een belangrijke rol. Intuititief zijn we 
gewend dat objecten gelokaliseerd zijn, betekenende dat objecten altijd op een bepaalde 
lokatie zijn (de Dom staat bijvoorbeeld ergens in Utrecht). In een quantumtheorie is dit 
anders. De positie van een bepaald deeltje is wel meetbaar, maar indien er geen meting 
wordt uitgevoerd, dan is het deeltje in zekere zin overal (in dat geval is de Dom overal 
in Utrecht). Dat lijkt misschien absurd (in het geval van de Dom), maar voor hele kleine 
deeltjes is dit wel degelijk het geval en experimented geverifieerd. 

Het feit dat deeltjes niet meer gelokaliseerd zijn, leidt tot effecten die in de klas- 
sieke theorie (zonder quantumeffecten) niet bestaan. Een voorbeeld is het twee-spleten- 
experiment met elektronen, zie figuur l6Jl Wanneer een elektronenbron een elektron per 
tijdseenheid produceert, dan blijkt er een interferentiepatroon te ontstaan op het scherm 
achter de twee spleten. Dit interferentiepatroon is te verklaren door aan te nemen dat 
het enkele elektron in plaats van door een spleet, door beide spleten tegelijk is gegaan. 
Het niet lokaal zijn van de deeltjes is dus de oorzaak van het effect. 

Doordat deeltjes niet gelokaliseerd zijn, is de beschrijving van de krachten op deelt- 
jes anders dan in een klassieke theorie. In hetzelfde twee-spleten-experiment veran- 
dert het interferentiepartoon bijvoorbeeld als er een lange staafmagneet wordt geplaatst 
tussen de twee spleten (in M). Ondanks dat het elektron klassiek gezien deze magneet 
niet voelt, verandert het interferentiepatroon toch. Dit quantumeffect is simpelweg on- 
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Figure 6.1: Een voorstelling van een elektron uit B dat een beeld produceert op het 
scherm. Doordat het elektron in een quantumtheorie niet lokaal is, golft het elektron 
als het ware door de ruimte naar het scherm toe (het scherm is bijvoorbeeld een fo- 
tografische plaat). Deze golven zijn aangegeven met halve cirkels en zijn vergelijkbaar 
met watergolven. Doordat deze golven door twee spleten reizen, ontstaan er efFectief 
twee golfbronnen aan de rechterkant van de twee spleten. Aan de rechterkant van de 
twee spleten lopen de golven daarom door elkaar heen en ontstaat er op het scherm 
een patroon. Dit patroon wordt gevormd door de snijpunten van twee golflijnen en 
wordt het interferentiepatroon genoemd. Het interferentiepatroon verandert wanneer er 
in M een lange staafmagneet wordt geplaatst die het vlak in- of uitwijst. Ondanks dat 
deze magneet geen elektrische of magnetische velden produceert die voelbaar zijn voor 
de elektrongolven, verandert het interferentiepatroon toch. Deze verandering is onder 
andere afhankelijk van de richting in welke de magneet wijst: het vlak in of uit. 
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denkbaar in de klassieke natuurkunde en heet het Aharonov-Bohm effect. Volgens de 
klassieke natuurkunde produceert de magneet een magnetisch veld maar in deze speci- 
fieke situatie kunnen de elektrongolven het magnetisch veld echter niet voelen. Om het 
effect te kunnen verklaren moet er dus nog lets anders zijn wat de elektrongolven wel 
beinvloedt. Dit andere veld wordt het ijkveld of fotonv eld genoemd en dit veld blijkt 
inderdaad op het quantumniveau met de elektrongolven te kunnen wisselwerken. Het 
effect hangt samen met een theoretisch begrip genaamd ijkinvariantie. 

Wanneer we de allerkleinste deeltjes bestuderen moeten we onder andere deze ef- 
fecten in kaart brengen. Een probleem bij de studie van de structuur van hadronen is 
dat de bouwelementen, de quarks en gluonen, alleen lijken te bestaan in de hadronen 
zelf. In tegenstelling tot alle andere deeltjes die wij kennen (en dat zijn er heel wat), 
zijn quarks en gluonen nog nooit vrij geobserveerd. Het zijn dus een soort legoblokjes 
die afzonderlijk niet lijken te bestaan. Dit roept natuurlijk vragen op die op dit moment 
nog niet volledig zijn begrepen, en dit maakt de verdeling van quarks en gluonen in een 
hadron natuurlijk des te interessanter 

De structuur van hadronen wordt bestudeerd in experimenten. Een van die ex- 
perimenten is elektron-hadron verstrooiing. In dat experiment botst een elektron met 
hoge sneUieid op een hadron. Door deze botsing breekt het hadron op in stukken en 
dat levert een scala aan andere deeltjes en andere soorten hadronen op, allemaal met 
bepaalde snelheden en richtingen. Het theoretische model wat wordt toegepast is een- 
voudig in oorsprong. Het idee is dat in bepaalde situaties de kans dat een bepaald 
deeltje of hadron na een botsing wordt gemeten evenredig is met de kans dat een quark 
in een hadron wordt geraakt, vermenigvuldigt met de kans dat ditzelfde quark vervalt 
in het gemeten hadron. Er wordt dus verondersteld dat het vervalproces van een quark 
onafhankelijk is van de structuur van het geraakte hadron. Deze aanname, genaamd 
factorizatie, kan gedeeltelijk worden onderbouwd en heeft in veel experimenten ook 
een goede beschrijving gegeven. 

Een observatie die tot voor kort nog onbegrepen was, is het optreden van enkel- 
voudige spin asymmetrieen in elektron-hadron verstrooiing en hadmn-hadron verstrooi- 
ing. In het laatste experiment worden in plaats van een elektron op een hadron te schi- 
eten twee hadronen keihard op elkaar geschoten en wordt er een ander soort hadron na 
de botsing gemeten. Hadronen kunnen in het algemeen een spin hebben. Dat betekent 
dat ze als het ware tollen om hun eigen as. Als in het experiment blijkt dat waneer 
de spin van het ingaande hadron wordt omgedraaid ook de snelheid of richting van het 
gemeten deeltje verandert, spreekt men van een enkelvoudige spin asymmetric. Deze 
asymmetrieen treden ook op in de eerder besproken elektron-hadron verstrooiing. In 
figuur n~3l op bladzijde^Jis voor dit proces de gemeten asymmetric (langs de verticale 
as) weergegeven. 

In dit proefschrift wordt de oorzaak van deze enkelvoudige spin asymmetrieen on- 
derzocht. In hoofdstuk|2lwordt een aantal hoge energie verstrooiingsprocessen geintro- 
duceerd. De basis van de theorie wordt uitgelegd en de distributie en fragmentatie 
functies worden gedefinieerd die later blijken op te duiken in experimentele groothe- 
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den. Deze functies beschrijven de structuur van het hadron en bevatten de informatie 
die we willen weten. In hoofstuk|2is een nieuwe klasse van functies ontdekt die ons 
begrip van de structuur kunnen vergroten. Een voorbeeld van deze nieuwe functies is 
g-'-, aangegeven in vergeliiking l2.65l (oD bladziidel34>. 

In hoofdstuk|3]wordt getracht een aantal verstrooiingsprocessen te beschrijven. In 
deze beschrijving blijkt dat ijkinvariantie een grote rol speelt. Ijkinvariantie is een 
wiskundige symmetric die aanwezig zou moeten zijn in de theoretische beschrijving. 
Tot voor kort was dit niet helemaal netjes meegenomen, maar men wist ook niet of 
het missen van de symmetric werkelijk een probleem zou kunnen zijn. Door in dit 
proefschrift de theorie nauwkeurig uit te werken, wordt er een bijdrage geleverd aan 
het herstellen van ijkinvariantie. Behalve dat dit een verbeterde theoretische beschrijv- 
ing oplevert, blijkt er een subtiel effect te zijn wat een natuurlijke verklaring voor het 
optreden van enkelvoudige spin asymmetrieen geeft. Een verklaring in elektron-hadron 
verstrooiing is de wisselwerking tussen de geraakte quark en het ijkveld igluonveld) 
van het hadron (de gluonen zijn geillustreerd in figuur [l~2l door gekrulde lijnen). Het 
effect is daardoor in een aantal opzichten vergelijkbaar met het Aharonov-Bohm ef- 
fect. In beide gevallen wordt het effect veroorzaakt doordat een deeltje (hetzij elektron 
of quark) wisselwerkt met een ijkveld in gebieden waar de elektrische en magnetische 
velden geen effect hebben. 

Deze effecten in kaart brengend, zien we in hoofdstuk|3ldat de functies van hoofd- 
stuk 121 kunnen worden gemeten in elektron-hadron verstrooiing. Vergelijking 13.641 
is bijvoorbeeld belangrijk om de verscheidene functies te kunnen meten. In hoofd- 
stuk0] wordt de theorie verder ontwikkeld om ook hadron-hadron verstrooiing te kun- 
nen beschrijven. Een van de belangrijkste resultaten, maar technisch van aard, wordt 
geillustreerd in de figuren l4!9l en l4. lOl In dit hoofdstuk wordt verder ook de consisten- 
tie van de gefactorizeerde aanpak besproken: Is het vervalproces van quarks in hadro- 
nen wel onafhankelijk van de structuur van het geraakte hadron? Het antwoord op 
deze vraag blijkt niet eenvoudig. In hoofdstuk |5] worden de ontwikkelde technieken 
van hoofdstuk |4] gebruikt om voorspellingen te doen voor asymmetrieen in hadron- 
hadron verstrooiingsprocessen. Experimenten kunnen deze voorspellingen verifieren 
en daarmee de geldigheid van de ontwikkelde ideeen toetsen. 



